math: fix nextafter and nexttoward on maxdbl and maxflt input

old code (return x+x;) returns correct value and raises correct
flags only if the result is stored as double (or float)

diff --git a/src/math/nextafter.c b/src/math/nextafter.c
index a3b42c99d197e4fc925653e4535c1ba95c74cb8b..9ee82518b8f1d97a542856d1b524c8bd7b8d75ac 100644 (file)
--- a/src/math/nextafter.c
+++ b/src/math/nextafter.c
@@ -27,7 +27,7 @@ double nextafter(double x, double y)
        e = ux.bits >> 52 & 0x7ff;
        /* raise overflow if ux.value is infinite and x is finite */
        if (e == 0x7ff)
-               return x + x;
+               FORCE_EVAL(x+x);
        /* raise underflow if ux.value is subnormal or zero */
        if (e == 0)
                FORCE_EVAL(x*x + ux.value*ux.value);

diff --git a/src/math/nextafterf.c b/src/math/nextafterf.c
index b703487bbca25211cffb54dd000ebe5c34d69189..22b61dce7be9fb914ded838b8ace2ed2c43ace05 100644 (file)
--- a/src/math/nextafterf.c
+++ b/src/math/nextafterf.c
@@ -26,7 +26,7 @@ float nextafterf(float x, float y)
        e = ux.bits & 0x7f800000;
        /* raise overflow if ux.value is infinite and x is finite */
        if (e == 0x7f800000)
-               return x + x;
+               FORCE_EVAL(x+x);
        /* raise underflow if ux.value is subnormal or zero */
        if (e == 0)
                FORCE_EVAL(x*x + ux.value*ux.value);

diff --git a/src/math/nexttoward.c b/src/math/nexttoward.c
index 7355f2f19d5afa71fb7a985d72d724ca74d59fd8..6f32eca47b6814b03149006ef6f6e063eeec1507 100644 (file)
--- a/src/math/nexttoward.c
+++ b/src/math/nexttoward.c
@@ -36,7 +36,7 @@ double nexttoward(double x, long double y)
        e = ux.bits>>52 & 0x7ff;
        /* raise overflow if ux.value is infinite and x is finite */
        if (e == 0x7ff)
-               return x + x;
+               FORCE_EVAL(x+x);
        /* raise underflow if ux.value is subnormal or zero */
        if (e == 0)
                FORCE_EVAL(x*x + ux.value*ux.value);

diff --git a/src/math/nexttowardf.c b/src/math/nexttowardf.c
index 8648be6a684a669a2fd3a37784ccdac3733f060a..9a693b1a2da566093014b09e5ed9901dc0405c54 100644 (file)
--- a/src/math/nexttowardf.c
+++ b/src/math/nexttowardf.c
@@ -28,7 +28,7 @@ float nexttowardf(float x, long double y)
        e = ux.bits & 0x7f800000;
        /* raise overflow if ux.value is infinite and x is finite */
        if (e == 0x7f800000)
-               return x + x;
+               FORCE_EVAL(x+x);
        /* raise underflow if ux.value is subnormal or zero */
        if (e == 0)
                FORCE_EVAL(x*x + ux.value*ux.value);