From 8c6c9236361fbc077769673c259828216403bc33 Mon Sep 17 00:00:00 2001 From: Joseph Myers Date: Wed, 30 Sep 2015 21:44:42 +0000 Subject: Fix i386 acosh (-qNaN) spurious "invalid" exception. The i386 versions of acoshf and acosh raise a spurious "invalid" exception for an argument that is a quiet NaN with the sign bit set. The integer arithmetic to detect arguments < 1 also detects -NaN, and then the computation 0 / 0 in that case raises the exception. This patch fixes this by using (x - x) / (x - x) as the computation in that case instead, which will always raise the exception for non-NaN arguments reaching that code, but not for quiet NaN arguments. Tested for x86_64 and x86. [BZ #19032] * sysdeps/i386/fpu/e_acosh.S (__ieee754_acosh): For arguments < 1, compute result as (x - x) / (x - x) not as 0 / 0. * sysdeps/i386/fpu/e_acoshf.S (__ieee754_acoshf): Likewise. * math/libm-test.inc (acosh_test_data): Add another test of acosh. --- ChangeLog | 6 ++++++ 1 file changed, 6 insertions(+) (limited to 'ChangeLog') diff --git a/ChangeLog b/ChangeLog index ea038f764f..958b4a3a9d 100644 --- a/ChangeLog +++ b/ChangeLog @@ -1,5 +1,11 @@ 2015-09-30 Joseph Myers + [BZ #19032] + * sysdeps/i386/fpu/e_acosh.S (__ieee754_acosh): For arguments < 1, + compute result as (x - x) / (x - x) not as 0 / 0. + * sysdeps/i386/fpu/e_acoshf.S (__ieee754_acoshf): Likewise. + * math/libm-test.inc (acosh_test_data): Add another test of acosh. + * math/auto-libm-test-in: Add more tests of acos, acosh, asin, atan, atan2, atanh, cbrt, cos, cosh, erf, erfc, exp, exp10, exp2 and expm1. -- cgit 1.4.1