math: fix nextafter and nexttoward on maxdbl and maxflt input

[musl] / src / math / nexttoward.c

diff --git a/src/math/nexttoward.c b/src/math/nexttoward.c
index 7355f2f..6f32eca 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);