math: fix nextafter and nexttoward on maxdbl and maxflt input

[musl] / src / math / nexttowardf.c

diff --git a/src/math/nexttowardf.c b/src/math/nexttowardf.c
index 821f72a..9a693b1 100644 (file)
--- a/src/math/nexttowardf.c
+++ b/src/math/nexttowardf.c
@@ -28,10 +28,9 @@ 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) {
-               volatile float z = x*x + ux.value*ux.value;
-       }
+       if (e == 0)
+               FORCE_EVAL(x*x + ux.value*ux.value);
        return ux.value;
 }