aboutsummaryrefslogtreecommitdiff
path: root/sysdeps/ieee754
diff options
context:
space:
mode:
authorAndreas Jaeger <aj@suse.de>2002-01-07 13:28:49 +0000
committerAndreas Jaeger <aj@suse.de>2002-01-07 13:28:49 +0000
commitc195dcdd754977eabf32c03bd8f2c764a00efc59 (patch)
tree3b40dd04e7127fcd5bda68f78063109f6f59a53f /sysdeps/ieee754
parentaf0c938d7b0ddc982fb6c81bdef52d243812edfa (diff)
downloadglibc-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/ieee754')
-rw-r--r--sysdeps/ieee754/ldbl-96/e_lgammal_r.c26
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;