277 	double x, y, ax, ay, rx, ry, B, sqrt_A2my2, new_y;  in casinh()  local
284 	ay = fabs(y);  in casinh()
304 	if (ax > RECIP_EPSILON || ay > RECIP_EPSILON) {  in casinh()
320 	if (ax < SQRT_6_EPSILON / 4 && ay < SQRT_6_EPSILON / 4)  in casinh()
323 	do_hard_work(ax, ay, &rx, &B_is_usable, &B, &sqrt_A2my2, &new_y);  in casinh()
358 	double x, y, ax, ay, rx, ry, B, sqrt_A2mx2, new_x;  in cacos()  local
368 	ay = fabs(y);  in cacos()
388 	if (ax > RECIP_EPSILON || ay > RECIP_EPSILON) {  in cacos()
405 	if (ax < SQRT_6_EPSILON / 4 && ay < SQRT_6_EPSILON / 4)  in cacos()
408 	do_hard_work(ay, ax, &ry, &B_is_usable, &B, &sqrt_A2mx2, &new_x);  in cacos()
458 	double ax, ay, t;  in clog_for_large_values()  local
463 	ay = fabs(y);  in clog_for_large_values()
464 	if (ax < ay) {  in clog_for_large_values()
466 		ax = ay;  in clog_for_large_values()
467 		ay = t;  in clog_for_large_values()
489 	if (ax > QUARTER_SQRT_MAX || ay < SQRT_MIN)  in clog_for_large_values()
492 	return (CMPLX(log(ax * ax + ay * ay) / 2, atan2(y, x)));  in clog_for_large_values()
574 	double x, y, ax, ay, rx, ry;  in catanh()  local
579 	ay = fabs(y);  in catanh()
605 	if (ax > RECIP_EPSILON || ay > RECIP_EPSILON)  in catanh()
609 	if (ax < SQRT_3_EPSILON / 2 && ay < SQRT_3_EPSILON / 2) {  in catanh()
619 	if (ax == 1 && ay < DBL_EPSILON)  in catanh()
620 		rx = (m_ln2 - log(ay)) / 2;  in catanh()
622 		rx = log1p(4 * ax / sum_squares(ax - 1, ay)) / 4;  in catanh()
625 		ry = atan2(2, -ay) / 2;  in catanh()
626 	else if (ay < DBL_EPSILON)  in catanh()
627 		ry = atan2(2 * ay, (1 - ax) * (1 + ax)) / 2;  in catanh()
629 		ry = atan2(2 * ay, (1 - ax) * (1 + ax) - ay * ay) / 2;  in catanh()