diff options
author | Joseph Myers <joseph@codesourcery.com> | 2013-12-05 14:01:41 +0000 |
---|---|---|
committer | Joseph Myers <joseph@codesourcery.com> | 2013-12-05 14:01:41 +0000 |
commit | e47cc4e0ed06f4167b6bd8ac39e12d094e0dc474 (patch) | |
tree | 1a1f08bd5751795b08e013f499f49df0830917fc /math/w_tgamma.c | |
parent | bbf37bdc12c0370a0ac2fa597cddd660bad6ec5e (diff) | |
download | glibc-e47cc4e0ed06f4167b6bd8ac39e12d094e0dc474.zip glibc-e47cc4e0ed06f4167b6bd8ac39e12d094e0dc474.tar.gz glibc-e47cc4e0ed06f4167b6bd8ac39e12d094e0dc474.tar.bz2 |
Fix tgamma errno setting on underflow (bug 6810).
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 */ } |