96 	unsigned int ix, iy;  in denormal_mulf()  local
100 	ix = hx & 0x7fffffff;  in denormal_mulf()
102 	if (iy < 0x00800000 || ix == 0)  in denormal_mulf()
106 	ix &= 0x007fffff;  in denormal_mulf()
108 	m = (unsigned long long)ix * iy;  in denormal_mulf()
116 		ix = ((int) (m >> (w - 23)) & 0x007fffff) | (exp << 23);  in denormal_mulf()
118 		ix = (int) (m >> (w - 22 - exp)) & 0x007fffff;  in denormal_mulf()
120 		ix = 0;  in denormal_mulf()
122 	ix |= (hx ^ hy) & 0x80000000;  in denormal_mulf()
123 	return ix;  in denormal_mulf()
163 	unsigned long long ix, iy;  in denormal_muld()  local
167 	ix = hx & 0x7fffffffffffffffLL;  in denormal_muld()
169 	if (iy < 0x0010000000000000LL || ix == 0)  in denormal_muld()
173 	ix &= 0x000fffffffffffffLL;  in denormal_muld()
175 	mult64(ix, iy, &mh, &ml);  in denormal_muld()
188 		ix = (rshift64(mh, ml, w - 52) & 0x000fffffffffffffLL)  in denormal_muld()
191 		ix = rshift64(mh, ml, w - 51 - exp) & 0x000fffffffffffffLL;  in denormal_muld()
193 		ix = 0;  in denormal_muld()
195 	ix |= (hx ^ hy) & 0x8000000000000000LL;  in denormal_muld()
196 	return ix;  in denormal_muld()
200 static int denormal_subf1(unsigned int ix, unsigned int iy)  in denormal_subf1()  argument
205 	if (ix < 0x00800000)  in denormal_subf1()
206 		return ix - iy;  in denormal_subf1()
208 	exp = (ix & 0x7f800000) >> 23;  in denormal_subf1()
210 		return ix;  in denormal_subf1()
213 		return ix;  in denormal_subf1()
215 	frac = (ix & 0x007fffff) | 0x00800000;  in denormal_subf1()
227 static int denormal_addf1(unsigned int ix, unsigned int iy)  in denormal_addf1()  argument
232 	if (ix < 0x00800000)  in denormal_addf1()
233 		return ix + iy;  in denormal_addf1()
235 	exp = (ix & 0x7f800000) >> 23;  in denormal_addf1()
237 		return ix;  in denormal_addf1()
240 	  return ix;  in denormal_addf1()
242 	frac = (ix & 0x007fffff) | 0x00800000;  in denormal_addf1()
254 	unsigned int ix, iy;  in denormal_addf()  local
259 		ix = hx & 0x7fffffff;  in denormal_addf()
262 			ix = denormal_subf1(ix, iy);  in denormal_addf()
263 			if ((int) ix < 0) {  in denormal_addf()
264 				ix = -ix;  in denormal_addf()
268 			ix = denormal_subf1(iy, ix);  in denormal_addf()
273 		ix = hx & 0x7fffffff;  in denormal_addf()
276 			ix = denormal_addf1(ix, iy);  in denormal_addf()
278 			ix = denormal_addf1(iy, ix);  in denormal_addf()
281 	return sign | ix;  in denormal_addf()
285 static long long denormal_subd1(unsigned long long ix, unsigned long long iy)  in denormal_subd1()  argument
290 	if (ix < 0x0010000000000000LL)  in denormal_subd1()
291 		return ix - iy;  in denormal_subd1()
293 	exp = (ix & 0x7ff0000000000000LL) >> 52;  in denormal_subd1()
295 		return ix;  in denormal_subd1()
298 		return ix;  in denormal_subd1()
300 	frac = (ix & 0x000fffffffffffffLL) | 0x0010000000000000LL;  in denormal_subd1()
312 static long long denormal_addd1(unsigned long long ix, unsigned long long iy)  in denormal_addd1()  argument
317 	if (ix < 0x0010000000000000LL)  in denormal_addd1()
318 		return ix + iy;  in denormal_addd1()
320 	exp = (ix & 0x7ff0000000000000LL) >> 52;  in denormal_addd1()
322 		return ix;  in denormal_addd1()
325 	  return ix;  in denormal_addd1()
327 	frac = (ix & 0x000fffffffffffffLL) | 0x0010000000000000LL;  in denormal_addd1()
339 	unsigned long long ix, iy;  in denormal_addd()  local
344 		ix = hx & 0x7fffffffffffffffLL;  in denormal_addd()
347 			ix = denormal_subd1(ix, iy);  in denormal_addd()
348 			if ((int) ix < 0) {  in denormal_addd()
349 				ix = -ix;  in denormal_addd()
353 			ix = denormal_subd1(iy, ix);  in denormal_addd()
358 		ix = hx & 0x7fffffffffffffffLL;  in denormal_addd()
361 			ix = denormal_addd1(ix, iy);  in denormal_addd()
363 			ix = denormal_addd1(iy, ix);  in denormal_addd()
366 	return sign | ix;  in denormal_addd()