From eb3bf57345cb1451a5d2578dc3f7fbe5aff3b8b5 Mon Sep 17 00:00:00 2001 From: Ulrich Drepper Date: Fri, 22 Oct 1999 15:30:17 +0000 Subject: Update. * sysdeps/ieee754/ldbl-96/s_nexttoward.c: If x == y, return y. * sysdeps/ieee754/dbl-64/e_lgamma_r.c: Don't handle -Inf special. * sysdeps/ieee754/dbl-64/e_gamma_r.c (__ieee754_gamma_r): Check for -Inf and return NaN. --- sysdeps/ieee754/dbl-64/e_gamma_r.c | 3 +++ 1 file changed, 3 insertions(+) (limited to 'sysdeps/ieee754/dbl-64/e_gamma_r.c') diff --git a/sysdeps/ieee754/dbl-64/e_gamma_r.c b/sysdeps/ieee754/dbl-64/e_gamma_r.c index bd802c24f1..eff40463be 100644 --- a/sysdeps/ieee754/dbl-64/e_gamma_r.c +++ b/sysdeps/ieee754/dbl-64/e_gamma_r.c @@ -45,6 +45,9 @@ __ieee754_gamma_r (double x, int *signgamp) *signgamp = 0; return (x - x) / (x - x); } + if ((unsigned int) hx == 0xfff00000 && lx==0) + /* x == -Inf. According to ISO this is NaN. */ + return x - x; /* XXX FIXME. */ return __ieee754_exp (__ieee754_lgamma_r (x, signgamp)); -- cgit v1.2.3