From f8da88b080504633c5442a4fcd234cb63912a7bf Mon Sep 17 00:00:00 2001 From: Ulrich Drepper Date: Mon, 4 Dec 2000 04:56:36 +0000 Subject: Update. * sysdeps/i386/fpu/s_nextafterl.c (__nextafterl): Hopefully fix the handling of denormalized numbers. * math/test-misc.c (main): Test for fpclassify and pseudo denormal numbers on x86. * sysdeps/i386/fpu/fpclassifyl.c: New file. --- sysdeps/i386/fpu/s_nextafterl.c | 26 ++++++++++++++++++++++---- 1 file changed, 22 insertions(+), 4 deletions(-) (limited to 'sysdeps') diff --git a/sysdeps/i386/fpu/s_nextafterl.c b/sysdeps/i386/fpu/s_nextafterl.c index bca7181..a905858 100644 --- a/sysdeps/i386/fpu/s_nextafterl.c +++ b/sysdeps/i386/fpu/s_nextafterl.c @@ -59,8 +59,17 @@ static char rcsid[] = "$NetBSD: $"; if(ix>iy||((ix==iy) && (hx>hy||((hx==hy)&&(lx>ly))))) { /* x > y, x -= ulp */ if(lx==0) { - if (hx==0) esx -= 1; - hx -= 1; + if (hx <= 0x80000000) { + if (esx == 0) + hx = 0; + else { + esx -= 1; + hx = hx - 1; + if (esx > 0) + hx |= 0x80000000; + } + } else + hx -= 1; } lx -= 1; } else { /* x < y, x += ulp */ @@ -75,8 +84,17 @@ static char rcsid[] = "$NetBSD: $"; if(esy>=0||(ix>iy||((ix==iy)&&(hx>hy||((hx==hy)&&(lx>ly)))))){ /* x < y, x -= ulp */ if(lx==0) { - if (hx==0) esx -= 1; - hx -= 1; + if (hx <= 0x80000000) { + if (esx == 0) + hx = 0; + else { + esx -= 1; + hx = hx - 1; + if (esx > 0) + hx |= 0x80000000; + } + } else + hx -= 1; } lx -= 1; } else { /* x > y, x += ulp */ -- cgit v1.1