aboutsummaryrefslogtreecommitdiff
path: root/sysdeps/ieee754
diff options
context:
space:
mode:
authorUlrich Drepper <drepper@redhat.com>2002-06-20 07:23:23 +0000
committerUlrich Drepper <drepper@redhat.com>2002-06-20 07:23:23 +0000
commit4c49989e4178fc032cfdf37ee898babc7e3fa593 (patch)
tree24f6855b68683881dfed04a4404a0e93a2b7fc19 /sysdeps/ieee754
parent1474b1c6c55b3d647e1f3ccc89821f8a35124dcb (diff)
downloadglibc-4c49989e4178fc032cfdf37ee898babc7e3fa593.zip
glibc-4c49989e4178fc032cfdf37ee898babc7e3fa593.tar.gz
glibc-4c49989e4178fc032cfdf37ee898babc7e3fa593.tar.bz2
Update.
2002-06-13 Stephen L Moshier <steve@moshier.net> * sysdeps/ieee754/ldbl-96/s_scalbnl.c: Fix cases in which argument or result is subnormal. 2002-06-19 Ulrich Drepper <drepper@redhat.com> * sysdeps/i386/fpu/s_remquo.S: Fix three mistakes in computing the quotient. * sysdeps/i386/fpu/s_remquof.S: Likewise. * sysdeps/i386/fpu/s_remquol.S: Likewise. * sysdeps/hppa/dl-machine.h: Don't use multi-line string literals. * sysdeps/unix/sysv/linux/alpha/sendfile64.c: New file. * sysdeps/unix/sysv/linux/ia64/sendfile64.c: New file. * sysdeps/unix/sysv/linux/sparc/sparc64/sendfile64.c: New file. * sysdeps/unix/sysv/linux/x86_64/sendfile64.c: New file.
Diffstat (limited to 'sysdeps/ieee754')
-rw-r--r--sysdeps/ieee754/ldbl-96/s_scalbnl.c16
1 files changed, 8 insertions, 8 deletions
diff --git a/sysdeps/ieee754/ldbl-96/s_scalbnl.c b/sysdeps/ieee754/ldbl-96/s_scalbnl.c
index 34c52e7..9f441fd 100644
--- a/sysdeps/ieee754/ldbl-96/s_scalbnl.c
+++ b/sysdeps/ieee754/ldbl-96/s_scalbnl.c
@@ -33,8 +33,8 @@ static const long double
#else
static long double
#endif
-two63 = 4.50359962737049600000e+15,
-twom63 = 1.08420217248550443400e-19,
+two64 = 1.8446744073709551616e19L,
+twom64 = 5.421010862427522170037e-20L,
huge = 1.0e+4900L,
tiny = 1.0e-4900L;
@@ -50,9 +50,9 @@ tiny = 1.0e-4900L;
k = es&0x7fff; /* extract exponent */
if (k==0) { /* 0 or subnormal x */
if ((lx|(hx&0x7fffffff))==0) return x; /* +-0 */
- x *= two63;
- GET_LDOUBLE_EXP(es,x);
- k = (hx&0x7fff) - 63;
+ x *= two64;
+ GET_LDOUBLE_EXP(hx,x);
+ k = (hx&0x7fff) - 64;
}
if (k==0x7fff) return x+x; /* NaN or Inf */
k = k+n;
@@ -62,10 +62,10 @@ tiny = 1.0e-4900L;
return tiny*__copysignl(tiny,x);
if (k > 0) /* normal result */
{SET_LDOUBLE_EXP(x,(es&0x8000)|k); return x;}
- if (k <= -63)
+ if (k <= -64)
return tiny*__copysignl(tiny,x); /*underflow*/
- k += 63; /* subnormal result */
+ k += 64; /* subnormal result */
SET_LDOUBLE_EXP(x,(es&0x8000)|k);
- return x*twom63;
+ return x*twom64;
}
weak_alias (__scalbnl, scalbnl)