diff options
| -rw-r--r-- | ChangeLog | 22 | ||||
| -rw-r--r-- | sysdeps/i386/fpu/s_remquo.S | 5 | ||||
| -rw-r--r-- | sysdeps/i386/fpu/s_remquof.S | 5 | ||||
| -rw-r--r-- | sysdeps/i386/fpu/s_remquol.S | 5 | ||||
| -rw-r--r-- | sysdeps/ieee754/ldbl-96/s_scalbnl.c | 16 |
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) |