aboutsummaryrefslogtreecommitdiff
path: root/sysdeps/ieee754
diff options
context:
space:
mode:
authorUlrich Drepper <drepper@redhat.com>1999-10-20 03:20:31 +0000
committerUlrich Drepper <drepper@redhat.com>1999-10-20 03:20:31 +0000
commitf30e0cd35eae3bfac07df89db16c04de1a39671d (patch)
tree60b216a7b4e062961d9774dfc57b80f3362b1974 /sysdeps/ieee754
parentfe559c5e1ce629f94f1621c76f5cfd59d6028830 (diff)
downloadglibc-f30e0cd35eae3bfac07df89db16c04de1a39671d.zip
glibc-f30e0cd35eae3bfac07df89db16c04de1a39671d.tar.gz
glibc-f30e0cd35eae3bfac07df89db16c04de1a39671d.tar.bz2
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')
-rw-r--r--sysdeps/ieee754/flt-32/e_gammaf_r.c3
-rw-r--r--sysdeps/ieee754/flt-32/e_lgammaf_r.c2
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 926905e..f0cd438 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 f744d53..5ebebb7 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);