diff options
| author | Ulrich Drepper <drepper@redhat.com> | 1999-10-20 03:20:31 +0000 |
|---|---|---|
| committer | Ulrich Drepper <drepper@redhat.com> | 1999-10-20 03:20:31 +0000 |
| commit | f30e0cd35eae3bfac07df89db16c04de1a39671d (patch) | |
| tree | 60b216a7b4e062961d9774dfc57b80f3362b1974 /sysdeps/ieee754/flt-32 | |
| parent | fe559c5e1ce629f94f1621c76f5cfd59d6028830 (diff) | |
| download | glibc-f30e0cd35eae3bfac07df89db16c04de1a39671d.tar.gz glibc-f30e0cd35eae3bfac07df89db16c04de1a39671d.tar.xz glibc-f30e0cd35eae3bfac07df89db16c04de1a39671d.zip | |
Update.
1999-10-19 Ulrich Drepper <drepper@cygnus.com> * 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.
Diffstat (limited to 'sysdeps/ieee754/flt-32')
| -rw-r--r-- | sysdeps/ieee754/flt-32/e_gammaf_r.c | 3 | ||||
| -rw-r--r-- | sysdeps/ieee754/flt-32/e_lgammaf_r.c | 2 |
2 files changed, 3 insertions, 2 deletions
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); |