about summary refs log tree commit diff
diff options
context:
space:
mode:
-rw-r--r--ChangeLog22
-rw-r--r--sysdeps/i386/fpu/s_remquo.S5
-rw-r--r--sysdeps/i386/fpu/s_remquof.S5
-rw-r--r--sysdeps/i386/fpu/s_remquol.S5
-rw-r--r--sysdeps/ieee754/ldbl-96/s_scalbnl.c16
5 files changed, 35 insertions, 18 deletions
diff --git a/ChangeLog b/ChangeLog
index 077cfd5d89..d1eb75e2e3 100644
--- a/ChangeLog
+++ b/ChangeLog
@@ -1,3 +1,17 @@
+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.
+
 2002-06-17  David Mosberger  <davidm@hpl.hp.com>
 
 	* sysdeps/unix/sysv/linux/ia64/sys/ucontext.h (uc_sigmask): Fix
@@ -142,10 +156,10 @@
 
 2002-06-07  Jakub Jelinek  <jakub@redhat.com>
 
-	* sysdeps/unix/sysv/linux/alpha/sendfile64.c: New.
-	* sysdeps/unix/sysv/linux/ia64/sendfile64.c: New.
-	* sysdeps/unix/sysv/linux/sparc/sparc64/sendfile64.c: New.
-	* sysdeps/unix/sysv/linux/x86_64/sendfile64.c: New.
+	* 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.
 	* sysdeps/unix/sysv/linux/syscalls.list (sendfile): Remove EXTRA.
 	(sendfile64): Likewise.
 	* sysdeps/unix/sysv/linux/alpha/syscalls.list (sendfile): Likewise.
diff --git a/sysdeps/i386/fpu/s_remquo.S b/sysdeps/i386/fpu/s_remquo.S
index 4857b98ab4..5056593214 100644
--- a/sysdeps/i386/fpu/s_remquo.S
+++ b/sysdeps/i386/fpu/s_remquo.S
@@ -31,9 +31,10 @@ ENTRY (BP_SYM (__remquo))
 	andl	$4, %ecx
 	andl	$3, %eax
 	orl	%eax, %ecx
-	movl	$0xef2960, %eax
+	leal	(%ecx,%ecx,2),%ecx
+	movl	$0xef2a60, %eax
 	shrl	%cl, %eax
-	andl	$3, %eax
+	andl	$7, %eax
 	movl	QUOP(%esp), %ecx
 	CHECK_BOUNDS_BOTH_WIDE (%ecx, QUOP(%esp), $4)
 	movl	DVDND+4(%esp), %edx
diff --git a/sysdeps/i386/fpu/s_remquof.S b/sysdeps/i386/fpu/s_remquof.S
index 90f46dcbc0..d3c5965be4 100644
--- a/sysdeps/i386/fpu/s_remquof.S
+++ b/sysdeps/i386/fpu/s_remquof.S
@@ -31,9 +31,10 @@ ENTRY (BP_SYM (__remquof))
 	andl	$4, %ecx
 	andl	$3, %eax
 	orl	%eax, %ecx
-	movl	$0xef2960, %eax
+	leal	(%ecx,%ecx,2),%ecx
+	movl	$0xef2a60, %eax
 	shrl	%cl, %eax
-	andl	$3, %eax
+	andl	$7, %eax
 	movl	QUOP(%esp), %ecx
 	CHECK_BOUNDS_BOTH_WIDE (%ecx, QUOP(%esp), $4)
 	movl	DVDND(%esp), %edx
diff --git a/sysdeps/i386/fpu/s_remquol.S b/sysdeps/i386/fpu/s_remquol.S
index c71c002243..65240adbe4 100644
--- a/sysdeps/i386/fpu/s_remquol.S
+++ b/sysdeps/i386/fpu/s_remquol.S
@@ -31,9 +31,10 @@ ENTRY (BP_SYM (__remquol))
 	andl	$4, %ecx
 	andl	$3, %eax
 	orl	%eax, %ecx
-	movl	$0xef2960, %eax
+	leal	(%ecx,%ecx,2),%ecx
+	movl	$0xef2a60, %eax
 	shrl	%cl, %eax
-	andl	$3, %eax
+	andl	$7, %eax
 	movl	QUOP(%esp), %ecx
 	CHECK_BOUNDS_BOTH_WIDE (%ecx, QUOP(%esp), $4)
 	movl	DVDND+8(%esp), %edx
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)