diff options
| author | dim <dim@FreeBSD.org> | 2015-08-09 10:00:13 +0000 |
|---|---|---|
| committer | dim <dim@FreeBSD.org> | 2015-08-09 10:00:13 +0000 |
| commit | bcea9d50b15e4f0027a5dd526e0e2a612238471e (patch) | |
| tree | 04301b54fbc56c06e09f6cae83bac5a34da59033 /lib/msun | |
| parent | 30a2b84ce8d55ac387fdd09308fec2ac49bbd03c (diff) | |
| download | FreeBSD-src-bcea9d50b15e4f0027a5dd526e0e2a612238471e.zip FreeBSD-src-bcea9d50b15e4f0027a5dd526e0e2a612238471e.tar.gz | |
In libm's exp2(3), avoid left-shifting a negative integer, which is
undefined. Replace it with the intended value, in a defined way.
Reviewed by: bde
MFC after: 3 days
Diffstat (limited to 'lib/msun')
| -rw-r--r-- | lib/msun/src/s_exp2.c | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/lib/msun/src/s_exp2.c b/lib/msun/src/s_exp2.c index fde11c2..5cf2e15 100644 --- a/lib/msun/src/s_exp2.c +++ b/lib/msun/src/s_exp2.c @@ -376,14 +376,14 @@ 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) + 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 >= -(1021 << 20)) { if (k == 1024 << 20) return (r * 2.0 * 0x1p1023); return (r * twopk); |
