diff options
Diffstat (limited to 'math/w_tgamma.c')
| -rw-r--r-- | math/w_tgamma.c | 5 |
1 files changed, 4 insertions, 1 deletions
diff --git a/math/w_tgamma.c b/math/w_tgamma.c index 6b6c7c5..504996e 100644 --- a/math/w_tgamma.c +++ b/math/w_tgamma.c @@ -15,6 +15,7 @@ * depending on the library mode. */ +#include <errno.h> #include <math.h> #include <math_private.h> @@ -24,13 +25,15 @@ __tgamma(double x) int local_signgam; double y = __ieee754_gamma_r(x,&local_signgam); - if(__builtin_expect(!__finite(y), 0) + if(__glibc_unlikely (!__finite (y) || y == 0) && (__finite (x) || __isinf (x) < 0) && _LIB_VERSION != _IEEE_) { if (x == 0.0) return __kernel_standard(x,x,50); /* tgamma pole */ else if(__floor(x)==x&&x<0.0) return __kernel_standard(x,x,41); /* tgamma domain */ + else if (y == 0) + __set_errno (ERANGE); /* tgamma underflow */ else return __kernel_standard(x,x,40); /* tgamma overflow */ } |