1.1 --- /dev/null Thu Jan 01 00:00:00 1970 +0000
1.2 +++ b/os/ossrv/genericopenlibs/liboil/src/ref/squaresum_f64_ref.c Fri Jun 15 03:10:57 2012 +0200
1.3 @@ -0,0 +1,145 @@
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 +/**
1.42 + * oil_squaresum_f64:
1.43 + * @d:
1.44 + * @s:
1.45 + * @n:
1.46 + *
1.47 + * Sums the square of each element in @s and places the result in
1.48 + * @dest.
1.49 + */
1.50 +OIL_DEFINE_CLASS (squaresum_f64, "double *d, double *s, int n");
1.51 +
1.52 +static void
1.53 +squaresum_f64_ref(double *dest, double *src, int n)
1.54 +{
1.55 + double sum2 = 0;
1.56 + double errsum = 0;
1.57 + double tmp;
1.58 + double x;
1.59 + int i;
1.60 +
1.61 + for(i=0;i<n;i++){
1.62 + tmp = sum2;
1.63 + x = src[i]*src[i];
1.64 + sum2 += x;
1.65 + errsum += (tmp - sum2) + x;
1.66 + }
1.67 +
1.68 + *dest = sum2 + errsum;
1.69 +}
1.70 +OIL_DEFINE_IMPL_REF (squaresum_f64_ref, squaresum_f64);
1.71 +
1.72 +/**
1.73 + * oil_squaresum_shifted_s16:
1.74 + * @d:
1.75 + * @s:
1.76 + * @n:
1.77 + *
1.78 + * Square each element in @s and divide by (1<<15), and sum the
1.79 + * results, placing the final result in @d.
1.80 + */
1.81 +OIL_DEFINE_CLASS (squaresum_shifted_s16, "uint32_t *d, int16_t *s, int n");
1.82 +
1.83 +static void
1.84 +squaresum_shifted_s16_ref(uint32_t *dest, int16_t *src, int n)
1.85 +{
1.86 + uint32_t sum2 = 0;
1.87 + uint32_t x;
1.88 + int i;
1.89 +
1.90 + for(i=0;i<n;i++){
1.91 + x = (src[i]*src[i] + (1<<14)) >> 15;
1.92 + sum2 += x;
1.93 + }
1.94 +
1.95 + *dest = sum2;
1.96 +}
1.97 +OIL_DEFINE_IMPL_REF (squaresum_shifted_s16_ref, squaresum_shifted_s16);
1.98 +
1.99 +
1.100 +
1.101 +#ifdef __SYMBIAN32__
1.102 +
1.103 +OilFunctionClass* __oil_function_class_squaresum_f64() {
1.104 + return &_oil_function_class_squaresum_f64;
1.105 +}
1.106 +#endif
1.107 +
1.108 +#ifdef __SYMBIAN32__
1.109 +
1.110 +OilFunctionClass* __oil_function_class_squaresum_shifted_s16() {
1.111 + return &_oil_function_class_squaresum_shifted_s16;
1.112 +}
1.113 +#endif
1.114 +
1.115 +
1.116 +
1.117 +#ifdef __SYMBIAN32__
1.118 +
1.119 +OilFunctionImpl* __oil_function_impl_squaresum_f64_ref() {
1.120 + return &_oil_function_impl_squaresum_f64_ref;
1.121 +}
1.122 +#endif
1.123 +
1.124 +#ifdef __SYMBIAN32__
1.125 +
1.126 +OilFunctionImpl* __oil_function_impl_squaresum_shifted_s16_ref() {
1.127 + return &_oil_function_impl_squaresum_shifted_s16_ref;
1.128 +}
1.129 +#endif
1.130 +
1.131 +
1.132 +
1.133 +#ifdef __SYMBIAN32__
1.134 +
1.135 +EXPORT_C void** _oil_function_class_ptr_squaresum_f64 () {
1.136 + oil_function_class_ptr_squaresum_f64 = __oil_function_class_squaresum_f64();
1.137 + return &oil_function_class_ptr_squaresum_f64->func;
1.138 + }
1.139 +#endif
1.140 +
1.141 +#ifdef __SYMBIAN32__
1.142 +
1.143 +EXPORT_C void** _oil_function_class_ptr_squaresum_shifted_s16 () {
1.144 + oil_function_class_ptr_squaresum_shifted_s16 = __oil_function_class_squaresum_shifted_s16();
1.145 + return &oil_function_class_ptr_squaresum_shifted_s16->func;
1.146 + }
1.147 +#endif
1.148 +