From 3d8b06bc6108aa20a6591f52e7d9cbbd63c20e0a Mon Sep 17 00:00:00 2001 From: Joseph Myers Date: Mon, 6 Jun 2016 22:21:11 +0000 Subject: Fix dbl-64 asin (sNaN) (bug 20213). The dbl-64 version of asin returns sNaN for sNaN arguments. This patch fixes it to add NaN arguments to themselves so that qNaN is returned in this case. Tested for x86_64 and x86. [BZ #20213] * sysdeps/ieee754/dbl-64/e_asin.c (__ieee754_asin): Add NaN argument to itself. * math/libm-test.inc (asin_test_data): Add sNaN tests. --- math/libm-test.inc | 2 ++ 1 file changed, 2 insertions(+) (limited to 'math/libm-test.inc') diff --git a/math/libm-test.inc b/math/libm-test.inc index 6981720..f9b2c8a 100644 --- a/math/libm-test.inc +++ b/math/libm-test.inc @@ -1853,6 +1853,8 @@ static const struct test_f_f_data asin_test_data[] = TEST_f_f (asin, minus_infty, qnan_value, INVALID_EXCEPTION|ERRNO_EDOM), TEST_f_f (asin, qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED), TEST_f_f (asin, -qnan_value, qnan_value, NO_INEXACT_EXCEPTION|ERRNO_UNCHANGED), + TEST_f_f (asin, snan_value, qnan_value, INVALID_EXCEPTION), + TEST_f_f (asin, -snan_value, qnan_value, INVALID_EXCEPTION), /* asin x == qNaN plus invalid exception for |x| > 1. */ TEST_f_f (asin, 1.125L, qnan_value, INVALID_EXCEPTION|ERRNO_EDOM), -- cgit v1.1