99 	unsigned int ix, iy;  in denormal_mulf()  local
103 	ix = hx & 0x7fffffff;  in denormal_mulf()
105 	if (iy < 0x00800000 || ix == 0)  in denormal_mulf()
109 	ix &= 0x007fffff;  in denormal_mulf()
111 	m = (unsigned long long)ix * iy;  in denormal_mulf()
119 		ix = ((int) (m >> (w - 23)) & 0x007fffff) | (exp << 23);  in denormal_mulf()
121 		ix = (int) (m >> (w - 22 - exp)) & 0x007fffff;  in denormal_mulf()
123 		ix = 0;  in denormal_mulf()
125 	ix |= (hx ^ hy) & 0x80000000;  in denormal_mulf()
126 	return ix;  in denormal_mulf()
166 	unsigned long long ix, iy;  in denormal_muld()  local
170 	ix = hx & 0x7fffffffffffffffLL;  in denormal_muld()
172 	if (iy < 0x0010000000000000LL || ix == 0)  in denormal_muld()
176 	ix &= 0x000fffffffffffffLL;  in denormal_muld()
178 	mult64(ix, iy, &mh, &ml);  in denormal_muld()
191 		ix = (rshift64(mh, ml, w - 52) & 0x000fffffffffffffLL)  in denormal_muld()
194 		ix = rshift64(mh, ml, w - 51 - exp) & 0x000fffffffffffffLL;  in denormal_muld()
196 		ix = 0;  in denormal_muld()
198 	ix |= (hx ^ hy) & 0x8000000000000000LL;  in denormal_muld()
199 	return ix;  in denormal_muld()
203 static int denormal_subf1(unsigned int ix, unsigned int iy)  in denormal_subf1()  argument
208 	if (ix < 0x00800000)  in denormal_subf1()
209 		return ix - iy;  in denormal_subf1()
211 	exp = (ix & 0x7f800000) >> 23;  in denormal_subf1()
213 		return ix;  in denormal_subf1()
216 		return ix;  in denormal_subf1()
218 	frac = (ix & 0x007fffff) | 0x00800000;  in denormal_subf1()
230 static int denormal_addf1(unsigned int ix, unsigned int iy)  in denormal_addf1()  argument
235 	if (ix < 0x00800000)  in denormal_addf1()
236 		return ix + iy;  in denormal_addf1()
238 	exp = (ix & 0x7f800000) >> 23;  in denormal_addf1()
240 		return ix;  in denormal_addf1()
243 	  return ix;  in denormal_addf1()
245 	frac = (ix & 0x007fffff) | 0x00800000;  in denormal_addf1()
257 	unsigned int ix, iy;  in denormal_addf()  local
262 		ix = hx & 0x7fffffff;  in denormal_addf()
265 			ix = denormal_subf1(ix, iy);  in denormal_addf()
266 			if ((int) ix < 0) {  in denormal_addf()
267 				ix = -ix;  in denormal_addf()
271 			ix = denormal_subf1(iy, ix);  in denormal_addf()
276 		ix = hx & 0x7fffffff;  in denormal_addf()
279 			ix = denormal_addf1(ix, iy);  in denormal_addf()
281 			ix = denormal_addf1(iy, ix);  in denormal_addf()
284 	return sign | ix;  in denormal_addf()
288 static long long denormal_subd1(unsigned long long ix, unsigned long long iy)  in denormal_subd1()  argument
293 	if (ix < 0x0010000000000000LL)  in denormal_subd1()
294 		return ix - iy;  in denormal_subd1()
296 	exp = (ix & 0x7ff0000000000000LL) >> 52;  in denormal_subd1()
298 		return ix;  in denormal_subd1()
301 		return ix;  in denormal_subd1()
303 	frac = (ix & 0x000fffffffffffffLL) | 0x0010000000000000LL;  in denormal_subd1()
315 static long long denormal_addd1(unsigned long long ix, unsigned long long iy)  in denormal_addd1()  argument
320 	if (ix < 0x0010000000000000LL)  in denormal_addd1()
321 		return ix + iy;  in denormal_addd1()
323 	exp = (ix & 0x7ff0000000000000LL) >> 52;  in denormal_addd1()
325 		return ix;  in denormal_addd1()
328 	  return ix;  in denormal_addd1()
330 	frac = (ix & 0x000fffffffffffffLL) | 0x0010000000000000LL;  in denormal_addd1()
342 	unsigned long long ix, iy;  in denormal_addd()  local
347 		ix = hx & 0x7fffffffffffffffLL;  in denormal_addd()
350 			ix = denormal_subd1(ix, iy);  in denormal_addd()
351 			if ((int) ix < 0) {  in denormal_addd()
352 				ix = -ix;  in denormal_addd()
356 			ix = denormal_subd1(iy, ix);  in denormal_addd()
361 		ix = hx & 0x7fffffffffffffffLL;  in denormal_addd()
364 			ix = denormal_addd1(ix, iy);  in denormal_addd()
366 			ix = denormal_addd1(iy, ix);  in denormal_addd()
369 	return sign | ix;  in denormal_addd()