include/vdt/sin.h

0001 /*
0002  * cos.h
0003  * The basic idea is to exploit Pade polynomials.
0004  * A lot of ideas were inspired by the cephes math library (by Stephen L. Moshier
0005  * moshier@na-net.ornl.gov) as well as actual code.
0006  * The Cephes library can be found here:  http://www.netlib.org/cephes/
0007  *
0008  *  Created on: Jun 23, 2012
0009  *      Author: Danilo Piparo, Thomas Hauth, Vincenzo Innocente
0010  */
0011
0012 /*
0013  * VDT is free software: you can redistribute it and/or modify
0014  * it under the terms of the GNU Lesser Public License as published by
0015  * the Free Software Foundation, either version 3 of the License, or
0016  * (at your option) any later version.
0017  *
0018  * This program is distributed in the hope that it will be useful,
0019  * but WITHOUT ANY WARRANTY; without even the implied warranty of
0020  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
0021  * GNU Lesser Public License for more details.
0022  *
0023  * You should have received a copy of the GNU Lesser Public License
0024  * along with this program.  If not, see <http://www.gnu.org/licenses/>.
0025  */
0026
0027 #ifndef SIN_H_
0028 #define SIN_H_
0029
0030 #include "sincos.h"
0031
0032 namespace vdt{
0033
0034 // Sin double precision --------------------------------------------------------
0035
0036 /// Double precision sine: just call sincos.
0037 inline double fast_sin(double x){double s,c;fast_sincos(x,s,c);return s;}
0038
0039 //------------------------------------------------------------------------------
0040
0041 inline float fast_sinf(float x){float s,c;fast_sincosf(x,s,c);return s;}
0042
0043 //------------------------------------------------------------------------------
0044 void sinv(const uint32_t size, double const * __restrict__ iarray, double* __restrict__ oarray);
0045 void fast_sinv(const uint32_t size, double const * __restrict__ iarray, double* __restrict__ oarray);
0046 void sinfv(const uint32_t size, float const * __restrict__ iarray, float* __restrict__ oarray);
0047 void fast_sinfv(const uint32_t size, float const * __restrict__ iarray, float* __restrict__ oarray);
0048
0049
0050 } //vdt namespace
0051
0052 #endif /* SIN_H_ */