aboutsummaryrefslogtreecommitdiff
path: root/sysdeps/ieee754
diff options
context:
space:
mode:
authorUlrich Drepper <drepper@redhat.com>2003-12-07 21:22:46 +0000
committerUlrich Drepper <drepper@redhat.com>2003-12-07 21:22:46 +0000
commit07449987c90864fceb487994b4aaa5170803ab24 (patch)
tree76771090b9cb1866d55b5986c677e51dd768c3af /sysdeps/ieee754
parent64b02fd2f5fd5d2e2b7187b5e982fa6b59587500 (diff)
downloadglibc-07449987c90864fceb487994b4aaa5170803ab24.tar
glibc-07449987c90864fceb487994b4aaa5170803ab24.tar.gz
glibc-07449987c90864fceb487994b4aaa5170803ab24.tar.bz2
glibc-07449987c90864fceb487994b4aaa5170803ab24.zip
Update.
2003-12-07 Ulrich Drepper <drepper@redhat.com> * sysdeps/i386/fpu/s_nexttowardf.c: Construct overflow value correctly. * sysdeps/i386/fpu/s_nexttoward.c: Likewise. * sysdeps/ieee754/ldbl-128/s_nexttoward.c: Likewise. * sysdeps/ieee754/ldbl-96/s_nexttoward.c: Likewise. * math/Makefile (tests): Add bug-nexttoward. * math/bug-nexttowward.c: New file. * sysdeps/generic/s_nextafter.c: Make sure overflow exception is set. * sysdeps/ieee754/flt-32/s_nextafterf.c: Likewise. * math/bug-nextafter.c (main): Add tests for overflow and negative values.
Diffstat (limited to 'sysdeps/ieee754')
-rw-r--r--sysdeps/ieee754/flt-32/s_nextafterf.c10
-rw-r--r--sysdeps/ieee754/ldbl-128/s_nexttoward.c9
-rw-r--r--sysdeps/ieee754/ldbl-96/s_nexttoward.c9
3 files changed, 20 insertions, 8 deletions
diff --git a/sysdeps/ieee754/flt-32/s_nextafterf.c b/sysdeps/ieee754/flt-32/s_nextafterf.c
index e144e56acd..e1568e24c9 100644
--- a/sysdeps/ieee754/flt-32/s_nextafterf.c
+++ b/sysdeps/ieee754/flt-32/s_nextafterf.c
@@ -59,12 +59,10 @@ static char rcsid[] = "$NetBSD: s_nextafterf.c,v 1.4 1995/05/10 20:48:01 jtc Exp
}
hy = hx&0x7f800000;
if(hy>=0x7f800000) {
-#if FLT_EVAL_METHOD == 0
- return x+x; /* overflow */
-#else
- /* Note that y is either +Inf or -Inf. */
- return x+y;
-#endif
+ x = x+x; /* overflow */
+ if (FLT_EVAL_METHOD != 0)
+ asm ("" : "=m"(x) : "m"(x));
+ return x; /* overflow */
}
if(hy<0x00800000) { /* underflow */
y = x*x;
diff --git a/sysdeps/ieee754/ldbl-128/s_nexttoward.c b/sysdeps/ieee754/ldbl-128/s_nexttoward.c
index 2c654fb973..553e401973 100644
--- a/sysdeps/ieee754/ldbl-128/s_nexttoward.c
+++ b/sysdeps/ieee754/ldbl-128/s_nexttoward.c
@@ -27,6 +27,7 @@ static char rcsid[] = "$NetBSD: $";
#include "math.h"
#include "math_private.h"
+#include <float.h>
#ifdef __STDC__
double __nexttoward(double x, long double y)
@@ -83,7 +84,13 @@ static char rcsid[] = "$NetBSD: $";
}
}
hy = hx&0x7ff00000;
- if(hy>=0x7ff00000) return x+x; /* overflow */
+ if(hy>=0x7ff00000) {
+ x = x+x; /* overflow */
+ if (FLT_EVAL_METHOD != 0 && FLT_EVAL_METHOD != 1)
+ /* Force conversion to float. */
+ asm ("" : "=m"(x) : "m"(x));
+ return x;
+ }
if(hy<0x00100000) { /* underflow */
double x2 = x*x;
if(x2!=x) { /* raise underflow flag */
diff --git a/sysdeps/ieee754/ldbl-96/s_nexttoward.c b/sysdeps/ieee754/ldbl-96/s_nexttoward.c
index c29bcd146d..8b273af929 100644
--- a/sysdeps/ieee754/ldbl-96/s_nexttoward.c
+++ b/sysdeps/ieee754/ldbl-96/s_nexttoward.c
@@ -27,6 +27,7 @@ static char rcsid[] = "$NetBSD: $";
#include "math.h"
#include "math_private.h"
+#include <float.h>
#ifdef __STDC__
double __nexttoward(double x, long double y)
@@ -80,7 +81,13 @@ static char rcsid[] = "$NetBSD: $";
}
}
hy = hx&0x7ff00000;
- if(hy>=0x7ff00000) return x+x; /* overflow */
+ if(hy>=0x7ff00000) {
+ x = x+x; /* overflow */
+ if (FLT_EVAL_METHOD != 0 && FLT_EVAL_METHOD != 1)
+ /* Force conversion to float. */
+ asm ("" : "=m"(x) : "m"(x));
+ return x;
+ }
if(hy<0x00100000) { /* underflow */
double x2 = x*x;
if(x2!=x) { /* raise underflow flag */