1.1 --- /dev/null Thu Jan 01 00:00:00 1970 +0000
1.2 +++ b/os/ossrv/genericopenlibs/liboil/src/simdpack/sincos_f64.c Fri Jun 15 03:10:57 2012 +0200
1.3 @@ -0,0 +1,88 @@
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/liboilclasses.h>
1.38 +#include <math.h>
1.39 +
1.40 +/* XXX This gets more inaccurate as N increases */
1.41 +static void
1.42 +sincos_f64_i20_fast(double *dest_sin, double *dest_cos, int n, double *offset,
1.43 + double *interval)
1.44 +{
1.45 + int i;
1.46 + double c,s,dc,ds,temp_s,temp_c;
1.47 +
1.48 + dest_sin[0] = s = sin(*offset);
1.49 + dest_cos[0] = c = cos(*offset);
1.50 + ds = sin(*interval);
1.51 + dc = cos(*interval);
1.52 + for(i=1;i<n;i++){
1.53 + temp_s = s*dc + c*ds;
1.54 + temp_c = c*dc - s*ds;
1.55 + dest_sin[i] = s = temp_s;
1.56 + dest_cos[i] = c = temp_c;
1.57 + }
1.58 +}
1.59 +OIL_DEFINE_IMPL (sincos_f64_i20_fast, sincos_f64);
1.60 +
1.61 +#ifdef HAVE_SINCOS
1.62 +static void
1.63 +sincos_f64_sincos (double *dest_sin, double *dest_cos, int n, double *offset,
1.64 + double *interval)
1.65 +{
1.66 + int i;
1.67 +
1.68 + for(i=0;i<n;i++){
1.69 + sincos (*offset + *interval * i, dest_sin + i, dest_cos + i);
1.70 + }
1.71 +}
1.72 +OIL_DEFINE_IMPL (sincos_f64_sincos, sincos_f64);
1.73 +#endif
1.74 +
1.75 +
1.76 +
1.77 +#ifdef __SYMBIAN32__
1.78 +
1.79 +OilFunctionImpl* __oil_function_impl_sincos_f64_i20_fast() {
1.80 + return &_oil_function_impl_sincos_f64_i20_fast;
1.81 +}
1.82 +#endif
1.83 +#if 0
1.84 +#ifdef __SYMBIAN32__
1.85 +
1.86 +OilFunctionImpl* __oil_function_impl_sincos_f64_sincos() {
1.87 + return &_oil_function_impl_sincos_f64_sincos;
1.88 +}
1.89 +#endif
1.90 +#endif
1.91 +