aboutsummaryrefslogtreecommitdiff
path: root/math
diff options
context:
space:
mode:
authorJoseph Myers <joseph@codesourcery.com>2016-06-03 21:30:12 +0000
committerJoseph Myers <joseph@codesourcery.com>2016-06-03 21:30:12 +0000
commit8cbd1453ec1c59e54652edcd88256f4023ff77b9 (patch)
treefd75adff737c0b9172a5e59d2307d36e5190afe8 /math
parentcfac4de69cc70fbd6364e58cf77891a7a30f2e09 (diff)
downloadglibc-8cbd1453ec1c59e54652edcd88256f4023ff77b9.zip
glibc-8cbd1453ec1c59e54652edcd88256f4023ff77b9.tar.gz
glibc-8cbd1453ec1c59e54652edcd88256f4023ff77b9.tar.bz2
Fix x86/x86_64 nextafterl incrementing negative subnormals (bug 20205).
The x86 / x86_64 implementation of nextafterl (also used for nexttowardl) produces incorrect results (NaNs) when negative subnormals, the low 32 bits of whose mantissa are zero, are incremented towards zero. This patch fixes this by disabling the logic to decrement the exponent in that case. Tested for x86_64 and x86. [BZ #20205] * sysdeps/i386/fpu/s_nextafterl.c (__nextafterl): Do not adjust exponent when incrementing negative subnormal with low mantissa word zero. * math/libm-test.inc (nextafter_test_data) [TEST_COND_intel96]: Add another test.
Diffstat (limited to 'math')
-rw-r--r--math/libm-test.inc4
1 files changed, 4 insertions, 0 deletions
diff --git a/math/libm-test.inc b/math/libm-test.inc
index 8c64c8e..aaa0148 100644
--- a/math/libm-test.inc
+++ b/math/libm-test.inc
@@ -9965,6 +9965,10 @@ static const struct test_ff_f_data nextafter_test_data[] =
TEST_ff_f (nextafter, -min_subnorm_value, 0, minus_zero, INEXACT_EXCEPTION|UNDERFLOW_EXCEPTION|ERRNO_ERANGE),
TEST_ff_f (nextafter, -min_subnorm_value, minus_zero, minus_zero, INEXACT_EXCEPTION|UNDERFLOW_EXCEPTION|ERRNO_ERANGE),
+#if TEST_COND_intel96
+ TEST_ff_f (nextafter, -0x0.fffffffep-16382L, 0.0L, -0x0.fffffffdfffffffep-16382L, INEXACT_EXCEPTION|UNDERFLOW_EXCEPTION|ERRNO_ERANGE),
+#endif
+
#if MANT_DIG >= 64
// XXX Enable once gcc is fixed.
//TEST_ff_f (nextafter, 0x0.00000040000000000000p-16385L, -0.1L, 0x0.0000003ffffffff00000p-16385L),