Correct handling of arguments with exponent >31 and < 63.

diff --git a/sysdeps/libm-ieee754/s_modfl.c b/sysdeps/libm-ieee754/s_modfl.c
index ad16ef65ec0e4ed5ba52ecb6eae6d7165f89c2be..5f2b4b5e2c44f09022ac1794ad28180930722dc3 100644 (file)
--- a/sysdeps/libm-ieee754/s_modfl.c
+++ b/sysdeps/libm-ieee754/s_modfl.c
@@ -71,7 +71,7 @@ static long double one = 1.0;
            SET_LDOUBLE_WORDS(x,se&0x8000,0,0); /* return +-0 */
            return x;
        } else {                        /* fraction part in low x */
-           i = ((u_int32_t)(0xffffffff))>>(j0-32);
+           i = ((u_int32_t)(0x7fffffff))>>(j0-32);
            if((i1&i)==0) {             /* x is integral */
                *iptr = x;
                INSERT_WORDS(x,se&0x8000,0);    /* return +-0 */