From e02920bc029019443326eecaa7b267b78ff2892e Mon Sep 17 00:00:00 2001 From: Joseph Myers Date: Mon, 29 Jun 2015 23:29:35 +0000 Subject: Improve tgamma accuracy (bug 18613). In non-default rounding modes, tgamma can be slightly less accurate than permitted by glibc's accuracy goals. Part of the problem is error accumulation, addressed in this patch by setting round-to-nearest for internal computations. However, there was also a bug in the code dealing with computing pow (x + n, x + n) where x + n is not exactly representable, providing another source of error even in round-to-nearest mode; it was necessary to address both bugs to get errors for all testcases within glibc's accuracy goals. Given this second fix, accuracy in round-to-nearest mode is also improved (hence regeneration of ulps for tgamma should be from scratch - truncate libm-test-ulps or at least remove existing tgamma entries - so that the expected ulps can be reduced). Some additional complications also arose. Certain tgamma tests should strictly, according to IEEE semantics, overflow or not depending on the rounding mode; this is beyond the scope of glibc's accuracy goals for any function without exactly-determined results, but gen-auto-libm-tests doesn't handle being lax there as it does for underflow. (libm-test.inc also doesn't handle being lax about whether the result in cases very close to the overflow threshold is infinity or a finite value close to overflow, but that doesn't cause problems in this case though I've seen it cause problems with random test generation for some functions.) Thus, spurious-overflow markings, with a comment, are added to auto-libm-test-in (no bug in Bugzilla because the issue is with the testsuite, not a user-visible bug in glibc). And on x86, after the patch I saw ERANGE issues as previously reported by Carlos (see my commentary in ), which needed addressing by ensuring excess range and precision were eliminated at various points if FLT_EVAL_METHOD != 0. I also noticed and fixed a cosmetic issue where 1.0f was used in long double functions and should have been 1.0L. This completes the move of all functions to testing in all rounding modes with ALL_RM_TEST, so gen-libm-have-vector-test.sh is updated to remove the workaround for some functions not using ALL_RM_TEST. Tested for x86_64, x86, mips64 and powerpc. [BZ #18613] * sysdeps/ieee754/dbl-64/e_gamma_r.c (gamma_positive): Take log of X_ADJ not X when adjusting exponent. (__ieee754_gamma_r): Do intermediate computations in round-to-nearest then adjust overflowing and underflowing results as needed. * sysdeps/ieee754/flt-32/e_gammaf_r.c (gammaf_positive): Take log of X_ADJ not X when adjusting exponent. (__ieee754_gammaf_r): Do intermediate computations in round-to-nearest then adjust overflowing and underflowing results as needed. * sysdeps/ieee754/ldbl-128/e_gammal_r.c (gammal_positive): Take log of X_ADJ not X when adjusting exponent. (__ieee754_gammal_r): Do intermediate computations in round-to-nearest then adjust overflowing and underflowing results as needed. Use 1.0L not 1.0f as numerator of division. * sysdeps/ieee754/ldbl-128ibm/e_gammal_r.c (gammal_positive): Take log of X_ADJ not X when adjusting exponent. (__ieee754_gammal_r): Do intermediate computations in round-to-nearest then adjust overflowing and underflowing results as needed. Use 1.0L not 1.0f as numerator of division. * sysdeps/ieee754/ldbl-96/e_gammal_r.c (gammal_positive): Take log of X_ADJ not X when adjusting exponent. (__ieee754_gammal_r): Do intermediate computations in round-to-nearest then adjust overflowing and underflowing results as needed. Use 1.0L not 1.0f as numerator of division. * math/libm-test.inc (tgamma_test_data): Remove one test. Moved to auto-libm-test-in. (tgamma_test): Use ALL_RM_TEST. * math/auto-libm-test-in: Add one test of tgamma. Mark some other tests of tgamma with spurious-overflow. * math/auto-libm-test-out: Regenerated. * math/gen-libm-have-vector-test.sh: Do not check for START. * sysdeps/i386/fpu/libm-test-ulps: Update. * sysdeps/x86_64/fpu/libm-test-ulps: Likewise. --- math/auto-libm-test-in | 10 +- math/auto-libm-test-out | 263 ++++++++++++++++++++++++-------------- math/gen-libm-have-vector-test.sh | 8 -- math/libm-test.inc | 5 +- 4 files changed, 174 insertions(+), 112 deletions(-) (limited to 'math') diff --git a/math/auto-libm-test-in b/math/auto-libm-test-in index f4313a2b99..f783c5b2e8 100644 --- a/math/auto-libm-test-in +++ b/math/auto-libm-test-in @@ -2682,19 +2682,22 @@ tgamma 0x1p-113 tgamma -0x1p-113 tgamma 0x1p-127 tgamma -0x1p-127 -tgamma 0x1p-128 +# IEEE semantics mean overflow very close to the threshold depends on +# the rounding mode; gen-auto-libm-tests does not reflect that glibc +# does not try to achieve this. +tgamma 0x1p-128 spurious-overflow:flt-32 tgamma -0x1p-128 tgamma 0x1p-149 tgamma -0x1p-149 tgamma 0x1p-1023 tgamma -0x1p-1023 -tgamma 0x1p-1024 +tgamma 0x1p-1024 spurious-overflow:dbl-64 spurious-overflow:ldbl-128ibm tgamma -0x1p-1024 tgamma 0x1p-1074 tgamma -0x1p-1074 tgamma 0x1p-16383 tgamma -0x1p-16383 -tgamma 0x1p-16384 +tgamma 0x1p-16384 spurious-overflow:ldbl-96-intel spurious-overflow:ldbl-96-m68k spurious-overflow:ldbl-128 tgamma -0x1p-16384 tgamma 0x1p-16445 tgamma -0x1p-16445 @@ -3075,6 +3078,7 @@ tgamma 0x6.db8c603359a971081bc4a2e9dfdp+8 tgamma 0x6.db8c603359a971081bc4a2e9dfd4p+8 tgamma 1e3 tgamma -100000.5 +tgamma max tgamma -0x3.06644cp+0 tgamma -0x6.fe4636e0c5064p+0 diff --git a/math/auto-libm-test-out b/math/auto-libm-test-out index b71429b4ac..14fbd03c94 100644 --- a/math/auto-libm-test-out +++ b/math/auto-libm-test-out @@ -197134,31 +197134,31 @@ tgamma -0x1p-127 = tgamma tonearest ldbl-128ibm -0x2p-128L : -0x8p+124L : inexact-ok = tgamma towardzero ldbl-128ibm -0x2p-128L : -0x8p+124L : inexact-ok = tgamma upward ldbl-128ibm -0x2p-128L : -0x8p+124L : inexact-ok -tgamma 0x1p-128 -= tgamma downward flt-32 0x1p-128f : 0xf.fffffp+124f : inexact-ok +tgamma 0x1p-128 spurious-overflow:flt-32 += tgamma downward flt-32 0x1p-128f : 0xf.fffffp+124f : inexact-ok overflow-ok:flt-32 = tgamma tonearest flt-32 0x1p-128f : plus_infty : inexact-ok overflow errno-erange -= tgamma towardzero flt-32 0x1p-128f : 0xf.fffffp+124f : inexact-ok += tgamma towardzero flt-32 0x1p-128f : 0xf.fffffp+124f : inexact-ok overflow-ok:flt-32 = tgamma upward flt-32 0x1p-128f : plus_infty : inexact-ok overflow errno-erange -= tgamma downward dbl-64 0x1p-128 : 0xf.ffffffffffff8p+124 : inexact-ok -= tgamma tonearest dbl-64 0x1p-128 : 0x1p+128 : inexact-ok -= tgamma towardzero dbl-64 0x1p-128 : 0xf.ffffffffffff8p+124 : inexact-ok -= tgamma upward dbl-64 0x1p-128 : 0x1p+128 : inexact-ok -= tgamma downward ldbl-96-intel 0x1p-128L : 0xf.fffffffffffffffp+124L : inexact-ok -= tgamma tonearest ldbl-96-intel 0x1p-128L : 0x1p+128L : inexact-ok -= tgamma towardzero ldbl-96-intel 0x1p-128L : 0xf.fffffffffffffffp+124L : inexact-ok -= tgamma upward ldbl-96-intel 0x1p-128L : 0x1p+128L : inexact-ok -= tgamma downward ldbl-96-m68k 0x1p-128L : 0xf.fffffffffffffffp+124L : inexact-ok -= tgamma tonearest ldbl-96-m68k 0x1p-128L : 0x1p+128L : inexact-ok -= tgamma towardzero ldbl-96-m68k 0x1p-128L : 0xf.fffffffffffffffp+124L : inexact-ok -= tgamma upward ldbl-96-m68k 0x1p-128L : 0x1p+128L : inexact-ok -= tgamma downward ldbl-128 0x1p-128L : 0xf.fffffffffffffffffffffffffff8p+124L : inexact-ok -= tgamma tonearest ldbl-128 0x1p-128L : 0x1p+128L : inexact-ok -= tgamma towardzero ldbl-128 0x1p-128L : 0xf.fffffffffffffffffffffffffff8p+124L : inexact-ok -= tgamma upward ldbl-128 0x1p-128L : 0x1p+128L : inexact-ok -= tgamma downward ldbl-128ibm 0x1p-128L : 0xf.fffffffffffffffffffffffffcp+124L : inexact-ok -= tgamma tonearest ldbl-128ibm 0x1p-128L : 0x1p+128L : inexact-ok -= tgamma towardzero ldbl-128ibm 0x1p-128L : 0xf.fffffffffffffffffffffffffcp+124L : inexact-ok -= tgamma upward ldbl-128ibm 0x1p-128L : 0x1p+128L : inexact-ok += tgamma downward dbl-64 0x1p-128 : 0xf.ffffffffffff8p+124 : inexact-ok overflow-ok:flt-32 += tgamma tonearest dbl-64 0x1p-128 : 0x1p+128 : inexact-ok overflow-ok:flt-32 += tgamma towardzero dbl-64 0x1p-128 : 0xf.ffffffffffff8p+124 : inexact-ok overflow-ok:flt-32 += tgamma upward dbl-64 0x1p-128 : 0x1p+128 : inexact-ok overflow-ok:flt-32 += tgamma downward ldbl-96-intel 0x1p-128L : 0xf.fffffffffffffffp+124L : inexact-ok overflow-ok:flt-32 += tgamma tonearest ldbl-96-intel 0x1p-128L : 0x1p+128L : inexact-ok overflow-ok:flt-32 += tgamma towardzero ldbl-96-intel 0x1p-128L : 0xf.fffffffffffffffp+124L : inexact-ok overflow-ok:flt-32 += tgamma upward ldbl-96-intel 0x1p-128L : 0x1p+128L : inexact-ok overflow-ok:flt-32 += tgamma downward ldbl-96-m68k 0x1p-128L : 0xf.fffffffffffffffp+124L : inexact-ok overflow-ok:flt-32 += tgamma tonearest ldbl-96-m68k 0x1p-128L : 0x1p+128L : inexact-ok overflow-ok:flt-32 += tgamma towardzero ldbl-96-m68k 0x1p-128L : 0xf.fffffffffffffffp+124L : inexact-ok overflow-ok:flt-32 += tgamma upward ldbl-96-m68k 0x1p-128L : 0x1p+128L : inexact-ok overflow-ok:flt-32 += tgamma downward ldbl-128 0x1p-128L : 0xf.fffffffffffffffffffffffffff8p+124L : inexact-ok overflow-ok:flt-32 += tgamma tonearest ldbl-128 0x1p-128L : 0x1p+128L : inexact-ok overflow-ok:flt-32 += tgamma towardzero ldbl-128 0x1p-128L : 0xf.fffffffffffffffffffffffffff8p+124L : inexact-ok overflow-ok:flt-32 += tgamma upward ldbl-128 0x1p-128L : 0x1p+128L : inexact-ok overflow-ok:flt-32 += tgamma downward ldbl-128ibm 0x1p-128L : 0xf.fffffffffffffffffffffffffcp+124L : inexact-ok overflow-ok:flt-32 += tgamma tonearest ldbl-128ibm 0x1p-128L : 0x1p+128L : inexact-ok overflow-ok:flt-32 += tgamma towardzero ldbl-128ibm 0x1p-128L : 0xf.fffffffffffffffffffffffffcp+124L : inexact-ok overflow-ok:flt-32 += tgamma upward ldbl-128ibm 0x1p-128L : 0x1p+128L : inexact-ok overflow-ok:flt-32 tgamma -0x1p-128 = tgamma downward flt-32 -0x1p-128f : minus_infty : inexact-ok overflow errno-erange = tgamma tonearest flt-32 -0x1p-128f : minus_infty : inexact-ok overflow errno-erange @@ -197324,47 +197324,47 @@ tgamma -0x1p-1023 = tgamma tonearest ldbl-128ibm -0x2p-1024L : -0x8p+1020L : inexact-ok = tgamma towardzero ldbl-128ibm -0x2p-1024L : -0x8p+1020L : inexact-ok = tgamma upward ldbl-128ibm -0x2p-1024L : -0x8p+1020L : inexact-ok -tgamma 0x1p-1024 +tgamma 0x1p-1024 spurious-overflow:dbl-64 spurious-overflow:ldbl-128ibm = tgamma downward flt-32 0x8p-152f : 0xf.fffffp+124f : inexact-ok overflow errno-erange-ok = tgamma tonearest flt-32 0x8p-152f : plus_infty : inexact-ok overflow errno-erange = tgamma towardzero flt-32 0x8p-152f : 0xf.fffffp+124f : inexact-ok overflow errno-erange-ok = tgamma upward flt-32 0x8p-152f : plus_infty : inexact-ok overflow errno-erange -= tgamma downward dbl-64 0x8p-152 : 0x1.fffffffffffffp+148 : inexact-ok -= tgamma tonearest dbl-64 0x8p-152 : 0x2p+148 : inexact-ok -= tgamma towardzero dbl-64 0x8p-152 : 0x1.fffffffffffffp+148 : inexact-ok -= tgamma upward dbl-64 0x8p-152 : 0x2p+148 : inexact-ok -= tgamma downward ldbl-96-intel 0x8p-152L : 0x1.fffffffffffffffep+148L : inexact-ok -= tgamma tonearest ldbl-96-intel 0x8p-152L : 0x2p+148L : inexact-ok -= tgamma towardzero ldbl-96-intel 0x8p-152L : 0x1.fffffffffffffffep+148L : inexact-ok -= tgamma upward ldbl-96-intel 0x8p-152L : 0x2p+148L : inexact-ok -= tgamma downward ldbl-96-m68k 0x8p-152L : 0x1.fffffffffffffffep+148L : inexact-ok -= tgamma tonearest ldbl-96-m68k 0x8p-152L : 0x2p+148L : inexact-ok -= tgamma towardzero ldbl-96-m68k 0x8p-152L : 0x1.fffffffffffffffep+148L : inexact-ok -= tgamma upward ldbl-96-m68k 0x8p-152L : 0x2p+148L : inexact-ok -= tgamma downward ldbl-128 0x8p-152L : 0x1.ffffffffffffffffffffffffffffp+148L : inexact-ok -= tgamma tonearest ldbl-128 0x8p-152L : 0x2p+148L : inexact-ok -= tgamma towardzero ldbl-128 0x8p-152L : 0x1.ffffffffffffffffffffffffffffp+148L : inexact-ok -= tgamma upward ldbl-128 0x8p-152L : 0x2p+148L : inexact-ok -= tgamma downward ldbl-128ibm 0x8p-152L : 0x1.ffffffffffffffffffffffffff8p+148L : inexact-ok -= tgamma tonearest ldbl-128ibm 0x8p-152L : 0x2p+148L : inexact-ok -= tgamma towardzero ldbl-128ibm 0x8p-152L : 0x1.ffffffffffffffffffffffffff8p+148L : inexact-ok -= tgamma upward ldbl-128ibm 0x8p-152L : 0x2p+148L : inexact-ok -= tgamma downward dbl-64 0x1p-1024 : 0xf.ffffffffffff8p+1020 : inexact-ok += tgamma downward dbl-64 0x8p-152 : 0x1.fffffffffffffp+148 : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma tonearest dbl-64 0x8p-152 : 0x2p+148 : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma towardzero dbl-64 0x8p-152 : 0x1.fffffffffffffp+148 : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma upward dbl-64 0x8p-152 : 0x2p+148 : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma downward ldbl-96-intel 0x8p-152L : 0x1.fffffffffffffffep+148L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma tonearest ldbl-96-intel 0x8p-152L : 0x2p+148L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma towardzero ldbl-96-intel 0x8p-152L : 0x1.fffffffffffffffep+148L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma upward ldbl-96-intel 0x8p-152L : 0x2p+148L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma downward ldbl-96-m68k 0x8p-152L : 0x1.fffffffffffffffep+148L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma tonearest ldbl-96-m68k 0x8p-152L : 0x2p+148L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma towardzero ldbl-96-m68k 0x8p-152L : 0x1.fffffffffffffffep+148L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma upward ldbl-96-m68k 0x8p-152L : 0x2p+148L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma downward ldbl-128 0x8p-152L : 0x1.ffffffffffffffffffffffffffffp+148L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma tonearest ldbl-128 0x8p-152L : 0x2p+148L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma towardzero ldbl-128 0x8p-152L : 0x1.ffffffffffffffffffffffffffffp+148L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma upward ldbl-128 0x8p-152L : 0x2p+148L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma downward ldbl-128ibm 0x8p-152L : 0x1.ffffffffffffffffffffffffff8p+148L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma tonearest ldbl-128ibm 0x8p-152L : 0x2p+148L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma towardzero ldbl-128ibm 0x8p-152L : 0x1.ffffffffffffffffffffffffff8p+148L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma upward ldbl-128ibm 0x8p-152L : 0x2p+148L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma downward dbl-64 0x1p-1024 : 0xf.ffffffffffff8p+1020 : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm = tgamma tonearest dbl-64 0x1p-1024 : plus_infty : inexact-ok overflow errno-erange -= tgamma towardzero dbl-64 0x1p-1024 : 0xf.ffffffffffff8p+1020 : inexact-ok += tgamma towardzero dbl-64 0x1p-1024 : 0xf.ffffffffffff8p+1020 : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm = tgamma upward dbl-64 0x1p-1024 : plus_infty : inexact-ok overflow errno-erange -= tgamma downward ldbl-96-intel 0x1p-1024L : 0xf.fffffffffffffffp+1020L : inexact-ok -= tgamma tonearest ldbl-96-intel 0x1p-1024L : 0x1p+1024L : inexact-ok -= tgamma towardzero ldbl-96-intel 0x1p-1024L : 0xf.fffffffffffffffp+1020L : inexact-ok -= tgamma upward ldbl-96-intel 0x1p-1024L : 0x1p+1024L : inexact-ok -= tgamma downward ldbl-96-m68k 0x1p-1024L : 0xf.fffffffffffffffp+1020L : inexact-ok -= tgamma tonearest ldbl-96-m68k 0x1p-1024L : 0x1p+1024L : inexact-ok -= tgamma towardzero ldbl-96-m68k 0x1p-1024L : 0xf.fffffffffffffffp+1020L : inexact-ok -= tgamma upward ldbl-96-m68k 0x1p-1024L : 0x1p+1024L : inexact-ok -= tgamma downward ldbl-128 0x1p-1024L : 0xf.fffffffffffffffffffffffffff8p+1020L : inexact-ok -= tgamma tonearest ldbl-128 0x1p-1024L : 0x1p+1024L : inexact-ok -= tgamma towardzero ldbl-128 0x1p-1024L : 0xf.fffffffffffffffffffffffffff8p+1020L : inexact-ok -= tgamma upward ldbl-128 0x1p-1024L : 0x1p+1024L : inexact-ok += tgamma downward ldbl-96-intel 0x1p-1024L : 0xf.fffffffffffffffp+1020L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma tonearest ldbl-96-intel 0x1p-1024L : 0x1p+1024L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma towardzero ldbl-96-intel 0x1p-1024L : 0xf.fffffffffffffffp+1020L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma upward ldbl-96-intel 0x1p-1024L : 0x1p+1024L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma downward ldbl-96-m68k 0x1p-1024L : 0xf.fffffffffffffffp+1020L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma tonearest ldbl-96-m68k 0x1p-1024L : 0x1p+1024L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma towardzero ldbl-96-m68k 0x1p-1024L : 0xf.fffffffffffffffp+1020L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma upward ldbl-96-m68k 0x1p-1024L : 0x1p+1024L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma downward ldbl-128 0x1p-1024L : 0xf.fffffffffffffffffffffffffff8p+1020L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma tonearest ldbl-128 0x1p-1024L : 0x1p+1024L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma towardzero ldbl-128 0x1p-1024L : 0xf.fffffffffffffffffffffffffff8p+1020L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm += tgamma upward ldbl-128 0x1p-1024L : 0x1p+1024L : inexact-ok overflow-ok:dbl-64 overflow-ok:ldbl-128ibm = tgamma downward ldbl-128ibm 0x1p-1024L : 0xf.ffffffffffffbffffffffffffcp+1020L : inexact-ok overflow errno-erange-ok = tgamma tonearest ldbl-128ibm 0x1p-1024L : plus_infty : inexact-ok overflow errno-erange = tgamma towardzero ldbl-128ibm 0x1p-1024L : 0xf.ffffffffffffbffffffffffffcp+1020L : inexact-ok overflow errno-erange-ok @@ -197618,62 +197618,62 @@ tgamma -0x1p-16383 = tgamma tonearest ldbl-128 -0x2p-16384L : -0x8p+16380L : inexact-ok = tgamma towardzero ldbl-128 -0x2p-16384L : -0x8p+16380L : inexact-ok = tgamma upward ldbl-128 -0x2p-16384L : -0x8p+16380L : inexact-ok -tgamma 0x1p-16384 +tgamma 0x1p-16384 spurious-overflow:ldbl-96-intel spurious-overflow:ldbl-96-m68k spurious-overflow:ldbl-128 = tgamma downward flt-32 0x8p-152f : 0xf.fffffp+124f : inexact-ok overflow errno-erange-ok = tgamma tonearest flt-32 0x8p-152f : plus_infty : inexact-ok overflow errno-erange = tgamma towardzero flt-32 0x8p-152f : 0xf.fffffp+124f : inexact-ok overflow errno-erange-ok = tgamma upward flt-32 0x8p-152f : plus_infty : inexact-ok overflow errno-erange -= tgamma downward dbl-64 0x8p-152 : 0x1.fffffffffffffp+148 : inexact-ok -= tgamma tonearest dbl-64 0x8p-152 : 0x2p+148 : inexact-ok -= tgamma towardzero dbl-64 0x8p-152 : 0x1.fffffffffffffp+148 : inexact-ok -= tgamma upward dbl-64 0x8p-152 : 0x2p+148 : inexact-ok -= tgamma downward ldbl-96-intel 0x8p-152L : 0x1.fffffffffffffffep+148L : inexact-ok -= tgamma tonearest ldbl-96-intel 0x8p-152L : 0x2p+148L : inexact-ok -= tgamma towardzero ldbl-96-intel 0x8p-152L : 0x1.fffffffffffffffep+148L : inexact-ok -= tgamma upward ldbl-96-intel 0x8p-152L : 0x2p+148L : inexact-ok -= tgamma downward ldbl-96-m68k 0x8p-152L : 0x1.fffffffffffffffep+148L : inexact-ok -= tgamma tonearest ldbl-96-m68k 0x8p-152L : 0x2p+148L : inexact-ok -= tgamma towardzero ldbl-96-m68k 0x8p-152L : 0x1.fffffffffffffffep+148L : inexact-ok -= tgamma upward ldbl-96-m68k 0x8p-152L : 0x2p+148L : inexact-ok -= tgamma downward ldbl-128 0x8p-152L : 0x1.ffffffffffffffffffffffffffffp+148L : inexact-ok -= tgamma tonearest ldbl-128 0x8p-152L : 0x2p+148L : inexact-ok -= tgamma towardzero ldbl-128 0x8p-152L : 0x1.ffffffffffffffffffffffffffffp+148L : inexact-ok -= tgamma upward ldbl-128 0x8p-152L : 0x2p+148L : inexact-ok -= tgamma downward ldbl-128ibm 0x8p-152L : 0x1.ffffffffffffffffffffffffff8p+148L : inexact-ok -= tgamma tonearest ldbl-128ibm 0x8p-152L : 0x2p+148L : inexact-ok -= tgamma towardzero ldbl-128ibm 0x8p-152L : 0x1.ffffffffffffffffffffffffff8p+148L : inexact-ok -= tgamma upward ldbl-128ibm 0x8p-152L : 0x2p+148L : inexact-ok += tgamma downward dbl-64 0x8p-152 : 0x1.fffffffffffffp+148 : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma tonearest dbl-64 0x8p-152 : 0x2p+148 : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma towardzero dbl-64 0x8p-152 : 0x1.fffffffffffffp+148 : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma upward dbl-64 0x8p-152 : 0x2p+148 : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma downward ldbl-96-intel 0x8p-152L : 0x1.fffffffffffffffep+148L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma tonearest ldbl-96-intel 0x8p-152L : 0x2p+148L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma towardzero ldbl-96-intel 0x8p-152L : 0x1.fffffffffffffffep+148L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma upward ldbl-96-intel 0x8p-152L : 0x2p+148L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma downward ldbl-96-m68k 0x8p-152L : 0x1.fffffffffffffffep+148L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma tonearest ldbl-96-m68k 0x8p-152L : 0x2p+148L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma towardzero ldbl-96-m68k 0x8p-152L : 0x1.fffffffffffffffep+148L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma upward ldbl-96-m68k 0x8p-152L : 0x2p+148L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma downward ldbl-128 0x8p-152L : 0x1.ffffffffffffffffffffffffffffp+148L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma tonearest ldbl-128 0x8p-152L : 0x2p+148L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma towardzero ldbl-128 0x8p-152L : 0x1.ffffffffffffffffffffffffffffp+148L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma upward ldbl-128 0x8p-152L : 0x2p+148L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma downward ldbl-128ibm 0x8p-152L : 0x1.ffffffffffffffffffffffffff8p+148L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma tonearest ldbl-128ibm 0x8p-152L : 0x2p+148L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma towardzero ldbl-128ibm 0x8p-152L : 0x1.ffffffffffffffffffffffffff8p+148L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma upward ldbl-128ibm 0x8p-152L : 0x2p+148L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 = tgamma downward dbl-64 0x4p-1076 : 0xf.ffffffffffff8p+1020 : inexact-ok overflow errno-erange-ok = tgamma tonearest dbl-64 0x4p-1076 : plus_infty : inexact-ok overflow errno-erange = tgamma towardzero dbl-64 0x4p-1076 : 0xf.ffffffffffff8p+1020 : inexact-ok overflow errno-erange-ok = tgamma upward dbl-64 0x4p-1076 : plus_infty : inexact-ok overflow errno-erange -= tgamma downward ldbl-96-intel 0x4p-1076L : 0x3.fffffffffffffffcp+1072L : inexact-ok -= tgamma tonearest ldbl-96-intel 0x4p-1076L : 0x4p+1072L : inexact-ok -= tgamma towardzero ldbl-96-intel 0x4p-1076L : 0x3.fffffffffffffffcp+1072L : inexact-ok -= tgamma upward ldbl-96-intel 0x4p-1076L : 0x4p+1072L : inexact-ok -= tgamma downward ldbl-96-m68k 0x4p-1076L : 0x3.fffffffffffffffcp+1072L : inexact-ok -= tgamma tonearest ldbl-96-m68k 0x4p-1076L : 0x4p+1072L : inexact-ok -= tgamma towardzero ldbl-96-m68k 0x4p-1076L : 0x3.fffffffffffffffcp+1072L : inexact-ok -= tgamma upward ldbl-96-m68k 0x4p-1076L : 0x4p+1072L : inexact-ok -= tgamma downward ldbl-128 0x4p-1076L : 0x3.fffffffffffffffffffffffffffep+1072L : inexact-ok -= tgamma tonearest ldbl-128 0x4p-1076L : 0x4p+1072L : inexact-ok -= tgamma towardzero ldbl-128 0x4p-1076L : 0x3.fffffffffffffffffffffffffffep+1072L : inexact-ok -= tgamma upward ldbl-128 0x4p-1076L : 0x4p+1072L : inexact-ok += tgamma downward ldbl-96-intel 0x4p-1076L : 0x3.fffffffffffffffcp+1072L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma tonearest ldbl-96-intel 0x4p-1076L : 0x4p+1072L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma towardzero ldbl-96-intel 0x4p-1076L : 0x3.fffffffffffffffcp+1072L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma upward ldbl-96-intel 0x4p-1076L : 0x4p+1072L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma downward ldbl-96-m68k 0x4p-1076L : 0x3.fffffffffffffffcp+1072L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma tonearest ldbl-96-m68k 0x4p-1076L : 0x4p+1072L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma towardzero ldbl-96-m68k 0x4p-1076L : 0x3.fffffffffffffffcp+1072L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma upward ldbl-96-m68k 0x4p-1076L : 0x4p+1072L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma downward ldbl-128 0x4p-1076L : 0x3.fffffffffffffffffffffffffffep+1072L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma tonearest ldbl-128 0x4p-1076L : 0x4p+1072L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma towardzero ldbl-128 0x4p-1076L : 0x3.fffffffffffffffffffffffffffep+1072L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 += tgamma upward ldbl-128 0x4p-1076L : 0x4p+1072L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 = tgamma downward ldbl-128ibm 0x4p-1076L : 0xf.ffffffffffffbffffffffffffcp+1020L : inexact-ok overflow errno-erange-ok = tgamma tonearest ldbl-128ibm 0x4p-1076L : plus_infty : inexact-ok overflow errno-erange = tgamma towardzero ldbl-128ibm 0x4p-1076L : 0xf.ffffffffffffbffffffffffffcp+1020L : inexact-ok overflow errno-erange-ok = tgamma upward ldbl-128ibm 0x4p-1076L : plus_infty : inexact-ok overflow errno-erange -= tgamma downward ldbl-96-intel 0x1p-16384L : 0xf.fffffffffffffffp+16380L : inexact-ok += tgamma downward ldbl-96-intel 0x1p-16384L : 0xf.fffffffffffffffp+16380L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 = tgamma tonearest ldbl-96-intel 0x1p-16384L : plus_infty : inexact-ok overflow errno-erange -= tgamma towardzero ldbl-96-intel 0x1p-16384L : 0xf.fffffffffffffffp+16380L : inexact-ok += tgamma towardzero ldbl-96-intel 0x1p-16384L : 0xf.fffffffffffffffp+16380L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 = tgamma upward ldbl-96-intel 0x1p-16384L : plus_infty : inexact-ok overflow errno-erange -= tgamma downward ldbl-96-m68k 0x1p-16384L : 0xf.fffffffffffffffp+16380L : inexact-ok += tgamma downward ldbl-96-m68k 0x1p-16384L : 0xf.fffffffffffffffp+16380L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 = tgamma tonearest ldbl-96-m68k 0x1p-16384L : plus_infty : inexact-ok overflow errno-erange -= tgamma towardzero ldbl-96-m68k 0x1p-16384L : 0xf.fffffffffffffffp+16380L : inexact-ok += tgamma towardzero ldbl-96-m68k 0x1p-16384L : 0xf.fffffffffffffffp+16380L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 = tgamma upward ldbl-96-m68k 0x1p-16384L : plus_infty : inexact-ok overflow errno-erange -= tgamma downward ldbl-128 0x1p-16384L : 0xf.fffffffffffffffffffffffffff8p+16380L : inexact-ok += tgamma downward ldbl-128 0x1p-16384L : 0xf.fffffffffffffffffffffffffff8p+16380L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 = tgamma tonearest ldbl-128 0x1p-16384L : plus_infty : inexact-ok overflow errno-erange -= tgamma towardzero ldbl-128 0x1p-16384L : 0xf.fffffffffffffffffffffffffff8p+16380L : inexact-ok += tgamma towardzero ldbl-128 0x1p-16384L : 0xf.fffffffffffffffffffffffffff8p+16380L : inexact-ok overflow-ok:ldbl-96-intel overflow-ok:ldbl-96-m68k overflow-ok:ldbl-128 = tgamma upward ldbl-128 0x1p-16384L : plus_infty : inexact-ok overflow errno-erange tgamma -0x1p-16384 = tgamma downward flt-32 -0x8p-152f : minus_infty : inexact-ok overflow errno-erange @@ -221443,6 +221443,75 @@ tgamma -100000.5 = tgamma tonearest ldbl-128ibm -0x1.86a08p+16L : -0x0p+0L : inexact-ok underflow errno-erange = tgamma towardzero ldbl-128ibm -0x1.86a08p+16L : -0x0p+0L : inexact-ok underflow errno-erange = tgamma upward ldbl-128ibm -0x1.86a08p+16L : -0x0p+0L : inexact-ok underflow errno-erange +tgamma max += tgamma downward flt-32 0xf.fffffp+124f : 0xf.fffffp+124f : inexact-ok overflow errno-erange-ok += tgamma tonearest flt-32 0xf.fffffp+124f : plus_infty : inexact-ok overflow errno-erange += tgamma towardzero flt-32 0xf.fffffp+124f : 0xf.fffffp+124f : inexact-ok overflow errno-erange-ok += tgamma upward flt-32 0xf.fffffp+124f : plus_infty : inexact-ok overflow errno-erange += tgamma downward dbl-64 0xf.fffffp+124 : 0xf.ffffffffffff8p+1020 : inexact-ok overflow errno-erange-ok += tgamma tonearest dbl-64 0xf.fffffp+124 : plus_infty : inexact-ok overflow errno-erange += tgamma towardzero dbl-64 0xf.fffffp+124 : 0xf.ffffffffffff8p+1020 : inexact-ok overflow errno-erange-ok += tgamma upward dbl-64 0xf.fffffp+124 : plus_infty : inexact-ok overflow errno-erange += tgamma downward ldbl-96-intel 0xf.fffffp+124L : 0xf.fffffffffffffffp+16380L : inexact-ok overflow errno-erange-ok += tgamma tonearest ldbl-96-intel 0xf.fffffp+124L : plus_infty : inexact-ok overflow errno-erange += tgamma towardzero ldbl-96-intel 0xf.fffffp+124L : 0xf.fffffffffffffffp+16380L : inexact-ok overflow errno-erange-ok += tgamma upward ldbl-96-intel 0xf.fffffp+124L : plus_infty : inexact-ok overflow errno-erange += tgamma downward ldbl-96-m68k 0xf.fffffp+124L : 0xf.fffffffffffffffp+16380L : inexact-ok overflow errno-erange-ok += tgamma tonearest ldbl-96-m68k 0xf.fffffp+124L : plus_infty : inexact-ok overflow errno-erange += tgamma towardzero ldbl-96-m68k 0xf.fffffp+124L : 0xf.fffffffffffffffp+16380L : inexact-ok overflow errno-erange-ok += tgamma upward ldbl-96-m68k 0xf.fffffp+124L : plus_infty : inexact-ok overflow errno-erange += tgamma downward ldbl-128 0xf.fffffp+124L : 0xf.fffffffffffffffffffffffffff8p+16380L : inexact-ok overflow errno-erange-ok += tgamma tonearest ldbl-128 0xf.fffffp+124L : plus_infty : inexact-ok overflow errno-erange += tgamma towardzero ldbl-128 0xf.fffffp+124L : 0xf.fffffffffffffffffffffffffff8p+16380L : inexact-ok overflow errno-erange-ok += tgamma upward ldbl-128 0xf.fffffp+124L : plus_infty : inexact-ok overflow errno-erange += tgamma downward ldbl-128ibm 0xf.fffffp+124L : 0xf.ffffffffffffbffffffffffffcp+1020L : inexact-ok overflow errno-erange-ok += tgamma tonearest ldbl-128ibm 0xf.fffffp+124L : plus_infty : inexact-ok overflow errno-erange += tgamma towardzero ldbl-128ibm 0xf.fffffp+124L : 0xf.ffffffffffffbffffffffffffcp+1020L : inexact-ok overflow errno-erange-ok += tgamma upward ldbl-128ibm 0xf.fffffp+124L : plus_infty : inexact-ok overflow errno-erange += tgamma downward dbl-64 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffff8p+1020 : inexact-ok overflow errno-erange-ok += tgamma tonearest dbl-64 0xf.ffffffffffff8p+1020 : plus_infty : inexact-ok overflow errno-erange += tgamma towardzero dbl-64 0xf.ffffffffffff8p+1020 : 0xf.ffffffffffff8p+1020 : inexact-ok overflow errno-erange-ok += tgamma upward dbl-64 0xf.ffffffffffff8p+1020 : plus_infty : inexact-ok overflow errno-erange += tgamma downward ldbl-96-intel 0xf.ffffffffffff8p+1020L : 0xf.fffffffffffffffp+16380L : inexact-ok overflow errno-erange-ok += tgamma tonearest ldbl-96-intel 0xf.ffffffffffff8p+1020L : plus_infty : inexact-ok overflow errno-erange += tgamma towardzero ldbl-96-intel 0xf.ffffffffffff8p+1020L : 0xf.fffffffffffffffp+16380L : inexact-ok overflow errno-erange-ok += tgamma upward ldbl-96-intel 0xf.ffffffffffff8p+1020L : plus_infty : inexact-ok overflow errno-erange += tgamma downward ldbl-96-m68k 0xf.ffffffffffff8p+1020L : 0xf.fffffffffffffffp+16380L : inexact-ok overflow errno-erange-ok += tgamma tonearest ldbl-96-m68k 0xf.ffffffffffff8p+1020L : plus_infty : inexact-ok overflow errno-erange += tgamma towardzero ldbl-96-m68k 0xf.ffffffffffff8p+1020L : 0xf.fffffffffffffffp+16380L : inexact-ok overflow errno-erange-ok += tgamma upward ldbl-96-m68k 0xf.ffffffffffff8p+1020L : plus_infty : inexact-ok overflow errno-erange += tgamma downward ldbl-128 0xf.ffffffffffff8p+1020L : 0xf.fffffffffffffffffffffffffff8p+16380L : inexact-ok overflow errno-erange-ok += tgamma tonearest ldbl-128 0xf.ffffffffffff8p+1020L : plus_infty : inexact-ok overflow errno-erange += tgamma towardzero ldbl-128 0xf.ffffffffffff8p+1020L : 0xf.fffffffffffffffffffffffffff8p+16380L : inexact-ok overflow errno-erange-ok += tgamma upward ldbl-128 0xf.ffffffffffff8p+1020L : plus_infty : inexact-ok overflow errno-erange += tgamma downward ldbl-128ibm 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffffbffffffffffffcp+1020L : inexact-ok overflow errno-erange-ok += tgamma tonearest ldbl-128ibm 0xf.ffffffffffff8p+1020L : plus_infty : inexact-ok overflow errno-erange += tgamma towardzero ldbl-128ibm 0xf.ffffffffffff8p+1020L : 0xf.ffffffffffffbffffffffffffcp+1020L : inexact-ok overflow errno-erange-ok += tgamma upward ldbl-128ibm 0xf.ffffffffffff8p+1020L : plus_infty : inexact-ok overflow errno-erange += tgamma downward ldbl-96-intel 0xf.fffffffffffffffp+16380L : 0xf.fffffffffffffffp+16380L : inexact-ok overflow errno-erange-ok += tgamma tonearest ldbl-96-intel 0xf.fffffffffffffffp+16380L : plus_infty : inexact-ok overflow errno-erange += tgamma towardzero ldbl-96-intel 0xf.fffffffffffffffp+16380L : 0xf.fffffffffffffffp+16380L : inexact-ok overflow errno-erange-ok += tgamma upward ldbl-96-intel 0xf.fffffffffffffffp+16380L : plus_infty : inexact-ok overflow errno-erange += tgamma downward ldbl-96-m68k 0xf.fffffffffffffffp+16380L : 0xf.fffffffffffffffp+16380L : inexact-ok overflow errno-erange-ok += tgamma tonearest ldbl-96-m68k 0xf.fffffffffffffffp+16380L : plus_infty : inexact-ok overflow errno-erange += tgamma towardzero ldbl-96-m68k 0xf.fffffffffffffffp+16380L : 0xf.fffffffffffffffp+16380L : inexact-ok overflow errno-erange-ok += tgamma upward ldbl-96-m68k 0xf.fffffffffffffffp+16380L : plus_infty : inexact-ok overflow errno-erange += tgamma downward ldbl-128 0xf.fffffffffffffffp+16380L : 0xf.fffffffffffffffffffffffffff8p+16380L : inexact-ok overflow errno-erange-ok += tgamma tonearest ldbl-128 0xf.fffffffffffffffp+16380L : plus_infty : inexact-ok overflow errno-erange += tgamma towardzero ldbl-128 0xf.fffffffffffffffp+16380L : 0xf.fffffffffffffffffffffffffff8p+16380L : inexact-ok overflow errno-erange-ok += tgamma upward ldbl-128 0xf.fffffffffffffffp+16380L : plus_infty : inexact-ok overflow errno-erange += tgamma downward ldbl-128 0xf.fffffffffffffffffffffffffff8p+16380L : 0xf.fffffffffffffffffffffffffff8p+16380L : inexact-ok overflow errno-erange-ok += tgamma tonearest ldbl-128 0xf.fffffffffffffffffffffffffff8p+16380L : plus_infty : inexact-ok overflow errno-erange += tgamma towardzero ldbl-128 0xf.fffffffffffffffffffffffffff8p+16380L : 0xf.fffffffffffffffffffffffffff8p+16380L : inexact-ok overflow errno-erange-ok += tgamma upward ldbl-128 0xf.fffffffffffffffffffffffffff8p+16380L : plus_infty : inexact-ok overflow errno-erange += tgamma downward ldbl-128 0xf.ffffffffffffbffffffffffffcp+1020L : 0xf.fffffffffffffffffffffffffff8p+16380L : inexact-ok overflow errno-erange-ok += tgamma tonearest ldbl-128 0xf.ffffffffffffbffffffffffffcp+1020L : plus_infty : inexact-ok overflow errno-erange += tgamma towardzero ldbl-128 0xf.ffffffffffffbffffffffffffcp+1020L : 0xf.fffffffffffffffffffffffffff8p+16380L : inexact-ok overflow errno-erange-ok += tgamma upward ldbl-128 0xf.ffffffffffffbffffffffffffcp+1020L : plus_infty : inexact-ok overflow errno-erange += tgamma downward ldbl-128ibm 0xf.ffffffffffffbffffffffffffcp+1020L : 0xf.ffffffffffffbffffffffffffcp+1020L : inexact-ok overflow errno-erange-ok += tgamma tonearest ldbl-128ibm 0xf.ffffffffffffbffffffffffffcp+1020L : plus_infty : inexact-ok overflow errno-erange += tgamma towardzero ldbl-128ibm 0xf.ffffffffffffbffffffffffffcp+1020L : 0xf.ffffffffffffbffffffffffffcp+1020L : inexact-ok overflow errno-erange-ok += tgamma upward ldbl-128ibm 0xf.ffffffffffffbffffffffffffcp+1020L : plus_infty : inexact-ok overflow errno-erange tgamma -0x3.06644cp+0 = tgamma downward flt-32 -0x3.06644cp+0f : 0x6.79a69p+0f : inexact-ok = tgamma tonearest flt-32 -0x3.06644cp+0f : 0x6.79a698p+0f : inexact-ok diff --git a/math/gen-libm-have-vector-test.sh b/math/gen-libm-have-vector-test.sh index 529ea0c1ac..6532deb249 100644 --- a/math/gen-libm-have-vector-test.sh +++ b/math/gen-libm-have-vector-test.sh @@ -50,11 +50,3 @@ for func in $(cat libm-test.inc | grep ALL_RM_TEST | grep RUN_TEST_LOOP_fFF_11 | print_defs ${func}f "_fFF" print_defs ${func}l "_fFF" done - -# When all functions will use ALL_RM_TEST instead of using START directly, -# this code can be removed. -for func in $(grep 'START.*;$' libm-test.inc | sed -r "s/.*\(//; s/,.*//"); do - print_defs ${func} - print_defs ${func}f - print_defs ${func}l -done diff --git a/math/libm-test.inc b/math/libm-test.inc index 9e402ab634..17a20c45ad 100644 --- a/math/libm-test.inc +++ b/math/libm-test.inc @@ -9484,7 +9484,6 @@ tanh_test (void) static const struct test_f_f_data tgamma_test_data[] = { TEST_f_f (tgamma, plus_infty, plus_infty), - TEST_f_f (tgamma, max_value, plus_infty, OVERFLOW_EXCEPTION|ERRNO_ERANGE), TEST_f_f (tgamma, 0, plus_infty, DIVIDE_BY_ZERO_EXCEPTION|ERRNO_ERANGE), TEST_f_f (tgamma, minus_zero, minus_infty, DIVIDE_BY_ZERO_EXCEPTION|ERRNO_ERANGE), /* tgamma (x) == qNaN plus invalid exception for integer x <= 0. */ @@ -9499,9 +9498,7 @@ static const struct test_f_f_data tgamma_test_data[] = static void tgamma_test (void) { - START (tgamma,, 0); - RUN_TEST_LOOP_f_f (tgamma, tgamma_test_data, ); - END; + ALL_RM_TEST (tgamma, 0, tgamma_test_data, RUN_TEST_LOOP_f_f, END); } -- cgit 1.4.1