parisc/math-emu/sfsub.c

42     register unsigned int left, right, result, extent;
73                     Sgl_makequietnan(result);
74 		    *dstptr = result;
186 		/* Left is not a zero and must be the result.  Trapped
210 	    /* Both operands are denormalized.  The result must be exact
215 		Sgl_subtract(left,/*minus*/right,/*into*/result);
216 		if(Sgl_iszero_mantissa(result))
220 			Sgl_setone_sign(result);
224 			Sgl_setzero_sign(result);
226 		    *dstptr = result;
232 		Sgl_addition(left,right,/*into*/result);
233 		if(Sgl_isone_hidden(result))
235 		    *dstptr = result;
241 		/* need to normalize result */
242 	    	sign_save = Sgl_signextendedsign(result);
243 		Sgl_leftshiftby1(result);
244 		Sgl_normalize(result,result_exponent);
245 		Sgl_set_sign(result,/*using*/sign_save);
246                 Sgl_setwrapped_exponent(result,result_exponent,unfl);
247 		*dstptr = result;
251 	    *dstptr = result;
287 	Sgl_subtract_withextension(left,/*minus*/right,/*with*/extent,/*into*/result);
288 	if(Sgl_iszero_hidden(result))
291 	    /* A straightforward algorithm would now shift the result
300 	    sign_save = Sgl_signextendedsign(result);
301             Sgl_leftshiftby1_withextent(result,extent,result);
303             /* Need to check for a zero result.  The sign and exponent
307     	    if(Sgl_iszero(result))
310 		if(Is_rounding_mode(ROUNDMINUS)) Sgl_setone_sign(result);
311 		*dstptr = result;
316 	    if(Sgl_isone_hidden(result))
327 		    Sgl_set_sign(result,/*using*/sign_save);
343 	     * modified binary search.  We have already shifted the result
347 	    while(Sgl_iszero_hiddenhigh7mantissa(result))
349 		Sgl_leftshiftby8(result);
354 	    if(Sgl_iszero_hiddenhigh3mantissa(result))
357 		Sgl_leftshiftby4(result);
363 	    if((jumpsize = Sgl_hiddenhigh3mantissa(result)) > 7)
367 		Sgl_set_sign(result,/*using*/sign_save);
368 		Sgl_set_exponent(result,/*using*/result_exponent);
369 		*dstptr = result;
372 	    Sgl_sethigh4bits(result,/*using*/sign_save);
377 		    Sgl_leftshiftby3(result);
384 		    Sgl_leftshiftby2(result);
393 		    Sgl_leftshiftby1(result);
400 		Sgl_set_exponent(result,/*using*/result_exponent);
401 		*dstptr = result;	/* Sign bit is already set */
408 		Sgl_set_sign(result,sign_save);
409                 Sgl_setwrapped_exponent(result,result_exponent,unfl);
410 		*dstptr = result;
415 	     * Since we cannot get an inexact denormalized result,
418 	    Sgl_right_align(result,/*by*/(1-result_exponent),extent);
419 	    Sgl_clear_signexponent(result);
420 	    Sgl_set_sign(result,sign_save);
421 	    *dstptr = result;
429 	Sgl_addition(left,right,/*to*/result);
430 	if(Sgl_isone_hiddenoverflow(result))
433 	    Sgl_rightshiftby1_withextent(result,extent,extent);
434 	    Sgl_arithrightshiftby1(result);
439     /* Round the result.  If the extension is all zeros,then the result is
454 		  Sgl_isone_lowmantissa(result))
457 		    Sgl_increment(result);
463 	    if(Sgl_iszero_sign(result))
466 		Sgl_increment(result);
471 	    if(Sgl_isone_sign(result))
474 		Sgl_increment(result);
480 	if(Sgl_isone_hiddenoverflow(result)) result_exponent++;
487 	    Sgl_setwrapped_exponent(result,result_exponent,ovfl);
488 	    *dstptr = result;
499 	    Sgl_setoverflow(result);
502     else Sgl_set_exponent(result,result_exponent);
503     *dstptr = result;