1 // SPDX-License-Identifier: GPL-2.0
2 /*---------------------------------------------------------------------------+
9  |                  E-mail   billm@suburbia.net                              |
12  +---------------------------------------------------------------------------*/
24 	0x0000000000000000LL,	/* Dummy (not for - 1.0) */
49 /*--- poly_atan() -----------------------------------------------------------+
51  +---------------------------------------------------------------------------*/
55 	u_char transformed, inverted, sign1, sign2;  in poly_atan()  local
66 		/* This gives non-compatible stack contents... */  in poly_atan()
71 		exponent -= exponent(st1_ptr);  in poly_atan()
73 		/* This gives non-compatible stack contents... */  in poly_atan()
75 		exponent -= exponent16(st1_ptr);  in poly_atan()
79 			       ((st0_ptr->sigh < st1_ptr->sigh) ||  in poly_atan()
80 				((st0_ptr->sigh == st1_ptr->sigh) &&  in poly_atan()
81 				 (st0_ptr->sigl < st1_ptr->sigl))))) {  in poly_atan()
82 		inverted = 1;  in poly_atan()
87 		inverted = 0;  in poly_atan()
88 		exponent = -exponent;  in poly_atan()
96 	if ((exponent >= -1)  in poly_atan()
97 	    || ((exponent == -2) && (argSignif.msw > 0xd413ccd0))) {  in poly_atan()
98 		/* The argument is greater than sqrt(2)-1 (=0.414213562...) */  in poly_atan()
111 			argSignif.msw = 0;	/* Make the transformed arg -> 0.0 */  in poly_atan()
116 			if (exponent < -1)  in poly_atan()
117 				shr_Xsig(&Numer, -1 - exponent);  in poly_atan()
120 			shr_Xsig(&Denom, -exponent);  in poly_atan()
125 			exponent = -1 + norm_Xsig(&argSignif);  in poly_atan()
144 	shr_Xsig(&argSq, 2 * (-1 - exponent - 1));  in poly_atan()
145 	shr_Xsig(&argSqSq, 4 * (-1 - exponent - 1));  in poly_atan()
153 			oddplterms, HIPOWERop - 1);  in poly_atan()
157 			HIPOWERon - 1);  in poly_atan()
162 	shr_Xsig(&accumulatore, 1 + 2 * (-1 - exponent));  in poly_atan()
175 		/* compute pi/4 - accumulator */  in poly_atan()
176 		shr_Xsig(&accumulator, -1 - exponent);  in poly_atan()
179 		exponent = -1;  in poly_atan()
182 	if (inverted) {  in poly_atan()
183 		/* compute pi/2 - accumulator */  in poly_atan()
184 		shr_Xsig(&accumulator, -exponent);  in poly_atan()
191 		/* compute pi - accumulator */  in poly_atan()
192 		shr_Xsig(&accumulator, 1 - exponent);  in poly_atan()