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