136 	unsigned int ix, iy;  in denormal_mulf()  local
140 	ix = hx & 0x7fffffff;  in denormal_mulf()
142 	if (iy < 0x00800000 || ix == 0)  in denormal_mulf()
146 	ix &= 0x007fffff;  in denormal_mulf()
148 	m = (unsigned long long)ix * iy;  in denormal_mulf()
156 		ix = ((int) (m >> (w - 23)) & 0x007fffff) | (exp << 23);  in denormal_mulf()
158 		ix = (int) (m >> (w - 22 - exp)) & 0x007fffff;  in denormal_mulf()
160 		ix = 0;  in denormal_mulf()
162 	ix |= (hx ^ hy) & 0x80000000;  in denormal_mulf()
163 	return ix;  in denormal_mulf()
203 	unsigned long long ix, iy;  in denormal_muld()  local
207 	ix = hx & 0x7fffffffffffffffLL;  in denormal_muld()
209 	if (iy < 0x0010000000000000LL || ix == 0)  in denormal_muld()
213 	ix &= 0x000fffffffffffffLL;  in denormal_muld()
215 	mult64(ix, iy, &mh, &ml);  in denormal_muld()
228 		ix = (rshift64(mh, ml, w - 52) & 0x000fffffffffffffLL)  in denormal_muld()
231 		ix = rshift64(mh, ml, w - 51 - exp) & 0x000fffffffffffffLL;  in denormal_muld()
233 		ix = 0;  in denormal_muld()
235 	ix |= (hx ^ hy) & 0x8000000000000000LL;  in denormal_muld()
236 	return ix;  in denormal_muld()
240 static int denormal_subf1(unsigned int ix, unsigned int iy)  in denormal_subf1()  argument
245 	if (ix < 0x00800000)  in denormal_subf1()
246 		return ix - iy;  in denormal_subf1()
248 	exp = (ix & 0x7f800000) >> 23;  in denormal_subf1()
250 		return ix;  in denormal_subf1()
253 		return ix;  in denormal_subf1()
255 	frac = (ix & 0x007fffff) | 0x00800000;  in denormal_subf1()
267 static int denormal_addf1(unsigned int ix, unsigned int iy)  in denormal_addf1()  argument
272 	if (ix < 0x00800000)  in denormal_addf1()
273 		return ix + iy;  in denormal_addf1()
275 	exp = (ix & 0x7f800000) >> 23;  in denormal_addf1()
277 		return ix;  in denormal_addf1()
280 	  return ix;  in denormal_addf1()
282 	frac = (ix & 0x007fffff) | 0x00800000;  in denormal_addf1()
294 	unsigned int ix, iy;  in denormal_addf()  local
299 		ix = hx & 0x7fffffff;  in denormal_addf()
302 			ix = denormal_subf1(ix, iy);  in denormal_addf()
303 			if ((int) ix < 0) {  in denormal_addf()
304 				ix = -ix;  in denormal_addf()
308 			ix = denormal_subf1(iy, ix);  in denormal_addf()
313 		ix = hx & 0x7fffffff;  in denormal_addf()
316 			ix = denormal_addf1(ix, iy);  in denormal_addf()
318 			ix = denormal_addf1(iy, ix);  in denormal_addf()
321 	return sign | ix;  in denormal_addf()
325 static long long denormal_subd1(unsigned long long ix, unsigned long long iy)  in denormal_subd1()  argument
330 	if (ix < 0x0010000000000000LL)  in denormal_subd1()
331 		return ix - iy;  in denormal_subd1()
333 	exp = (ix & 0x7ff0000000000000LL) >> 52;  in denormal_subd1()
335 		return ix;  in denormal_subd1()
338 		return ix;  in denormal_subd1()
340 	frac = (ix & 0x000fffffffffffffLL) | 0x0010000000000000LL;  in denormal_subd1()
352 static long long denormal_addd1(unsigned long long ix, unsigned long long iy)  in denormal_addd1()  argument
357 	if (ix < 0x0010000000000000LL)  in denormal_addd1()
358 		return ix + iy;  in denormal_addd1()
360 	exp = (ix & 0x7ff0000000000000LL) >> 52;  in denormal_addd1()
362 		return ix;  in denormal_addd1()
365 	  return ix;  in denormal_addd1()
367 	frac = (ix & 0x000fffffffffffffLL) | 0x0010000000000000LL;  in denormal_addd1()
379 	unsigned long long ix, iy;  in denormal_addd()  local
384 		ix = hx & 0x7fffffffffffffffLL;  in denormal_addd()
387 			ix = denormal_subd1(ix, iy);  in denormal_addd()
388 			if ((int) ix < 0) {  in denormal_addd()
389 				ix = -ix;  in denormal_addd()
393 			ix = denormal_subd1(iy, ix);  in denormal_addd()
398 		ix = hx & 0x7fffffffffffffffLL;  in denormal_addd()
401 			ix = denormal_addd1(ix, iy);  in denormal_addd()
403 			ix = denormal_addd1(iy, ix);  in denormal_addd()
406 	return sign | ix;  in denormal_addd()