diff options
Diffstat (limited to 'sysdeps/ieee754/ldbl-96')
| -rw-r--r-- | sysdeps/ieee754/ldbl-96/e_gammal_r.c | 15 |
1 files changed, 9 insertions, 6 deletions
diff --git a/sysdeps/ieee754/ldbl-96/e_gammal_r.c b/sysdeps/ieee754/ldbl-96/e_gammal_r.c index 376d2b9d94..9b8cad8e6d 100644 --- a/sysdeps/ieee754/ldbl-96/e_gammal_r.c +++ b/sysdeps/ieee754/ldbl-96/e_gammal_r.c @@ -38,18 +38,21 @@ __ieee754_gammal_r (long double x, int *signgamp) *signgamp = 0; return x / x; } - if ((es & 0x8000) != 0 && x < 0xffffffff && __rintl (x) == x) - { - /* Return value for integer x < 0 is NaN with invalid exception. */ - *signgamp = 0; - return (x - x) / (x - x); - } if (es == 0xffffffff && ((hx & 0x7fffffff) | lx) == 0) { /* x == -Inf. According to ISO this is NaN. */ *signgamp = 0; return x - x; } + if ((es & 0x7fff) == 0x7fff && ((hx & 0x7fffffff) | lx) != 0) + /* NaN, return it. */ + return x; + if ((es & 0x8000) != 0 && x < 0xffffffff && __rintl (x) == x) + { + /* Return value for integer x < 0 is NaN with invalid exception. */ + *signgamp = 0; + return (x - x) / (x - x); + } /* XXX FIXME. */ return __ieee754_expl (__ieee754_lgammal_r (x, signgamp)); |