msun/src/catrig.c

148  * If B_is_usable is set to 0, sqrt_A2my2 = sqrt(A*A - y*y), and new_y = y.
149  * If returning sqrt_A2my2 has potential to result in an underflow, it is
154     double *sqrt_A2my2, double *new_y)
215 		*sqrt_A2my2 = A * (2 / DBL_EPSILON);
228 		 * sqrt_A2my2 = sqrt(Amy*(A+y))
235 			*sqrt_A2my2 = sqrt(x) * sqrt((A + y) / 2);
244 			*sqrt_A2my2 = sqrt(Amy * (A + y));
253 			*sqrt_A2my2 = x * (4 / DBL_EPSILON / DBL_EPSILON) * y /
261 			*sqrt_A2my2 = sqrt((1 - y) * (1 + y));
277 	double x, y, ax, ay, rx, ry, B, sqrt_A2my2, new_y;
323 	do_hard_work(ax, ay, &rx, &B_is_usable, &B, &sqrt_A2my2, &new_y);
327 		ry = atan2(new_y, sqrt_A2my2);