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