From 276ae1f29c34cfbc579e7959e8f0bc65916672c0 Mon Sep 17 00:00:00 2001 From: Ulrich Drepper Date: Sat, 17 Mar 2001 03:41:58 +0000 Subject: Update. * sysdeps/ieee754/ldbl-96/e_gammal_r.c (__ieee754_gammal_r): Handle NaN here. --- sysdeps/ieee754/ldbl-96/e_gammal_r.c | 15 +++++++++------ 1 file changed, 9 insertions(+), 6 deletions(-) (limited to 'sysdeps/ieee754/ldbl-96') 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)); -- cgit 1.4.1