Handle numbers > 1.0 correctly.

diff --git a/sysdeps/libm-ieee754/s_modfl.c b/sysdeps/libm-ieee754/s_modfl.c
index f97b368d6d4a34c4b3f9039d7991c6a94d5dc6f8..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)(0x7fffffff))>>(j0-20);
+           i = ((u_int32_t)(0x7fffffff))>>(j0-32);
            if((i1&i)==0) {             /* x is integral */
                *iptr = x;
                INSERT_WORDS(x,se&0x8000,0);    /* return +-0 */