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