From 526af54142db14c1edcd2d80dc1b56d33ff4e8ce Mon Sep 17 00:00:00 2001 From: Joseph Myers Date: Tue, 19 May 2015 23:05:22 +0000 Subject: Fix i386 atanhl spurious underflows (bug 18049). The i386 implementation of atanhl, for small arguments, does a calculation that involves computing twice the square of the argument, resulting in spurious underflows for some arguments. This patch fixes this by just returning the argument when its exponent is below -32, with underflow being forced as needed for subnormal arguments. Tested for x86 and x86_64. [BZ #18049] * sysdeps/i386/fpu/e_atanhl.S (__ieee754_atanhl): For exponents below -32, return the argument, with underflow if subnormal. * math/auto-libm-test-in: Add more tests of atanh. * math/auto-libm-test-out: Regenerated. --- sysdeps/i386/fpu/e_atanhl.S | 10 ++++++++++ 1 file changed, 10 insertions(+) (limited to 'sysdeps/i386/fpu') diff --git a/sysdeps/i386/fpu/e_atanhl.S b/sysdeps/i386/fpu/e_atanhl.S index 88a970378d..1ad6851058 100644 --- a/sysdeps/i386/fpu/e_atanhl.S +++ b/sysdeps/i386/fpu/e_atanhl.S @@ -56,6 +56,16 @@ ENTRY(__ieee754_atanhl) andl $0x7fff, %eax cmpl $0x7fff, %eax je 5f + cmpl $0x3fdf, %eax + jge 7f + // Exponent below -32; return x, with underflow if subnormal. + fldt 4(%esp) + cmpl $0, %eax + jne 8f + fld %st(0) + fmul %st(0) + fstp %st(0) +8: ret 7: #ifdef PIC -- cgit 1.4.1