From 43576de04afc6a0896a3ecc094e1581069a0652a Mon Sep 17 00:00:00 2001 From: Paul Zimmermann Date: Fri, 2 Apr 2021 08:21:06 +0200 Subject: Improve the accuracy of tgamma (BZ #26983) With this patch, the maximal known error for tgamma is now reduced to 9 ulps for dbl-64, for all rounding modes. Since exhaustive testing is not possible for dbl-64, it might be that there are still cases with an error larger than 9 ulps, but all known cases are fixed (intensive tests were done to find cases with large errors). Tested on x86_64 and powerpc (and by Adhemerval Zanella on aarch64, arm, s390x, sparc, and i686). Reviewed-by: Adhemerval Zanella --- sysdeps/x86_64/fpu/libm-test-ulps | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) (limited to 'sysdeps/x86_64/fpu/libm-test-ulps') diff --git a/sysdeps/x86_64/fpu/libm-test-ulps b/sysdeps/x86_64/fpu/libm-test-ulps index 0edb95e14c..c621a03268 100644 --- a/sysdeps/x86_64/fpu/libm-test-ulps +++ b/sysdeps/x86_64/fpu/libm-test-ulps @@ -1734,16 +1734,16 @@ float128: 4 ldouble: 5 Function: "tgamma_downward": -double: 8 +double: 9 float: 7 float128: 5 -ldouble: 5 +ldouble: 6 Function: "tgamma_towardzero": double: 9 float: 7 float128: 5 -ldouble: 5 +ldouble: 6 Function: "tgamma_upward": double: 9 -- cgit 1.4.1