summaryrefslogtreecommitdiff
path: root/libm/e_exp.c
diff options
context:
space:
mode:
Diffstat (limited to 'libm/e_exp.c')
-rw-r--r--libm/e_exp.c25
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)