| 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 |