From e47cc4e0ed06f4167b6bd8ac39e12d094e0dc474 Mon Sep 17 00:00:00 2001 From: Joseph Myers Date: Thu, 5 Dec 2013 14:01:41 +0000 Subject: Fix tgamma errno setting on underflow (bug 6810). --- math/w_tgammaf.c | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) (limited to 'math/w_tgammaf.c') diff --git a/math/w_tgammaf.c b/math/w_tgammaf.c index 8bb553e..fb15b6a 100644 --- a/math/w_tgammaf.c +++ b/math/w_tgammaf.c @@ -13,6 +13,7 @@ * ==================================================== */ +#include #include #include @@ -22,7 +23,7 @@ __tgammaf(float x) int local_signgam; float y = __ieee754_gammaf_r(x,&local_signgam); - if(__builtin_expect(!__finitef(y), 0) + if(__glibc_unlikely (!__finitef (y) || y == 0) && (__finitef (x) || __isinff (x) < 0) && _LIB_VERSION != _IEEE_) { if (x == (float)0.0) @@ -31,6 +32,9 @@ __tgammaf(float x) else if(__floorf(x)==x&&x<0.0f) /* tgammaf domain */ return __kernel_standard_f(x, x, 141); + else if (y == 0) + /* tgammaf underflow */ + __set_errno (ERANGE); else /* tgammaf overflow */ return __kernel_standard_f(x, x, 140); -- cgit v1.1