diff options
Diffstat (limited to 'math')
| -rw-r--r-- | math/libm-test.inc | 6 | ||||
| -rw-r--r-- | math/w_tgamma.c | 3 | ||||
| -rw-r--r-- | math/w_tgammaf.c | 3 | ||||
| -rw-r--r-- | math/w_tgammal.c | 3 |
4 files changed, 9 insertions, 6 deletions
diff --git a/math/libm-test.inc b/math/libm-test.inc index c508af6045..1ff59e0f97 100644 --- a/math/libm-test.inc +++ b/math/libm-test.inc @@ -13615,9 +13615,9 @@ tgamma_test (void) TEST_f_f (tgamma, 0, plus_infty, DIVIDE_BY_ZERO_EXCEPTION); TEST_f_f (tgamma, minus_zero, minus_infty, DIVIDE_BY_ZERO_EXCEPTION); /* tgamma (x) == qNaN plus invalid exception for integer x <= 0. */ - TEST_f_f (tgamma, -2, qnan_value, INVALID_EXCEPTION); - TEST_f_f (tgamma, -max_value, qnan_value, INVALID_EXCEPTION); - TEST_f_f (tgamma, minus_infty, qnan_value, INVALID_EXCEPTION); + TEST_f_f (tgamma, -2, qnan_value, INVALID_EXCEPTION|ERRNO_EDOM); + TEST_f_f (tgamma, -max_value, qnan_value, INVALID_EXCEPTION|ERRNO_EDOM); + TEST_f_f (tgamma, minus_infty, qnan_value, INVALID_EXCEPTION|ERRNO_EDOM); TEST_f_f (tgamma, qnan_value, qnan_value); TEST_f_f (tgamma, 0.5, M_SQRT_PIl); diff --git a/math/w_tgamma.c b/math/w_tgamma.c index 976b5fb95c..6b6c7c5241 100644 --- a/math/w_tgamma.c +++ b/math/w_tgamma.c @@ -24,7 +24,8 @@ __tgamma(double x) int local_signgam; double y = __ieee754_gamma_r(x,&local_signgam); - if(__builtin_expect(!__finite(y), 0)&&__finite(x) + if(__builtin_expect(!__finite(y), 0) + && (__finite (x) || __isinf (x) < 0) && _LIB_VERSION != _IEEE_) { if (x == 0.0) return __kernel_standard(x,x,50); /* tgamma pole */ diff --git a/math/w_tgammaf.c b/math/w_tgammaf.c index 48141354e2..8bb553e47b 100644 --- a/math/w_tgammaf.c +++ b/math/w_tgammaf.c @@ -22,7 +22,8 @@ __tgammaf(float x) int local_signgam; float y = __ieee754_gammaf_r(x,&local_signgam); - if(__builtin_expect(!__finitef(y), 0) && __finitef(x) + if(__builtin_expect(!__finitef(y), 0) + && (__finitef (x) || __isinff (x) < 0) && _LIB_VERSION != _IEEE_) { if (x == (float)0.0) /* tgammaf pole */ diff --git a/math/w_tgammal.c b/math/w_tgammal.c index 86adab29ea..72b247db2a 100644 --- a/math/w_tgammal.c +++ b/math/w_tgammal.c @@ -27,7 +27,8 @@ __tgammal(long double x) int local_signgam; long double y = __ieee754_gammal_r(x,&local_signgam); - if(__builtin_expect(!__finitel(y), 0) && __finitel(x) + if(__builtin_expect(!__finitel(y), 0) + && (__finitel (x) || __isinfl (x) < 0) && _LIB_VERSION != _IEEE_) { if(x==0.0) return __kernel_standard_l(x,x,250); /* tgamma pole */ |