34 	sx = hx&0x80000000;		/* sign of x */  in fmod()
36 	hy &= 0x7fffffff;	/* |y| */  in fmod()
39 	if((hy|ly)==0||(hx>=0x7ff00000)||	/* y=0,or x not finite */  in fmod()
40 	  ((hy|((ly|-ly)>>31))>0x7ff00000))	/* or y is NaN */  in fmod()
45 		return Zero[(u_int32_t)sx>>31];	/* |x|=|y| return x*0*/  in fmod()
49 	if(hx<0x00100000) {	/* subnormal x */  in fmod()
50 	    if(hx==0) {  in fmod()
51 		for (ix = -1043, i=lx; i>0; i<<=1) ix -=1;  in fmod()
53 		for (ix = -1022,i=(hx<<11); i>0; i<<=1) ix -=1;  in fmod()
58 	if(hy<0x00100000) {	/* subnormal y */  in fmod()
59 	    if(hy==0) {  in fmod()
60 		for (iy = -1043, i=ly; i>0; i<<=1) iy -=1;  in fmod()
62 		for (iy = -1022,i=(hy<<11); i>0; i<<=1) iy -=1;  in fmod()
68 	    hx = 0x00100000|(0x000fffff&hx);  in fmod()
76 		lx = 0;  in fmod()
80 	    hy = 0x00100000|(0x000fffff&hy);  in fmod()
88 		ly = 0;  in fmod()
96 	    if(hz<0){hx = hx+hx+(lx>>31); lx = lx+lx;}  in fmod()
98 	    	if((hz|lz)==0) 		/* return sign(x)*0 */  in fmod()
104 	if(hz>=0) {hx=hz;lx=lz;}  in fmod()
107 	if((hx|lx)==0) 			/* return sign(x)*0 */  in fmod()
109 	while(hx<0x00100000) {		/* normalize x */  in fmod()
114 	    hx = ((hx-0x00100000)|((iy+1023)<<20));  in fmod()