diff options
Diffstat (limited to 'libm/e_exp.c')
-rw-r--r-- | libm/e_exp.c | 25 |
1 files changed, 25 insertions, 0 deletions
diff --git a/libm/e_exp.c b/libm/e_exp.c index ebf2e5af0..abfb17675 100644 --- a/libm/e_exp.c +++ b/libm/e_exp.c @@ -155,3 +155,28 @@ double attribute_hidden __ieee754_exp(double x) /* default IEEE double exp */ return y*twom1000; } } + +/* + * wrapper exp(x) + */ +#ifndef _IEEE_LIBM +double exp(double x) +{ + static const double o_threshold = 7.09782712893383973096e+02; /* 0x40862E42, 0xFEFA39EF */ + static const double u_threshold = -7.45133219101941108420e+02; /* 0xc0874910, 0xD52D3051 */ + + double z = __ieee754_exp(x); + if (_LIB_VERSION == _IEEE_) + return z; + if (isfinite(x)) { + if (x > o_threshold) + return __kernel_standard(x, x, 6); /* exp overflow */ + if (x < u_threshold) + return __kernel_standard(x, x, 7); /* exp underflow */ + } + return z; +} +#else +strong_alias(__ieee754_exp, exp) +#endif +libm_hidden_def(exp) |