diff options
| author | christos <christos@NetBSD.org> | 2013-11-06 16:49:21 +0000 |
|---|---|---|
| committer | christos <christos@NetBSD.org> | 2013-11-06 16:49:21 +0000 |
| commit | 2b85e8d6b7e585b2effa57ff41a830d803c510cf (patch) | |
| tree | 082d32395ca421cf1c2440a67af4c77307e059cd /lib/libm/src | |
| parent | 8c38ca6aed5cd3f7554edbe966130e10c7efc218 (diff) | |
simplify.
Diffstat (limited to 'lib/libm/src')
| -rw-r--r-- | lib/libm/src/s_exp2.c | 15 |
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); } } |
