parisc/math-emu/fmpyfadd.c

1 // SPDX-License-Identifier: GPL-2.0-or-later
3  * Linux/PA-RISC Project (http://www.parisc-linux.org/)
5  * Floating-point emulation code
6  *  Copyright (C) 2001 Hewlett-Packard (Paul Bame) <bame@debian.org>
15  *	Double Floating-point Multiply Fused Add
16  *	Double Floating-point Multiply Negate Fused Add
17  *	Single Floating-point Multiply Fused Add
18  *	Single Floating-point Multiply Negate Fused Add
41  *  Double Floating-point Multiply Fused Add
77 	mpy_exponent = Dbl_exponent(opnd1p1) + Dbl_exponent(opnd2p1) - DBL_BIAS;  in dbl_fmpyfadd()
369 	 * so lets call it a Double extended.  in dbl_fmpyfadd()
377 	for (count = DBL_P-1; count >= 0; count -= 4) {  in dbl_fmpyfadd()
436 			 * double precision format from double extended.  in dbl_fmpyfadd()
458 	 * Copy opnd3 to the double extended variable called right.  in dbl_fmpyfadd()
484 		diff_exponent = add_exponent - mpy_exponent;  in dbl_fmpyfadd()
488 		diff_exponent = mpy_exponent - add_exponent;  in dbl_fmpyfadd()
513 		 * bit and force a post normalization phase.  in dbl_fmpyfadd()
540 				/* Must have been "x-x" or "x+(-x)". */  in dbl_fmpyfadd()
546 			result_exponent--;  in dbl_fmpyfadd()
562 				result_exponent -= 8;  in dbl_fmpyfadd()
569 				result_exponent -= 4;  in dbl_fmpyfadd()
578 				result_exponent -= 3;  in dbl_fmpyfadd()
584 				result_exponent -= 2;  in dbl_fmpyfadd()
592 				result_exponent -= 1;  in dbl_fmpyfadd()
703  *  Double Floating-point Multiply Negate Fused Add
737 	mpy_exponent = Dbl_exponent(opnd1p1) + Dbl_exponent(opnd2p1) - DBL_BIAS;
1029 	 * so lets call it a Double extended.
1037 	for (count = DBL_P-1; count >= 0; count -= 4) {
1096 			 * double precision format from double extended.
1118 	 * Copy opnd3 to the double extended variable called right.
1144 		diff_exponent = add_exponent - mpy_exponent;
1148 		diff_exponent = mpy_exponent - add_exponent;
1173 		 * bit and force a post normalization phase.
1200 				/* Must have been "x-x" or "x+(-x)". */
1206 			result_exponent--;
1222 				result_exponent -= 8;
1229 				result_exponent -= 4;
1238 				result_exponent -= 3;
1244 				result_exponent -= 2;
1252 				result_exponent -= 1;
1361  *  Single Floating-point Multiply Fused Add
1394 	mpy_exponent = Sgl_exponent(opnd1) + Sgl_exponent(opnd2) - SGL_BIAS;
1685 	 * so lets call it a Double extended.
1693 	for (count = SGL_P-1; count >= 0; count -= 4) {
1750 			 * double precision format from double extended.
1771 	 * Copy opnd3 to the double extended variable called right.
1796 		diff_exponent = add_exponent - mpy_exponent;
1800 		diff_exponent = mpy_exponent - add_exponent;
1824 		 * bit and force a post normalization phase.
1849 				/* Must have been "x-x" or "x+(-x)". */
1855 			result_exponent--;
1871 				result_exponent -= 8;
1878 				result_exponent -= 4;
1886 				result_exponent -= 3;
1891 				result_exponent -= 2;
1898 				result_exponent -= 1;
2002  *  Single Floating-point Multiply Negate Fused Add
2036 	mpy_exponent = Sgl_exponent(opnd1) + Sgl_exponent(opnd2) - SGL_BIAS;
2327 	 * so lets call it a Double extended.
2335 	for (count = SGL_P-1; count >= 0; count -= 4) {
2392 			 * double precision format from double extended.
2413 	 * Copy opnd3 to the double extended variable called right.
2438 		diff_exponent = add_exponent - mpy_exponent;
2442 		diff_exponent = mpy_exponent - add_exponent;
2466 		 * bit and force a post normalization phase.
2491 				/* Must have been "x-x" or "x+(-x)". */
2497 			result_exponent--;
2513 				result_exponent -= 8;
2520 				result_exponent -= 4;
2528 				result_exponent -= 3;
2533 				result_exponent -= 2;
2540 				result_exponent -= 1;