From f30e0cd35eae3bfac07df89db16c04de1a39671d Mon Sep 17 00:00:00 2001 From: Ulrich Drepper Date: Wed, 20 Oct 1999 03:20:31 +0000 Subject: Update. 1999-10-19 Ulrich Drepper * sysdeps/i386/fpu/s_nextafterl.c: Add __nextafterl and nextafterl aliases. * sysdeps/ieee754/flt-32/e_lgammaf_r.c: Don't handle -Inf special. * sysdeps/ieee754/flt-32/e_gammaf_r.c (__ieee754_gammaf_r): Check for -Inf and return NaN. * math/gen-libm-test.pl: Fix program name in help message. * math/libm-test.inc (check_complex): It's Imaginary, not Complex. * math/libm-test.inc (gamma_test): Result of gamma(-inf) is +inf. * sysdeps/i386/Implies: Correct order of libm directories. --- sysdeps/ieee754/flt-32/e_gammaf_r.c | 3 +++ sysdeps/ieee754/flt-32/e_lgammaf_r.c | 2 -- 2 files changed, 3 insertions(+), 2 deletions(-) (limited to 'sysdeps/ieee754/flt-32') diff --git a/sysdeps/ieee754/flt-32/e_gammaf_r.c b/sysdeps/ieee754/flt-32/e_gammaf_r.c index 926905e7cd..f0cd43840b 100644 --- a/sysdeps/ieee754/flt-32/e_gammaf_r.c +++ b/sysdeps/ieee754/flt-32/e_gammaf_r.c @@ -44,6 +44,9 @@ __ieee754_gammaf_r (float x, int *signgamp) *signgamp = 0; return (x - x) / (x - x); } + if (hx == 0xff800000) + /* x == -Inf. According to ISO this is NaN. */ + return x - x; /* XXX FIXME. */ return __ieee754_expf (__ieee754_lgammaf_r (x, signgamp)); diff --git a/sysdeps/ieee754/flt-32/e_lgammaf_r.c b/sysdeps/ieee754/flt-32/e_lgammaf_r.c index f744d5320e..5ebebb77af 100644 --- a/sysdeps/ieee754/flt-32/e_lgammaf_r.c +++ b/sysdeps/ieee754/flt-32/e_lgammaf_r.c @@ -162,8 +162,6 @@ static float zero= 0.0000000000e+00; /* purge off +-inf, NaN, +-0, and negative arguments */ *signgamp = 1; - if ((unsigned int)hx==0xff800000) - return x-x; ix = hx&0x7fffffff; if(ix>=0x7f800000) return x*x; if(ix==0) return one/fabsf(x); -- cgit v1.2.3-70-g09d2