24 static const double one = 1.0, Zero[] = {0.0, -0.0,};
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()
55 	} else ix = (hx>>20)-1023;  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()
64 	} else iy = (hy>>20)-1023;  in fmod()
67 	if(ix >= -1022)   in fmod()
68 	    hx = 0x00100000|(0x000fffff&hx);  in fmod()
70 	    n = -1022-ix;  in fmod()
72 	        hx = (hx<<n)|(lx>>(32-n));  in fmod()
75 		hx = lx<<(n-32);  in fmod()
76 		lx = 0;  in fmod()
79 	if(iy >= -1022)   in fmod()
80 	    hy = 0x00100000|(0x000fffff&hy);  in fmod()
82 	    n = -1022-iy;  in fmod()
84 	        hy = (hy<<n)|(ly>>(32-n));  in fmod()
87 		hy = ly<<(n-32);  in fmod()
88 		ly = 0;  in fmod()
93 	n = ix - iy;  in fmod()
94 	while(n--) {  in fmod()
95 	    hz=hx-hy;lz=lx-ly; if(lx<ly) hz -= 1;  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()
103 	hz=hx-hy;lz=lx-ly; if(lx<ly) hz -= 1;  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()
111 	    iy -= 1;  in fmod()
113 	if(iy>= -1022) {	/* normalize output */  in fmod()
114 	    hx = ((hx-0x00100000)|((iy+1023)<<20));  in fmod()
117 	    n = -1022 - iy;  in fmod()
119 		lx = (lx>>n)|((u_int32_t)hx<<(32-n));  in fmod()
122 		lx = (hx<<(32-n))|(lx>>n); hx = sx;  in fmod()
124 		lx = hx>>(n-32); hx = sx;  in fmod()