Fix i386 exp missing underflows (bug 18961).
authorJoseph Myers <joseph@codesourcery.com>
Mon, 14 Sep 2015 22:40:05 +0000 (22:40 +0000)
committerJoseph Myers <joseph@codesourcery.com>
Mon, 14 Sep 2015 22:40:05 +0000 (22:40 +0000)
commitde5e81691c80468ff97c5260b4020aeaecfe418d
tree210c246a20f943215e12e05ba41f105f3c628596
parent903af5af9a0404c1dfeb7b4db6d5a23a1e45593f
Fix i386 exp missing underflows (bug 18961).

On i386, the double version of exp can miss underflow exceptions if
the result is in the subnormal range for double but the last 11 bits
of the 64-bit extended-precision mantissa happen to be zero.  This
patch forces the exception in a similar way to previous fixes.

As with the exp2 fixes, the expf changes may in fact not be needed to
ensure underflow exceptions, but are included for consistency and to
fix the exp part of bug 18875 by ensuring that excess range and
precision is removed from underflowing return values.

Tested for x86_64 and x86.

[BZ #18875]
[BZ #18961]
* sysdeps/i386/fpu/e_exp.S (dbl_min): New object.
(MO): New macro.
(__ieee754_exp): For small results, force underflow exception and
remove excess range and precision from return value.
(__exp_finite): Likewise.
* sysdeps/i386/fpu/e_expf.S (flt_min): New object.
(MO): New macro.
(__ieee754_expf): For small results, force underflow exception and
remove excess range and precision from return value.
(__expf_finite): Likewise.
* math/auto-libm-test-in: Add more tests of exp.
* math/auto-libm-test-out: Regenerated.
ChangeLog
NEWS
math/auto-libm-test-in
math/auto-libm-test-out
sysdeps/i386/fpu/e_exp.S
sysdeps/i386/fpu/e_expf.S