about summary refs log tree commit diff
path: root/math/auto-libm-test-in
diff options
context:
space:
mode:
authorPaul Zimmermann <Paul.Zimmermann@inria.fr>2021-04-02 08:21:06 +0200
committerPaul Zimmermann <Paul.Zimmermann@inria.fr>2021-04-07 13:23:39 +0200
commit43576de04afc6a0896a3ecc094e1581069a0652a (patch)
tree42b35efc19ae2a9f22c354176d83173f74818268 /math/auto-libm-test-in
parentd1a3dcabf2f89233a99a4a9be08f9f407da0b6b4 (diff)
downloadglibc-43576de04afc6a0896a3ecc094e1581069a0652a.tar.gz
glibc-43576de04afc6a0896a3ecc094e1581069a0652a.tar.xz
glibc-43576de04afc6a0896a3ecc094e1581069a0652a.zip
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  <adhemerval.zanella@linaro.org>
Diffstat (limited to 'math/auto-libm-test-in')
-rw-r--r--math/auto-libm-test-in6
1 files changed, 5 insertions, 1 deletions
diff --git a/math/auto-libm-test-in b/math/auto-libm-test-in
index 9fbd0c626b..2045a0f8d1 100644
--- a/math/auto-libm-test-in
+++ b/math/auto-libm-test-in
@@ -8251,8 +8251,12 @@ tgamma -0x6.ce9158p+0
 tgamma -0xd.cbf53d0e7d06p+0
 # the next value generates larger error bounds on x86_64 (binary32)
 tgamma -0x3.0aa534p+0
-# the next value generates larger error bounds on x86_64 (binary64)
+# the next values generate large error bounds on x86_64 (binary64)
 tgamma -0x1.62b8c36c7180bp+4
+tgamma -0x1.62c4d519e8677p+3
+tgamma -0x1.c033cc426752fp+2
+tgamma -0x1.62cfd0d34ade2p+3
+tgamma -0x1.8814da6eb7dbp+5
 
 y0 0.125
 y0 0.75