summaryrefslogtreecommitdiff
path: root/lib/libm/src
diff options
context:
space:
mode:
authorchristos <christos@NetBSD.org>2013-11-06 16:49:21 +0000
committerchristos <christos@NetBSD.org>2013-11-06 16:49:21 +0000
commit2b85e8d6b7e585b2effa57ff41a830d803c510cf (patch)
tree082d32395ca421cf1c2440a67af4c77307e059cd /lib/libm/src
parent8c38ca6aed5cd3f7554edbe966130e10c7efc218 (diff)
simplify.
Diffstat (limited to 'lib/libm/src')
-rw-r--r--lib/libm/src/s_exp2.c15
1 files changed, 8 insertions, 7 deletions
diff --git a/lib/libm/src/s_exp2.c b/lib/libm/src/s_exp2.c
index 8d3d2907c7d..05655f8d17d 100644
--- a/lib/libm/src/s_exp2.c
+++ b/lib/libm/src/s_exp2.c
@@ -25,7 +25,7 @@
*/
#include <sys/cdefs.h>
-__RCSID("$NetBSD: s_exp2.c,v 1.2 2010/01/11 23:38:24 christos Exp $");
+__RCSID("$NetBSD: s_exp2.c,v 1.3 2013/11/06 16:49:21 christos Exp $");
#ifdef __FBSDID
__FBSDID("$FreeBSD: src/lib/msun/src/s_exp2.c,v 1.7 2008/02/22 02:27:34 das Exp $");
#endif
@@ -343,9 +343,9 @@ static const double tbl[TBLSIZE * 2] = {
double
exp2(double x)
{
- double r, t, twopk, twopkp1000, z;
+ double r, t, twopk, z;
uint32_t hx, ix, lx, i0;
- int k;
+ int k, big;
/* Filter out exceptional cases. */
GET_HIGH_WORD(hx,x);
@@ -378,19 +378,20 @@ exp2(double x)
/* Compute r = exp2(y) = exp2t[i0] * p(z - eps[i]). */
t = tbl[i0]; /* exp2t[i0] */
z -= tbl[i0 + 1]; /* eps[i0] */
- if (k >= -1021 << 20)
+ big = k >= -1021 << 20;
+ if (big)
INSERT_WORDS(twopk, 0x3ff00000 + k, 0);
else
- INSERT_WORDS(twopkp1000, 0x3ff00000 + k + (1000 << 20), 0);
+ INSERT_WORDS(twopk, 0x3ff00000 + k + (1000 << 20), 0);
r = t + t * z * (P1 + z * (P2 + z * (P3 + z * (P4 + z * P5))));
/* Scale by 2**(k>>20). */
- if(k >= -1021 << 20) {
+ if (big) {
if (k == 1024 << 20)
return (r * 2.0 * 0x1p1023);
return (r * twopk);
} else {
- return (r * twopkp1000 * twom1000);
+ return (r * twopk * twom1000);
}
}