32 	int32_t n,hx,hy,hz,ix,iy,sx,i;  in __ieee754_fmod()  local
35 	EXTRACT_WORDS(hx,lx,x);  in __ieee754_fmod()
37 	sx = hx&0x80000000;		/* sign of x */  in __ieee754_fmod()
38 	hx ^=sx;		/* |x| */  in __ieee754_fmod()
42 	if((hy|ly)==0||(hx>=0x7ff00000)||	/* y=0,or x not finite */  in __ieee754_fmod()
45 	if(hx<=hy) {  in __ieee754_fmod()
46 	    if((hx<hy)||(lx<ly)) return x;	/* |x|<|y| return x */  in __ieee754_fmod()
52 	if(hx<0x00100000) {	/* subnormal x */  in __ieee754_fmod()
53 	    if(hx==0) {  in __ieee754_fmod()
56 		for (ix = -1022,i=(hx<<11); i>0; i<<=1) ix -=1;  in __ieee754_fmod()
58 	} else ix = (hx>>20)-1023;  in __ieee754_fmod()
71 	    hx = 0x00100000|(0x000fffff&hx);  in __ieee754_fmod()
75 	        hx = (hx<<n)|(lx>>(32-n));  in __ieee754_fmod()
78 		hx = lx<<(n-32);  in __ieee754_fmod()
98 	    hz=hx-hy;lz=lx-ly; if(lx<ly) hz -= 1;  in __ieee754_fmod()
99 	    if(hz<0){hx = hx+hx+(lx>>31); lx = lx+lx;}  in __ieee754_fmod()
103 	    	hx = hz+hz+(lz>>31); lx = lz+lz;  in __ieee754_fmod()
106 	hz=hx-hy;lz=lx-ly; if(lx<ly) hz -= 1;  in __ieee754_fmod()
107 	if(hz>=0) {hx=hz;lx=lz;}  in __ieee754_fmod()
110 	if((hx|lx)==0) 			/* return sign(x)*0 */  in __ieee754_fmod()
112 	while(hx<0x00100000) {		/* normalize x */  in __ieee754_fmod()
113 	    hx = hx+hx+(lx>>31); lx = lx+lx;  in __ieee754_fmod()
117 	    hx = ((hx-0x00100000)|((iy+1023)<<20));  in __ieee754_fmod()
118 	    INSERT_WORDS(x,hx|sx,lx);  in __ieee754_fmod()
122 		lx = (lx>>n)|((u_int32_t)hx<<(32-n));  in __ieee754_fmod()
123 		hx >>= n;  in __ieee754_fmod()
125 		lx = (hx<<(32-n))|(lx>>n); hx = sx;  in __ieee754_fmod()
127 		lx = hx>>(n-32); hx = sx;  in __ieee754_fmod()
129 	    INSERT_WORDS(x,hx|sx,lx);  in __ieee754_fmod()