return one + (t*t)/(w+w);
}
- /* |x| in [0.5*ln2,22], return (exp(|x|)+1/exp(|x|)/2; */
+ /* |x| in [0.5*ln2,22], return (exp(|x|)+1/exp(|x|))/2; */
if (ix < 0x40360000) {
t = exp(fabs(x));
return half*t + half/t;
return half*exp(fabs(x));
/* |x| in [log(maxdouble), overflowthresold] */
- if (ix<=0x408633CE)
+ if (ix <= 0x408633CE)
return __ldexp_exp(fabs(x), -1);
/* |x| > overflowthresold, cosh(x) overflow */