diff options
| author | Andreas Jaeger <aj@suse.de> | 2002-01-07 13:28:49 +0000 |
|---|---|---|
| committer | Andreas Jaeger <aj@suse.de> | 2002-01-07 13:28:49 +0000 |
| commit | c195dcdd754977eabf32c03bd8f2c764a00efc59 (patch) | |
| tree | 3b40dd04e7127fcd5bda68f78063109f6f59a53f /sysdeps | |
| parent | af0c938d7b0ddc982fb6c81bdef52d243812edfa (diff) | |
| download | glibc-c195dcdd754977eabf32c03bd8f2c764a00efc59.zip glibc-c195dcdd754977eabf32c03bd8f2c764a00efc59.tar.gz glibc-c195dcdd754977eabf32c03bd8f2c764a00efc59.tar.bz2 | |
Update.
2002-01-07 Stephen L Moshier <moshier@mediaone.net>
* sysdeps/ieee754/ldbl-96/s_erfl.c (erfcl): Fix K&R header.
* sysdeps/ieee754/ldbl-96/e_lgammal_r.c (sin_pi): Fix typo in test
for x < 0.25 and restore original range reduction method.
(__ieee754_lgammal_r): Make sure signgam is set before returning.
Diffstat (limited to 'sysdeps')
| -rw-r--r-- | sysdeps/ieee754/ldbl-96/e_lgammal_r.c | 26 |
1 files changed, 19 insertions, 7 deletions
diff --git a/sysdeps/ieee754/ldbl-96/e_lgammal_r.c b/sysdeps/ieee754/ldbl-96/e_lgammal_r.c index f39ef35..ee051c9 100644 --- a/sysdeps/ieee754/ldbl-96/e_lgammal_r.c +++ b/sysdeps/ieee754/ldbl-96/e_lgammal_r.c @@ -84,6 +84,7 @@ static long double half = 0.5L, one = 1.0L, pi = 3.14159265358979323846264L, + two63 = 9.223372036854775808e18L, /* lgam(1+x) = 0.5 x + x a(x)/b(x) -0.268402099609375 <= x <= 0 @@ -206,8 +207,7 @@ sin_pi (x) GET_LDOUBLE_WORDS (se, i0, i1, x); ix = se & 0x7fff; - - i1 = (ix << 16) | (i0 >> 16); + ix = (ix << 16) | (i0 >> 16); if (ix < 0x3ffd8000) /* 0.25 */ return __sinl (pi * x); y = -x; /* x is assume negative */ @@ -219,13 +219,25 @@ sin_pi (x) z = __floorl (y); if (z != y) { /* inexact anyway */ - y *= half; - y = 2.0 * (y - __floorl (y)); /* y = |x| mod 2.0 */ - n = (int) (y * 4.0); + y *= 0.5; + y = 2.0*(y - __floorl(y)); /* y = |x| mod 2.0 */ + n = (int) (y*4.0); } else { - return (zero + zero); + if (ix >= 0x403f8000) /* 2^64 */ + { + y = zero; n = 0; /* y must be even */ + } + else + { + if (ix < 0x403e8000) /* 2^63 */ + z = y + two63; /* exact */ + GET_LDOUBLE_WORDS (se, i0, i1, z); + n = i1 & 1; + y = n; + n <<= 2; + } } switch (n) @@ -267,6 +279,7 @@ __ieee754_lgammal_r (x, signgamp) int i, ix; u_int32_t se, i0, i1; + *signgamp = 1; GET_LDOUBLE_WORDS (se, i0, i1, x); ix = se & 0x7fff; @@ -276,7 +289,6 @@ __ieee754_lgammal_r (x, signgamp) ix = (ix << 16) | (i0 >> 16); /* purge off +-inf, NaN, +-0, and negative arguments */ - *signgamp = 1; if (ix >= 0x7fff0000) return x * x; |