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