support ld80 pseudo-denormal invalid bit patterns; treat them as nan

[oweals/musl.git] / src / math / nextafter.c

diff --git a/src/math/nextafter.c b/src/math/nextafter.c
index e4bfb022cb8d81471764a393442e024188fa8813..a3b42c99d197e4fc925653e4535c1ba95c74cb8b 100644 (file)
--- a/src/math/nextafter.c
+++ b/src/math/nextafter.c
@@ -29,9 +29,7 @@ double nextafter(double x, double y)
        if (e == 0x7ff)
                return x + x;
        /* raise underflow if ux.value is subnormal or zero */
-       if (e == 0) {
-               volatile double z;
-               z = x*x + ux.value*ux.value;
-       }
+       if (e == 0)
+               FORCE_EVAL(x*x + ux.value*ux.value);
        return ux.value;
 }