1*2fe8fb19SBen Gras /* k_sinf.c -- float version of k_sin.c
2*2fe8fb19SBen Gras * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
3*2fe8fb19SBen Gras */
4*2fe8fb19SBen Gras
5*2fe8fb19SBen Gras /*
6*2fe8fb19SBen Gras * ====================================================
7*2fe8fb19SBen Gras * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
8*2fe8fb19SBen Gras *
9*2fe8fb19SBen Gras * Developed at SunPro, a Sun Microsystems, Inc. business.
10*2fe8fb19SBen Gras * Permission to use, copy, modify, and distribute this
11*2fe8fb19SBen Gras * software is freely granted, provided that this notice
12*2fe8fb19SBen Gras * is preserved.
13*2fe8fb19SBen Gras * ====================================================
14*2fe8fb19SBen Gras */
15*2fe8fb19SBen Gras
16*2fe8fb19SBen Gras #include <sys/cdefs.h>
17*2fe8fb19SBen Gras #if defined(LIBM_SCCS) && !defined(lint)
18*2fe8fb19SBen Gras __RCSID("$NetBSD: k_sinf.c,v 1.7 2002/05/26 22:01:53 wiz Exp $");
19*2fe8fb19SBen Gras #endif
20*2fe8fb19SBen Gras
21*2fe8fb19SBen Gras #include "math.h"
22*2fe8fb19SBen Gras #include "math_private.h"
23*2fe8fb19SBen Gras
24*2fe8fb19SBen Gras static const float
25*2fe8fb19SBen Gras half = 5.0000000000e-01,/* 0x3f000000 */
26*2fe8fb19SBen Gras S1 = -1.6666667163e-01, /* 0xbe2aaaab */
27*2fe8fb19SBen Gras S2 = 8.3333337680e-03, /* 0x3c088889 */
28*2fe8fb19SBen Gras S3 = -1.9841270114e-04, /* 0xb9500d01 */
29*2fe8fb19SBen Gras S4 = 2.7557314297e-06, /* 0x3638ef1b */
30*2fe8fb19SBen Gras S5 = -2.5050759689e-08, /* 0xb2d72f34 */
31*2fe8fb19SBen Gras S6 = 1.5896910177e-10; /* 0x2f2ec9d3 */
32*2fe8fb19SBen Gras
33*2fe8fb19SBen Gras float
__kernel_sinf(float x,float y,int iy)34*2fe8fb19SBen Gras __kernel_sinf(float x, float y, int iy)
35*2fe8fb19SBen Gras {
36*2fe8fb19SBen Gras float z,r,v;
37*2fe8fb19SBen Gras int32_t ix;
38*2fe8fb19SBen Gras GET_FLOAT_WORD(ix,x);
39*2fe8fb19SBen Gras ix &= 0x7fffffff; /* high word of x */
40*2fe8fb19SBen Gras if(ix<0x32000000) /* |x| < 2**-27 */
41*2fe8fb19SBen Gras {if((int)x==0) return x;} /* generate inexact */
42*2fe8fb19SBen Gras z = x*x;
43*2fe8fb19SBen Gras v = z*x;
44*2fe8fb19SBen Gras r = S2+z*(S3+z*(S4+z*(S5+z*S6)));
45*2fe8fb19SBen Gras if(iy==0) return x+v*(S1+z*r);
46*2fe8fb19SBen Gras else return x-((z*(half*y-v*r)-y)-v*S1);
47*2fe8fb19SBen Gras }
48