aboutsummaryrefslogtreecommitdiff
path: root/sysdeps/i386/fpu/s_frexpl.S
diff options
context:
space:
mode:
authorJoseph Myers <joseph@codesourcery.com>2016-06-13 17:27:19 +0000
committerJoseph Myers <joseph@codesourcery.com>2016-06-13 17:27:19 +0000
commit88283451b26de29d08ba33148cf2f8cb2680f6f0 (patch)
tree334945ab68b411d9bac8696fbcaba40a6598ed2b /sysdeps/i386/fpu/s_frexpl.S
parentb39b6e0c902f0bf5f99d593b8f5953943777c610 (diff)
downloadglibc-88283451b26de29d08ba33148cf2f8cb2680f6f0.zip
glibc-88283451b26de29d08ba33148cf2f8cb2680f6f0.tar.gz
glibc-88283451b26de29d08ba33148cf2f8cb2680f6f0.tar.bz2
Fix frexp (NaN) (bug 20250).
Various implementations of frexp functions return sNaN for sNaN input. This patch fixes them to add such arguments to themselves so that qNaN is returned. Tested for x86_64, x86, mips64 and powerpc. [BZ #20250] * sysdeps/i386/fpu/s_frexpl.S (__frexpl): Add non-finite input to itself. * sysdeps/ieee754/dbl-64/s_frexp.c (__frexp): Add non-finite or zero input to itself. * sysdeps/ieee754/dbl-64/wordsize-64/s_frexp.c (__frexp): Likewise. * sysdeps/ieee754/flt-32/s_frexpf.c (__frexpf): Likewise. * sysdeps/ieee754/ldbl-128/s_frexpl.c (__frexpl): Likewise. * sysdeps/ieee754/ldbl-128ibm/s_frexpl.c (__frexpl): Likewise. * sysdeps/ieee754/ldbl-96/s_frexpl.c (__frexpl): Likewise. * math/libm-test.inc (frexp_test_data): Add sNaN tests.
Diffstat (limited to 'sysdeps/i386/fpu/s_frexpl.S')
-rw-r--r--sysdeps/i386/fpu/s_frexpl.S9
1 files changed, 8 insertions, 1 deletions
diff --git a/sysdeps/i386/fpu/s_frexpl.S b/sysdeps/i386/fpu/s_frexpl.S
index f687ab3..b33a8a8 100644
--- a/sysdeps/i386/fpu/s_frexpl.S
+++ b/sysdeps/i386/fpu/s_frexpl.S
@@ -50,7 +50,7 @@ ENTRY (__frexpl)
jz 1f
xorl %ecx, %ecx
cmpl $0x7fff, %eax
- je 1f
+ je 3f
cmpl $0, %eax
jne 2f
@@ -81,5 +81,12 @@ ENTRY (__frexpl)
movl %ecx, (%eax)
ret
+
+ /* Infinity or NaN; ensure signaling NaNs are quieted. */
+3: movl EXPP(%esp), %eax
+ fldt VAL0(%esp)
+ fadd %st
+ movl %ecx, (%eax)
+ ret
END (__frexpl)
weak_alias (__frexpl, frexpl)