22     int32_t n,hx,hy,hz,ix,iy,sx,i;  in bsd__ieee754_fmodf()  local
24     GET_FLOAT_WORD(hx,x);  in bsd__ieee754_fmodf()
26     sx = hx&0x80000000;        /* sign of x */  in bsd__ieee754_fmodf()
27     hx ^=sx;                   /* |x| */  in bsd__ieee754_fmodf()
31     if(hy==0||(hx>=0x7f800000)||           /* y=0,or x not finite */  in bsd__ieee754_fmodf()
34     if(hx<hy) return x;                    /* |x|<|y| return x */  in bsd__ieee754_fmodf()
35     if(hx==hy)  in bsd__ieee754_fmodf()
39     if(hx<0x00800000) {    /* subnormal x */  in bsd__ieee754_fmodf()
40         for (ix = -126,i=(hx<<8); i>0; i<<=1) ix -=1;  in bsd__ieee754_fmodf()
41     } else ix = (hx>>23)-127;  in bsd__ieee754_fmodf()
50         hx = 0x00800000|(0x007fffff&hx);  in bsd__ieee754_fmodf()
53         hx = hx<<n;  in bsd__ieee754_fmodf()
65         hz=hx-hy;  in bsd__ieee754_fmodf()
66         if(hz<0){hx = hx+hx;}  in bsd__ieee754_fmodf()
70             hx = hz+hz;  in bsd__ieee754_fmodf()
73     hz=hx-hy;  in bsd__ieee754_fmodf()
74     if(hz>=0) {hx=hz;}  in bsd__ieee754_fmodf()
77     if(hx==0)                     /* return sign(x)*0 */  in bsd__ieee754_fmodf()
79     while(hx<0x00800000) {        /* normalize x */  in bsd__ieee754_fmodf()
80         hx = hx+hx;  in bsd__ieee754_fmodf()
84         hx = ((hx-0x00800000)|((iy+127)<<23));  in bsd__ieee754_fmodf()
85         SET_FLOAT_WORD(x,hx|sx);  in bsd__ieee754_fmodf()
88         hx >>= n;  in bsd__ieee754_fmodf()
89         SET_FLOAT_WORD(x,hx|sx);  in bsd__ieee754_fmodf()