math/ld/ld128: Use 'lsl' for left shifts in fmodl, hypotl and remquol

These all iterate until a value flips sign by shifting a bit left, and
that only works if you do the shift using unsigned arithmetic and then
cast back to signed.

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

Author: keith-packard

newlib/libm/ld/ld128/e_fmodl.c

@@ -50,18 +50,18 @@ fmodl(long double x, long double y)
     /* determine ix = ilogb(x) */
 	if(hx<0x0001000000000000LL) {	/* subnormal x */
 	    if(hx==0) {
-		for (ix = -16431, i=lx; i>0; i<<=1) ix -=1;
+		for (ix = -16431, i=lx; i>0; i=lsl(i, 1)) ix -=1;
 	    } else {
-		for (ix = -16382, i=hx<<15; i>0; i<<=1) ix -=1;
+		for (ix = -16382, i=lsl(hx, 15); i>0; i=lsl(i, 1)) ix -=1;
 	    }
 	} else ix = (hx>>48)-0x3fff;
 
     /* determine iy = ilogb(y) */
 	if(hy<0x0001000000000000LL) {	/* subnormal y */
 	    if(hy==0) {
-		for (iy = -16431, i=ly; i>0; i<<=1) iy -=1;
+		for (iy = -16431, i=ly; i>0; i=lsl(i, 1)) iy -=1;
 	    } else {
-		for (iy = -16382, i=hy<<15; i>0; i<<=1) iy -=1;
+		for (iy = -16382, i=lsl(hy, 15); i>0; i=lsl(i, 1)) iy -=1;
 	    }
 	} else iy = (hy>>48)-0x3fff;
 

newlib/libm/ld/ld128/e_hypotl.c

@@ -116,7 +116,7 @@ hypotl(long double x, long double y)
 	    u_int64_t high;
 	    t1 = 1.0L;
 	    GET_LDOUBLE_MSW64(high,t1);
-	    SET_LDOUBLE_MSW64(t1,high+(k<<48));
+	    SET_LDOUBLE_MSW64(t1,high+lsl(k, 48));
 	    return check_oflowl(t1*w);
 	} else return w;
 }

newlib/libm/ld/ld128/s_remquol.c

@@ -115,7 +115,7 @@ remquol(long double x, long double y, int *quo)
 	    hz=_hx-hy;lz=lx-ly; if(lx<ly) hz -= 1;
 	    if(hz<0){_hx = _hx+_hx+(lx>>MANL_SHIFT); lx = lx+lx;}
 	    else {_hx = hz+hz+(lz>>MANL_SHIFT); lx = lz+lz; q++;}
-	    q <<= 1;
+	    q = lsl(q, 1);
 	}
 	hz=_hx-hy;lz=lx-ly; if(lx<ly) hz -= 1;
 	if(hz>=0) {_hx=hz;lx=lz;q++;}