parisc/math-emu/sfadd.c

42     register unsigned int left, right, result, extent;
75                     Sgl_makequietnan(result);
76 		    *dstptr = result;
185 		/* Left is not a zero and must be the result.  Trapped
208 	    /* Both operands are denormalized.  The result must be exact
213 		Sgl_subtract(left,/*minus*/right,/*into*/result);
214 		if(Sgl_iszero_mantissa(result))
218 			Sgl_setone_sign(result);
222 			Sgl_setzero_sign(result);
224 		    *dstptr = result;
230 		Sgl_addition(left,right,/*into*/result);
231 		if(Sgl_isone_hidden(result))
233 		    *dstptr = result;
239 		/* need to normalize result */
240 	    	sign_save = Sgl_signextendedsign(result);
241 		Sgl_leftshiftby1(result);
242 		Sgl_normalize(result,result_exponent);
243 		Sgl_set_sign(result,/*using*/sign_save);
244                 Sgl_setwrapped_exponent(result,result_exponent,unfl);
245 		*dstptr = result;
248 	    *dstptr = result;
284 	Sgl_subtract_withextension(left,/*minus*/right,/*with*/extent,/*into*/result);
285 	if(Sgl_iszero_hidden(result))
288 	    /* A straightforward algorithm would now shift the result
297 	    sign_save = Sgl_signextendedsign(result);
298             Sgl_leftshiftby1_withextent(result,extent,result);
300             /* Need to check for a zero result.  The sign and exponent
304     	    if(Sgl_iszero(result))
307 		if(Is_rounding_mode(ROUNDMINUS)) Sgl_setone_sign(result);
308 		*dstptr = result;
313 	    if(Sgl_isone_hidden(result))
324 		    Sgl_set_sign(result,/*using*/sign_save);
340 	     * modified binary search.  We have already shifted the result
344 	    while(Sgl_iszero_hiddenhigh7mantissa(result))
346 		Sgl_leftshiftby8(result);
351 	    if(Sgl_iszero_hiddenhigh3mantissa(result))
354 		Sgl_leftshiftby4(result);
360 	    if((jumpsize = Sgl_hiddenhigh3mantissa(result)) > 7)
364 		Sgl_set_sign(result,/*using*/sign_save);
365 		Sgl_set_exponent(result,/*using*/result_exponent);
366 		*dstptr = result;
369 	    Sgl_sethigh4bits(result,/*using*/sign_save);
374 		    Sgl_leftshiftby3(result);
381 		    Sgl_leftshiftby2(result);
390 		    Sgl_leftshiftby1(result);
397 		Sgl_set_exponent(result,/*using*/result_exponent);
398 		*dstptr = result;
405 		Sgl_set_sign(result,sign_save);
406                 Sgl_setwrapped_exponent(result,result_exponent,unfl);
407 		*dstptr = result;
412 	     * Since we cannot get an inexact denormalized result,
415 	    Sgl_right_align(result,/*by*/(1-result_exponent),extent);
416 	    Sgl_clear_signexponent(result);
417 	    Sgl_set_sign(result,sign_save);
418 	    *dstptr = result;
426 	Sgl_addition(left,right,/*to*/result);
427 	if(Sgl_isone_hiddenoverflow(result))
430 	    Sgl_rightshiftby1_withextent(result,extent,extent);
431 	    Sgl_arithrightshiftby1(result);
436     /* Round the result.  If the extension is all zeros,then the result is
451 		  Sgl_isone_lowmantissa(result))
454 		    Sgl_increment(result);
460 	    if(Sgl_iszero_sign(result))
463 		Sgl_increment(result);
468 	    if(Sgl_isone_sign(result))
471 		Sgl_increment(result);
477 	if(Sgl_isone_hiddenoverflow(result)) result_exponent++;
484 	    Sgl_setwrapped_exponent(result,result_exponent,ovfl);
485 	    *dstptr = result;
496 	    Sgl_setoverflow(result);
499     else Sgl_set_exponent(result,result_exponent);
500     *dstptr = result;