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 sinhl(long double x) | 4 long double sinhl(long double x) { |
5 { | 5 return sinh(x); |
6 » return sinh(x); | |
7 } | 6 } |
8 #elif LDBL_MANT_DIG == 64 && LDBL_MAX_EXP == 16384 | 7 #elif LDBL_MANT_DIG == 64 && LDBL_MAX_EXP == 16384 |
9 long double sinhl(long double x) | 8 long double sinhl(long double x) { |
10 { | 9 union ldshape u = {x}; |
11 » union ldshape u = {x}; | 10 unsigned ex = u.i.se & 0x7fff; |
12 » unsigned ex = u.i.se & 0x7fff; | 11 long double h, t, absx; |
13 » long double h, t, absx; | |
14 | 12 |
15 » h = 0.5; | 13 h = 0.5; |
16 » if (u.i.se & 0x8000) | 14 if (u.i.se & 0x8000) |
17 » » h = -h; | 15 h = -h; |
18 » /* |x| */ | 16 /* |x| */ |
19 » u.i.se = ex; | 17 u.i.se = ex; |
20 » absx = u.f; | 18 absx = u.f; |
21 | 19 |
22 » /* |x| < log(LDBL_MAX) */ | 20 /* |x| < log(LDBL_MAX) */ |
23 » if (ex < 0x3fff+13 || (ex == 0x3fff+13 && u.i.m>>32 < 0xb17217f7)) { | 21 if (ex < 0x3fff + 13 || (ex == 0x3fff + 13 && u.i.m >> 32 < 0xb17217f7)) { |
24 » » t = expm1l(absx); | 22 t = expm1l(absx); |
25 » » if (ex < 0x3fff) { | 23 if (ex < 0x3fff) { |
26 » » » if (ex < 0x3fff-32) | 24 if (ex < 0x3fff - 32) |
27 » » » » return x; | 25 return x; |
28 » » » return h*(2*t - t*t/(1+t)); | 26 return h * (2 * t - t * t / (1 + t)); |
29 » » } | 27 } |
30 » » return h*(t + t/(t+1)); | 28 return h * (t + t / (t + 1)); |
31 » } | 29 } |
32 | 30 |
33 » /* |x| > log(LDBL_MAX) or nan */ | 31 /* |x| > log(LDBL_MAX) or nan */ |
34 » t = expl(0.5*absx); | 32 t = expl(0.5 * absx); |
35 » return h*t*t; | 33 return h * t * t; |
36 } | 34 } |
37 #elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384 | 35 #elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384 |
38 // TODO: broken implementation to make things compile | 36 // TODO: broken implementation to make things compile |
39 long double sinhl(long double x) | 37 long double sinhl(long double x) { |
40 { | 38 return sinh(x); |
41 » return sinh(x); | |
42 } | 39 } |
43 #endif | 40 #endif |
OLD | NEW |