87 do_hard_work(float x, float y, float *rx, int *B_is_usable, float *B,  in do_hard_work()  argument
93 	R = hypotf(x, y + 1);  in do_hard_work()
94 	S = hypotf(x, y - 1);  in do_hard_work()
101 		if (y == 1 && x < FLT_EPSILON * FLT_EPSILON / 128) {  in do_hard_work()
103 		} else if (x >= FLT_EPSILON * fabsf(y - 1)) {  in do_hard_work()
104 			Am1 = f(x, 1 + y, R) + f(x, 1 - y, S);  in do_hard_work()
106 		} else if (y < 1) {  in do_hard_work()
107 			*rx = x / sqrtf((1 - y) * (1 + y));  in do_hard_work()
109 			*rx = log1pf((y - 1) + sqrtf((y - 1) * (y + 1)));  in do_hard_work()
115 	*new_y = y;  in do_hard_work()
117 	if (y < FOUR_SQRT_MIN) {  in do_hard_work()
120 		*new_y = y * (2 / FLT_EPSILON);  in do_hard_work()
124 	*B = y / A;  in do_hard_work()
129 		if (y == 1 && x < FLT_EPSILON / 128) {  in do_hard_work()
130 			*sqrt_A2my2 = sqrtf(x) * sqrtf((A + y) / 2);  in do_hard_work()
131 		} else if (x >= FLT_EPSILON * fabsf(y - 1)) {  in do_hard_work()
132 			Amy = f(x, y + 1, R) + f(x, y - 1, S);  in do_hard_work()
133 			*sqrt_A2my2 = sqrtf(Amy * (A + y));  in do_hard_work()
134 		} else if (y > 1) {  in do_hard_work()
135 			*sqrt_A2my2 = x * (4 / FLT_EPSILON / FLT_EPSILON) * y /  in do_hard_work()
136 			    sqrtf((y + 1) * (y - 1));  in do_hard_work()
137 			*new_y = y * (4 / FLT_EPSILON / FLT_EPSILON);  in do_hard_work()
139 			*sqrt_A2my2 = sqrtf((1 - y) * (1 + y));  in do_hard_work()
147 	float x, y, ax, ay, rx, ry, B, sqrt_A2my2, new_y;  in casinhf()  local
152 	y = cimagf(z);  in casinhf()
154 	ay = fabsf(y);  in casinhf()
156 	if (isnan(x) || isnan(y)) {  in casinhf()
158 			return (CMPLXF(x, y + y));  in casinhf()
159 		if (isinf(y))  in casinhf()
160 			return (CMPLXF(y, x + x));  in casinhf()
161 		if (y == 0)  in casinhf()
162 			return (CMPLXF(x + x, y));  in casinhf()
163 		return (CMPLXF(nan_mix(x, y), nan_mix(x, y)));  in casinhf()
172 		    copysignf(cimagf(w), y)));  in casinhf()
175 	if (x == 0 && y == 0)  in casinhf()
188 	return (CMPLXF(copysignf(rx, x), copysignf(ry, y)));  in casinhf()
202 	float x, y, ax, ay, rx, ry, B, sqrt_A2mx2, new_x;  in cacosf()  local
208 	y = cimagf(z);  in cacosf()
210 	sy = signbit(y);  in cacosf()
212 	ay = fabsf(y);  in cacosf()
214 	if (isnan(x) || isnan(y)) {  in cacosf()
216 			return (CMPLXF(y + y, -INFINITY));  in cacosf()
217 		if (isinf(y))  in cacosf()
218 			return (CMPLXF(x + x, -y));  in cacosf()
220 			return (CMPLXF(pio2_hi + pio2_lo, y + y));  in cacosf()
221 		return (CMPLXF(nan_mix(x, y), nan_mix(x, y)));  in cacosf()
233 	if (x == 1 && y == 0)  in cacosf()
234 		return (CMPLXF(0, -y));  in cacosf()
239 		return (CMPLXF(pio2_hi - (x - pio2_lo), -y));  in cacosf()
279 	float x, y;  in clog_for_large_values()  local
283 	y = cimagf(z);  in clog_for_large_values()
285 	ay = fabsf(y);  in clog_for_large_values()
293 		return (CMPLXF(logf(hypotf(x / m_e, y / m_e)) + 1,  in clog_for_large_values()
294 		    atan2f(y, x)));  in clog_for_large_values()
297 		return (CMPLXF(logf(hypotf(x, y)), atan2f(y, x)));  in clog_for_large_values()
299 	return (CMPLXF(logf(ax * ax + ay * ay) / 2, atan2f(y, x)));  in clog_for_large_values()
303 sum_squares(float x, float y)  in sum_squares()  argument
306 	if (y < SQRT_MIN)  in sum_squares()
309 	return (x * x + y * y);  in sum_squares()
313 real_part_reciprocal(float x, float y)  in real_part_reciprocal()  argument
321 	GET_FLOAT_WORD(hy, y);  in real_part_reciprocal()
328 		return (x / y / y);  in real_part_reciprocal()
330 		return (x / (x * x + y * y));  in real_part_reciprocal()
333 	y *= scale;  in real_part_reciprocal()
334 	return (x / (x * x + y * y) * scale);  in real_part_reciprocal()
340 	float x, y, ax, ay, rx, ry;  in catanhf()  local
343 	y = cimagf(z);  in catanhf()
345 	ay = fabsf(y);  in catanhf()
347 	if (y == 0 && ax <= 1)  in catanhf()
348 		return (CMPLXF(atanhf(x), y));  in catanhf()
351 		return (CMPLXF(x, atanf(y)));  in catanhf()
353 	if (isnan(x) || isnan(y)) {  in catanhf()
355 			return (CMPLXF(copysignf(0, x), y + y));  in catanhf()
356 		if (isinf(y))  in catanhf()
358 			    copysignf(pio2_hi + pio2_lo, y)));  in catanhf()
359 		return (CMPLXF(nan_mix(x, y), nan_mix(x, y)));  in catanhf()
363 		return (CMPLXF(real_part_reciprocal(x, y),  in catanhf()
364 		    copysignf(pio2_hi + pio2_lo, y)));  in catanhf()
383 	return (CMPLXF(copysignf(rx, x), copysignf(ry, y)));  in catanhf()