1*c9945492SAndroid Build Coastguard Worker #include "libm.h"
2*c9945492SAndroid Build Coastguard Worker
3*c9945492SAndroid Build Coastguard Worker #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
sinl(long double x)4*c9945492SAndroid Build Coastguard Worker long double sinl(long double x)
5*c9945492SAndroid Build Coastguard Worker {
6*c9945492SAndroid Build Coastguard Worker return sin(x);
7*c9945492SAndroid Build Coastguard Worker }
8*c9945492SAndroid Build Coastguard Worker #elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384
sinl(long double x)9*c9945492SAndroid Build Coastguard Worker long double sinl(long double x)
10*c9945492SAndroid Build Coastguard Worker {
11*c9945492SAndroid Build Coastguard Worker union ldshape u = {x};
12*c9945492SAndroid Build Coastguard Worker unsigned n;
13*c9945492SAndroid Build Coastguard Worker long double y[2], hi, lo;
14*c9945492SAndroid Build Coastguard Worker
15*c9945492SAndroid Build Coastguard Worker u.i.se &= 0x7fff;
16*c9945492SAndroid Build Coastguard Worker if (u.i.se == 0x7fff)
17*c9945492SAndroid Build Coastguard Worker return x - x;
18*c9945492SAndroid Build Coastguard Worker if (u.f < M_PI_4) {
19*c9945492SAndroid Build Coastguard Worker if (u.i.se < 0x3fff - LDBL_MANT_DIG/2) {
20*c9945492SAndroid Build Coastguard Worker /* raise inexact if x!=0 and underflow if subnormal */
21*c9945492SAndroid Build Coastguard Worker FORCE_EVAL(u.i.se == 0 ? x*0x1p-120f : x+0x1p120f);
22*c9945492SAndroid Build Coastguard Worker return x;
23*c9945492SAndroid Build Coastguard Worker }
24*c9945492SAndroid Build Coastguard Worker return __sinl(x, 0.0, 0);
25*c9945492SAndroid Build Coastguard Worker }
26*c9945492SAndroid Build Coastguard Worker n = __rem_pio2l(x, y);
27*c9945492SAndroid Build Coastguard Worker hi = y[0];
28*c9945492SAndroid Build Coastguard Worker lo = y[1];
29*c9945492SAndroid Build Coastguard Worker switch (n & 3) {
30*c9945492SAndroid Build Coastguard Worker case 0:
31*c9945492SAndroid Build Coastguard Worker return __sinl(hi, lo, 1);
32*c9945492SAndroid Build Coastguard Worker case 1:
33*c9945492SAndroid Build Coastguard Worker return __cosl(hi, lo);
34*c9945492SAndroid Build Coastguard Worker case 2:
35*c9945492SAndroid Build Coastguard Worker return -__sinl(hi, lo, 1);
36*c9945492SAndroid Build Coastguard Worker case 3:
37*c9945492SAndroid Build Coastguard Worker default:
38*c9945492SAndroid Build Coastguard Worker return -__cosl(hi, lo);
39*c9945492SAndroid Build Coastguard Worker }
40*c9945492SAndroid Build Coastguard Worker }
41*c9945492SAndroid Build Coastguard Worker #endif
42