35 	is = -0x7fffffff - 1,
36 	im = 0x0000ffff,
37 	iu = 0x00010000;
47 #define	__H3(x) *(0 + (int *) &x)
49 #define	__H0(x) *(0 + (int *) &x)
60 	int x0, y0, z0, carry;  in fmodl()  local
72 	sx = hx & 0x80000000;  in fmodl()
73 	x0 = hx ^ sx;  in fmodl()
74 	y0 &= 0x7fffffff;  in fmodl()
77 	if (x0 >= 0x7fff0000 ||	/* !finitel(x) */  in fmodl()
78 	    (y0 > 0x7fff0000) || (y0 == 0x7fff0000 && ((y1 | y2 | y3) != 0)) ||  in fmodl()
79 	    (y0 | y1 | y2 | y3) == 0)	/* isnanl(y) || y = 0 */  in fmodl()
90 	if (x0 < iu) {		/* subnormal x */  in fmodl()
92 		while (x0 == 0) {  in fmodl()
94 			x0 = x1 >> 16;  in fmodl()
99 		while (x0 < iu) {  in fmodl()
101 			x0 = (x0 << 1) | (x1 >> 31);  in fmodl()
107 		ix = (x0 >> 16) - 16383;  in fmodl()
108 		x0 = iu | (x0 & im);  in fmodl()
114 		while (y0 == 0) {  in fmodl()
136 		while (x0 == 0 && n >= 16) {  in fmodl()
138 			x0 = x1 >> 16;  in fmodl()
143 		while (x0 < iu && n >= 1) {  in fmodl()
145 			x0 = (x0 << 1) | (x1 >> 31);  in fmodl()
150 		carry = 0;  in fmodl()
153 		if (carry == 0) {  in fmodl()
160 		if (carry == 0) {  in fmodl()
167 		z0 = x0 - y0 - carry;  in fmodl()
168 		if (z0 < 0) {	/* double x */  in fmodl()
169 			x0 = x0 + x0 + ((x1 & is) != 0);  in fmodl()
170 			x1 = x1 + x1 + ((x2 & is) != 0);  in fmodl()
171 			x2 = x2 + x2 + ((x3 & is) != 0);  in fmodl()
174 			if (z0 == 0) {  in fmodl()
175 				if ((z1 | z2 | z3) == 0) {	/* 0: done */  in fmodl()
177 					__H1(a) = __H2(a) = __H3(a) = 0;  in fmodl()
182 			z0 = z0 + z0 + ((z1 & is) != 0);  in fmodl()
183 			z1 = z1 + z1 + ((z2 & is) != 0);  in fmodl()
184 			z2 = z2 + z2 + ((z3 & is) != 0);  in fmodl()
186 			x0 = z0;  in fmodl()
193 	carry = 0;  in fmodl()
196 	if (carry == 0) {  in fmodl()
203 	if (carry == 0) {  in fmodl()
210 	z0 = x0 - y0 - carry;  in fmodl()
211 	if (z0 >= 0) {  in fmodl()
212 		x0 = z0;  in fmodl()
218 	if ((x0 | x1 | x2 | x3) == 0) {  in fmodl()
220 		__H1(a) = __H2(a) = __H3(a) = 0;  in fmodl()
223 	while (x0 < iu) {  in fmodl()
224 		if (x0 == 0) {  in fmodl()
226 			x0 = x1 >> 16;  in fmodl()
231 			x0 = x0 + x0 + ((x1 & is) != 0);  in fmodl()
232 			x1 = x1 + x1 + ((x2 & is) != 0);  in fmodl()
233 			x2 = x2 + x2 + ((x3 & is) != 0);  in fmodl()
241 		__H0(a) = sx | (x0 - iu) | ((iy + 16383) << 16);  in fmodl()
248 		if (k != 0) {  in fmodl()
252 				x1 = (x0 << (32 - k)) | (x1 >> k);  in fmodl()
253 				x0 >>= k;  in fmodl()
257 				x1 = (x0 << (32 - k)) | (x1 >> k);  in fmodl()
258 				x0 = 0;  in fmodl()
265 			x1 = x0;  in fmodl()
266 			x0 = 0;  in fmodl()
268 		__H0(a) = x0 | sx;  in fmodl()