155 do_hard_work(double x, double y, double *rx, int *B_is_usable, double *B,  in do_hard_work()  argument
161 	R = hypot(x, y + 1);		/* |z+I| */  in do_hard_work()
162 	S = hypot(x, y - 1);		/* |z-I| */  in do_hard_work()
179 		if (y == 1 && x < DBL_EPSILON * DBL_EPSILON / 128) {  in do_hard_work()
185 		} else if (x >= DBL_EPSILON * fabs(y - 1)) {  in do_hard_work()
190 			Am1 = f(x, 1 + y, R) + f(x, 1 - y, S);  in do_hard_work()
192 		} else if (y < 1) {  in do_hard_work()
197 			*rx = x / sqrt((1 - y) * (1 + y));  in do_hard_work()
202 			*rx = log1p((y - 1) + sqrt((y - 1) * (y + 1)));  in do_hard_work()
208 	*new_y = y;  in do_hard_work()
210 	if (y < FOUR_SQRT_MIN) {  in do_hard_work()
218 		*new_y = y * (2 / DBL_EPSILON);  in do_hard_work()
223 	*B = y / A;  in do_hard_work()
232 		if (y == 1 && x < DBL_EPSILON / 128) {  in do_hard_work()
237 			*sqrt_A2my2 = sqrt(x) * sqrt((A + y) / 2);  in do_hard_work()
238 		} else if (x >= DBL_EPSILON * fabs(y - 1)) {  in do_hard_work()
245 			Amy = f(x, y + 1, R) + f(x, y - 1, S);  in do_hard_work()
246 			*sqrt_A2my2 = sqrt(Amy * (A + y));  in do_hard_work()
247 		} else if (y > 1) {  in do_hard_work()
255 			*sqrt_A2my2 = x * (4 / DBL_EPSILON / DBL_EPSILON) * y /  in do_hard_work()
256 			    sqrt((y + 1) * (y - 1));  in do_hard_work()
257 			*new_y = y * (4 / DBL_EPSILON / DBL_EPSILON);  in do_hard_work()
263 			*sqrt_A2my2 = sqrt((1 - y) * (1 + y));  in do_hard_work()
279 	double x, y, ax, ay, rx, ry, B, sqrt_A2my2, new_y;  in casinh()  local
284 	y = cimag(z);  in casinh()
286 	ay = fabs(y);  in casinh()
288 	if (isnan(x) || isnan(y)) {  in casinh()
291 			return (CMPLX(x, y + y));  in casinh()
293 		if (isinf(y))  in casinh()
294 			return (CMPLX(y, x + x));  in casinh()
296 		if (y == 0)  in casinh()
297 			return (CMPLX(x + x, y));  in casinh()
303 		return (CMPLX(nan_mix(x, y), nan_mix(x, y)));  in casinh()
312 		return (CMPLX(copysign(creal(w), x), copysign(cimag(w), y)));  in casinh()
316 	if (x == 0 && y == 0)  in casinh()
330 	return (CMPLX(copysign(rx, x), copysign(ry, y)));  in casinh()
360 	double x, y, ax, ay, rx, ry, B, sqrt_A2mx2, new_x;  in cacos()  local
366 	y = cimag(z);  in cacos()
368 	sy = signbit(y);  in cacos()
370 	ay = fabs(y);  in cacos()
372 	if (isnan(x) || isnan(y)) {  in cacos()
375 			return (CMPLX(y + y, -INFINITY));  in cacos()
377 		if (isinf(y))  in cacos()
378 			return (CMPLX(x + x, -y));  in cacos()
381 			return (CMPLX(pio2_hi + pio2_lo, y + y));  in cacos()
387 		return (CMPLX(nan_mix(x, y), nan_mix(x, y)));  in cacos()
401 	if (x == 1 && y == 0)  in cacos()
402 		return (CMPLX(0, -y));  in cacos()
408 		return (CMPLX(pio2_hi - (x - pio2_lo), -y));  in cacos()
459 	double x, y;  in clog_for_large_values()  local
463 	y = cimag(z);  in clog_for_large_values()
465 	ay = fabs(y);  in clog_for_large_values()
485 		return (CMPLX(log(hypot(x / m_e, y / m_e)) + 1, atan2(y, x)));  in clog_for_large_values()
492 		return (CMPLX(log(hypot(x, y)), atan2(y, x)));  in clog_for_large_values()
494 	return (CMPLX(log(ax * ax + ay * ay) / 2, atan2(y, x)));  in clog_for_large_values()
511 sum_squares(double x, double y)  in sum_squares()  argument
515 	if (y < SQRT_MIN)  in sum_squares()
518 	return (x * x + y * y);  in sum_squares()
531 real_part_reciprocal(double x, double y)  in real_part_reciprocal()  argument
543 	GET_HIGH_WORD(hy, y);  in real_part_reciprocal()
551 		return (x / y / y);	/* should avoid double div, but hard */  in real_part_reciprocal()
553 		return (x / (x * x + y * y));  in real_part_reciprocal()
557 	y *= scale;  in real_part_reciprocal()
558 	return (x / (x * x + y * y) * scale);  in real_part_reciprocal()
576 	double x, y, ax, ay, rx, ry;  in catanh()  local
579 	y = cimag(z);  in catanh()
581 	ay = fabs(y);  in catanh()
584 	if (y == 0 && ax <= 1)  in catanh()
585 		return (CMPLX(atanh(x), y));  in catanh()
589 		return (CMPLX(x, atan(y)));  in catanh()
591 	if (isnan(x) || isnan(y)) {  in catanh()
594 			return (CMPLX(copysign(0, x), y + y));  in catanh()
596 		if (isinf(y))  in catanh()
598 			    copysign(pio2_hi + pio2_lo, y)));  in catanh()
604 		return (CMPLX(nan_mix(x, y), nan_mix(x, y)));  in catanh()
608 		return (CMPLX(real_part_reciprocal(x, y),  in catanh()
609 		    copysign(pio2_hi + pio2_lo, y)));  in catanh()
633 	return (CMPLX(copysign(rx, x), copysign(ry, y)));  in catanh()