aboutsummaryrefslogtreecommitdiff
path: root/sysdeps/ieee754/ldbl-128
diff options
context:
space:
mode:
authorJoseph Myers <joseph@codesourcery.com>2014-01-06 18:20:20 +0000
committerJoseph Myers <joseph@codesourcery.com>2014-01-06 18:20:20 +0000
commiteb3fc44b56e5e780be5e2830c72d20b9e74fef8a (patch)
tree77225ddae4384596a6c9e46ae8340ee9e2698032 /sysdeps/ieee754/ldbl-128
parent5fd3eccbea76216091e147ebd84ed64a839f76eb (diff)
downloadglibc-eb3fc44b56e5e780be5e2830c72d20b9e74fef8a.tar
glibc-eb3fc44b56e5e780be5e2830c72d20b9e74fef8a.tar.gz
glibc-eb3fc44b56e5e780be5e2830c72d20b9e74fef8a.tar.bz2
glibc-eb3fc44b56e5e780be5e2830c72d20b9e74fef8a.zip
Fix ldbl-128 / ldbl-128ibm lgammal spurious underflow (bug 16400).
This patch fixes bug 16400, spurious underflow exceptions for ldbl-128 / ldbl-128ibm lgammal with small positive arguments, by just using -__logl (x) as the result in the problem cases (similar to the previous fix for problems with small negative arguments). Tested powerpc32, and also tested on mips64 that this does not require ulps regeneration for the ldbl-128 case. * sysdeps/ieee754/ldbl-128/e_lgammal_r.c (__ieee754_lgammal_r): Return -__logl (x) for small positive arguments without evaluating a polynomial.
Diffstat (limited to 'sysdeps/ieee754/ldbl-128')
-rw-r--r--sysdeps/ieee754/ldbl-128/e_lgammal_r.c4
1 files changed, 3 insertions, 1 deletions
diff --git a/sysdeps/ieee754/ldbl-128/e_lgammal_r.c b/sysdeps/ieee754/ldbl-128/e_lgammal_r.c
index 23ab9b9a43..1961355a73 100644
--- a/sysdeps/ieee754/ldbl-128/e_lgammal_r.c
+++ b/sysdeps/ieee754/ldbl-128/e_lgammal_r.c
@@ -805,7 +805,9 @@ __ieee754_lgammal_r (long double x, int *signgamp)
{
case 0:
/* log gamma (x + 1) = log(x) + log gamma(x) */
- if (x <= 0.125)
+ if (x < 0x1p-120L)
+ return -__logl (x);
+ else if (x <= 0.125)
{
p = x * neval (x, RN1, NRN1) / deval (x, RD1, NRD1);
}