Fix ldbl-128 / ldbl-128ibm lgammal spurious underflow (bug 16400).

[thirdparty/glibc.git] / sysdeps / ieee754 / ldbl-128 / e_lgammal_r.c

diff --git a/sysdeps/ieee754/ldbl-128/e_lgammal_r.c b/sysdeps/ieee754/ldbl-128/e_lgammal_r.c
index 23ab9b9a43977259c8f1c6f002a8c8fd347b7c38..1961355a73e4016db104a1b172a0e81db2d19930 100644 (file)
--- 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);
            }