liboil/ref/sincos_f64.c

   1 /*
   2  * LIBOIL - Library of Optimized Inner Loops
   3  * Copyright (c) 2003,2004 David A. Schleef <ds@schleef.org>
   4  * All rights reserved.
   5  *
   6  * Redistribution and use in source and binary forms, with or without
   7  * modification, are permitted provided that the following conditions
   8  * are met:
   9  * 1. Redistributions of source code must retain the above copyright
  10  *    notice, this list of conditions and the following disclaimer.
  11  * 2. Redistributions in binary form must reproduce the above copyright
  12  *    notice, this list of conditions and the following disclaimer in the
  13  *    documentation and/or other materials provided with the distribution.
  14  *
  15  * THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR
  16  * IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
  17  * WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  18  * ARE DISCLAIMED.  IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT,
  19  * INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
  20  * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
  21  * SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
  22  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT,
  23  * STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING
  24  * IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
  25  * POSSIBILITY OF SUCH DAMAGE.
  26  */
  27
  28 #ifdef HAVE_CONFIG_H
  29 #include "config.h"
  30 #endif
  31
  32 #include <liboil/liboilfunction.h>
  33 #include <math.h>
  34
  35 /**
  36  * oil_sincos_f64:
  37  * @dest1:
  38  * @dest2:
  39  * @n:
  40  * @s1_1:
  41  * @s2_1:
  42  *
  43  * Calculates sin(x) and cos(x) and places the results in @dest1 and
  44  * @dest2 respectively.  Values for x start at @s1_1 and are incremented
  45  * by @s2_1 for each destination element.
  46  */
  47 OIL_DEFINE_CLASS (sincos_f64,
  48     "double *dest1, double *dest2, int n, double *s1_1, double *s2_1");
  49
  50 static void
  51 sincos_f64_ref (double *dest_sin, double *dest_cos, int n, double *offset,
  52     double *interval)
  53 {
  54         int i;
  55         double x;
  56
  57         for(i=0;i<n;i++){
  58           x = *offset + *interval * i;
  59                 dest_sin[i] = sin(x);
  60                 dest_cos[i] = cos(x);
  61         }
  62 }
  63 OIL_DEFINE_IMPL_REF (sincos_f64_ref, sincos_f64);
  64