math: Use 'lsl' for left shifts in fmod/fmod

These two functions are counting loops by waiting for
a single bit stored as a signed value to go negative by
shifting into the top bit. Use our lsl macro, which converts
this to defined behavior.

Signed-off-by: Keith Packard <[email protected]>

Author: keith-packard

newlib/libm/math/s_fmod.c

@@ -57,10 +57,10 @@ fmod64(__float64 x, __float64 y)
     /* determine ix = ilogb(x) */
     if (hx < 0x00100000) { /* subnormal x */
         if (hx == 0) {
-            for (ix = -1043, i = lx; i > 0; i <<= 1)
+            for (ix = -1043, i = lx; i > 0; i = lsl(i, 1))
                 ix -= 1;
         } else {
-            for (ix = -1022, i = (hx << 11); i > 0; i <<= 1)
+            for (ix = -1022, i = lsl(hx, 11); i > 0; i = lsl(i, 1))
                 ix -= 1;
         }
     } else
@@ -69,10 +69,10 @@ fmod64(__float64 x, __float64 y)
     /* determine iy = ilogb(y) */
     if (hy < 0x00100000) { /* subnormal y */
         if (hy == 0) {
-            for (iy = -1043, i = ly; i > 0; i <<= 1)
+            for (iy = -1043, i = ly; i > 0; i = lsl(i, 1))
                 iy -= 1;
         } else {
-            for (iy = -1022, i = (hy << 11); i > 0; i <<= 1)
+            for (iy = -1022, i = lsl(hy, 11); i > 0; i = lsl(i, 1))
                 iy -= 1;
         }
     } else

newlib/libm/math/sf_fmod.c

@@ -56,14 +56,14 @@ fmodf(float x, float y)
 
     /* determine ix = ilogb(x) */
     if (FLT_UWORD_IS_SUBNORMAL(hx)) { /* subnormal x */
-        for (ix = -126, i = (hx << 8); i > 0; i <<= 1)
+        for (ix = -126, i = lsl(hx, 8); i > 0; i = lsl(i, 1))
             ix -= 1;
     } else
         ix = (hx >> 23) - 127;
 
     /* determine iy = ilogb(y) */
     if (FLT_UWORD_IS_SUBNORMAL(hy)) { /* subnormal y */
-        for (iy = -126, i = (hy << 8); i >= 0; i <<= 1)
+        for (iy = -126, i = lsl(hy, 8); i >= 0; i = lsl(i, 1))
             iy -= 1;
     } else
         iy = (hy >> 23) - 127;