41   int n, hx, hy, hz, ix, iy, sx, i;  in fmod()  local
44   hx = __HI (x); /* high word of x */  in fmod()
48   sx = hx & 0x80000000; /* sign of x */  in fmod()
49   hx ^= sx; /* |x| */  in fmod()
53   if ((hy | ly) == 0 || (hx >= 0x7ff00000) || /* y = 0, or x not finite */  in fmod()
58   if (hx <= hy)  in fmod()
60     if ((hx < hy) || (lx < ly)) /* |x| < |y| return x */  in fmod()
71   if (hx < 0x00100000) /* subnormal x */  in fmod()
73     if (hx == 0)  in fmod()
82       for (ix = -1022, i = (hx << 11); i > 0; i <<= 1)  in fmod()
90     ix = (hx >> 20) - 1023;  in fmod()
119     hx = 0x00100000 | (0x000fffff & hx);  in fmod()
126       hx = (((unsigned int) hx) << n) | (lx >> (32 - n));  in fmod()
131       hx = lx << (n - 32);  in fmod()
158     hz = hx - hy;  in fmod()
166       hx = hx + hx + (lx >> 31);  in fmod()
175       hx = hz + hz + (lz >> 31);  in fmod()
179   hz = hx - hy;  in fmod()
187     hx = hz;  in fmod()
192   if ((hx | lx) == 0) /* return sign(x) * 0 */  in fmod()
196   while (hx < 0x00100000) /* normalize x */  in fmod()
198     hx = hx + hx + (lx >> 31);  in fmod()
206     hx = ((hx - 0x00100000) | ((iy + 1023) << 20));  in fmod()
207     ret.as_int.hi = hx | sx;  in fmod()
215       lx = (lx >> n) | ((unsigned) hx << (32 - n));  in fmod()
216       hx >>= n;  in fmod()
220       lx = (hx << (32 - n)) | (lx >> n);  in fmod()
221       hx = sx;  in fmod()
225       lx = hx >> (n - 32);  in fmod()
226       hx = sx;  in fmod()
228     ret.as_int.hi = hx | sx;  in fmod()