diff options
Diffstat (limited to 'lib/msun/src/s_exp2.c')
| -rw-r--r-- | lib/msun/src/s_exp2.c | 18 |
1 files changed, 9 insertions, 9 deletions
diff --git a/lib/msun/src/s_exp2.c b/lib/msun/src/s_exp2.c index 56554f1..63b8997 100644 --- a/lib/msun/src/s_exp2.c +++ b/lib/msun/src/s_exp2.c @@ -340,7 +340,7 @@ static const double tbl[TBLSIZE * 2] = { double exp2(double x) { - double r, t, z; + double r, t, twopk, twopkp1000, z; uint32_t hx, hr, ix, lx, i0; int k; @@ -375,19 +375,19 @@ 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) + INSERT_WORDS(twopk, 0x3ff00000 + k, 0); + else + INSERT_WORDS(twopkp1000, 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 (k != 0) { - GET_HIGH_WORD(hr, r); - SET_HIGH_WORD(r, hr + k); - } - return (r); + if (k == 1024 << 20) + return (r * 2.0 * 0x1p1023); + return (r * twopk); } else { - GET_HIGH_WORD(hr, r); - SET_HIGH_WORD(r, hr + (k + (1000 << 20))); - return (r * twom1000); + return (r * twopkp1000 * twom1000); } } |
