Fix log1p missing underflows (bug 16339).
authorJoseph Myers <joseph@codesourcery.com>
Thu, 14 May 2015 23:38:07 +0000 (23:38 +0000)
committerJoseph Myers <joseph@codesourcery.com>
Thu, 14 May 2015 23:38:07 +0000 (23:38 +0000)
commit0b7a5f920163d03806d7c5d9d1c83b16942c9496
treea108347db6f5f094bf917214f562b9a3a7ae8b04
parent95b07fbcc702a83b421cb5d0488ba2817d3123c1
Fix log1p missing underflows (bug 16339).

Similar to various other bugs in this area, some log1p implementations
do not raise the underflow exception for subnormal arguments, when the
result is tiny and inexact.  This patch forces the exception in a
similar way to previous fixes.  (The ldbl-128ibm implementation
doesn't currently need any change as it already generates this
exception, albeit through code that would generate spurious exceptions
in other cases; special code for this issue will only be needed there
when fixing the spurious exceptions.)

Tested for x86_64, x86, powerpc and mips64.

[BZ #16339]
* sysdeps/i386/fpu/s_log1p.S (dbl_min): New object.
(__log1p): Force underflow exception for results with small
absolute value.
* sysdeps/i386/fpu/s_log1pf.S (flt_min): New object.
(__log1pf): Force underflow exception for results with small
absolute value.
* sysdeps/ieee754/dbl-64/s_log1p.c: Include <float.h>.
(__log1p): Force underflow exception for results with small
absolute value.
* sysdeps/ieee754/flt-32/s_log1pf.c: Include <float.h>.
(__log1pf): Force underflow exception for results with small
absolute value.
* sysdeps/ieee754/ldbl-128/s_log1pl.c: Include <float.h>.
(__log1pl): Force underflow exception for results with small
absolute value.
* math/auto-libm-test-in: Do not allow missing underflow
exceptions from log1p.
* math/auto-libm-test-out: Regenerated.
ChangeLog
NEWS
math/auto-libm-test-in
math/auto-libm-test-out
sysdeps/i386/fpu/s_log1p.S
sysdeps/i386/fpu/s_log1pf.S
sysdeps/ieee754/dbl-64/s_log1p.c
sysdeps/ieee754/flt-32/s_log1pf.c
sysdeps/ieee754/ldbl-128/s_log1pl.c