1.1 --- /dev/null Thu Jan 01 00:00:00 1970 +0000
1.2 +++ b/os/ossrv/genericopenlibs/liboil/src/ref/multsum_ref.c Fri Jun 15 03:10:57 2012 +0200
1.3 @@ -0,0 +1,139 @@
1.4 +/*
1.5 + * LIBOIL - Library of Optimized Inner Loops
1.6 + * Copyright (c) 2003,2004 David A. Schleef <ds@schleef.org>
1.7 + * All rights reserved.
1.8 + *
1.9 + * Redistribution and use in source and binary forms, with or without
1.10 + * modification, are permitted provided that the following conditions
1.11 + * are met:
1.12 + * 1. Redistributions of source code must retain the above copyright
1.13 + * notice, this list of conditions and the following disclaimer.
1.14 + * 2. Redistributions in binary form must reproduce the above copyright
1.15 + * notice, this list of conditions and the following disclaimer in the
1.16 + * documentation and/or other materials provided with the distribution.
1.17 + *
1.18 + * THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR
1.19 + * IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
1.20 + * WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
1.21 + * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT,
1.22 + * INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
1.23 + * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
1.24 + * SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
1.25 + * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT,
1.26 + * STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING
1.27 + * IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
1.28 + * POSSIBILITY OF SUCH DAMAGE.
1.29 + */
1.30 +//Portions Copyright (c) 2008-2009 Nokia Corporation and/or its subsidiary(-ies). All rights reserved.
1.31 +
1.32 +#ifdef HAVE_CONFIG_H
1.33 +#include "config.h"
1.34 +#endif
1.35 +
1.36 +#include <liboil/liboilfunction.h>
1.37 +#include "liboil/simdpack/simdpack.h"
1.38 +#include <math.h>
1.39 +
1.40 +
1.41 +#define MULTSUM_DEFINE_REF(type) \
1.42 +static void multsum_ ## type ## _ref( \
1.43 + oil_type_ ## type *dest, \
1.44 + oil_type_ ## type *src1, int sstr1, \
1.45 + oil_type_ ## type *src2, int sstr2, \
1.46 + int n) \
1.47 +{ \
1.48 + int i; \
1.49 + double sum = 0; \
1.50 + double errsum = 0; \
1.51 + for(i=0;i<n;i++){ \
1.52 + oil_type_ ## type x; \
1.53 + oil_type_ ## type tmp; \
1.54 + x = OIL_GET(src1,sstr1*i,oil_type_ ## type) * OIL_GET(src2,sstr2*i,oil_type_ ## type); \
1.55 + tmp = sum; \
1.56 + sum += x; \
1.57 + errsum += (tmp - sum) + x; \
1.58 + } \
1.59 + *dest = sum + errsum; \
1.60 +} \
1.61 +OIL_DEFINE_IMPL_REF (multsum_ ## type ## _ref, multsum_ ## type); \
1.62 +OIL_DEFINE_CLASS (multsum_ ## type, \
1.63 + "oil_type_" #type " *dest, " \
1.64 + "oil_type_" #type " *src1, int sstr1, " \
1.65 + "oil_type_" #type " *src2, int sstr2, " \
1.66 + "int n")
1.67 +
1.68 +/**
1.69 + * oil_multsum_f32:
1.70 + * @dest:
1.71 + * @src1:
1.72 + * @sstr1:
1.73 + * @src2:
1.74 + * @sstr2:
1.75 + * @n:
1.76 + *
1.77 + * Multiplies each element in @src1 and @src2 and sums the results
1.78 + * over the entire array, and places the sum into @dest.
1.79 + */
1.80 +MULTSUM_DEFINE_REF(f32);
1.81 +/**
1.82 + * oil_multsum_f64:
1.83 + * @dest:
1.84 + * @src1:
1.85 + * @sstr1:
1.86 + * @src2:
1.87 + * @sstr2:
1.88 + * @n:
1.89 + *
1.90 + * Multiplies each element in @src1 and @src2 and sums the results
1.91 + * over the entire array, and places the sum into @dest.
1.92 + */
1.93 +MULTSUM_DEFINE_REF(f64);
1.94 +
1.95 +#ifdef __SYMBIAN32__
1.96 +
1.97 +OilFunctionClass* __oil_function_class_multsum_f64() {
1.98 + return &_oil_function_class_multsum_f64;
1.99 +}
1.100 +#endif
1.101 +
1.102 +#ifdef __SYMBIAN32__
1.103 +
1.104 +OilFunctionClass* __oil_function_class_multsum_f32() {
1.105 + return &_oil_function_class_multsum_f32;
1.106 +}
1.107 +#endif
1.108 +
1.109 +
1.110 +
1.111 +#ifdef __SYMBIAN32__
1.112 +
1.113 +OilFunctionImpl* __oil_function_impl_multsum_f32_ref() {
1.114 + return &_oil_function_impl_multsum_f32_ref;
1.115 +}
1.116 +#endif
1.117 +
1.118 +#ifdef __SYMBIAN32__
1.119 +
1.120 +OilFunctionImpl* __oil_function_impl_multsum_f64_ref() {
1.121 + return &_oil_function_impl_multsum_f64_ref;
1.122 +}
1.123 +#endif
1.124 +
1.125 +
1.126 +
1.127 +#ifdef __SYMBIAN32__
1.128 +
1.129 +EXPORT_C void** _oil_function_class_ptr_multsum_f64 () {
1.130 + oil_function_class_ptr_multsum_f64 = __oil_function_class_multsum_f64();
1.131 + return &oil_function_class_ptr_multsum_f64->func;
1.132 + }
1.133 +#endif
1.134 +
1.135 +#ifdef __SYMBIAN32__
1.136 +
1.137 +EXPORT_C void** _oil_function_class_ptr_multsum_f32 () {
1.138 + oil_function_class_ptr_multsum_f32 = __oil_function_class_multsum_f32();
1.139 + return &oil_function_class_ptr_multsum_f32->func;
1.140 + }
1.141 +#endif
1.142 +