about summary refs log tree commit diff
path: root/sysdeps/ieee754/ldbl-96/s_scalbnl.c
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/ldbl-96/s_scalbnl.c
parent1474b1c6c55b3d647e1f3ccc89821f8a35124dcb (diff)
downloadglibc-4c49989e4178fc032cfdf37ee898babc7e3fa593.tar.gz
glibc-4c49989e4178fc032cfdf37ee898babc7e3fa593.tar.xz
glibc-4c49989e4178fc032cfdf37ee898babc7e3fa593.zip
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/ldbl-96/s_scalbnl.c')
-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 34c52e773a..9f441fd1e9 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)