53 #define isinf(x)	(fabsf(x) == INFINITY)  argument
55 #define isnan(x)	((x) != (x))  argument
58 #define signbit(x)	(__builtin_signbitf(x))  argument
90 do_hard_work(float x, float y, float *rx, int *B_is_usable, float *B,  in do_hard_work()  argument
96 	R = hypotf(x, y + 1);  in do_hard_work()
97 	S = hypotf(x, y - 1);  in do_hard_work()
104 		if (y == 1 && x < FLT_EPSILON * FLT_EPSILON / 128) {  in do_hard_work()
105 			*rx = sqrtf(x);  in do_hard_work()
106 		} else if (x >= FLT_EPSILON * fabsf(y - 1)) {  in do_hard_work()
107 			Am1 = f(x, 1 + y, R) + f(x, 1 - y, S);  in do_hard_work()
110 			*rx = x / sqrtf((1 - y) * (1 + y));  in do_hard_work()
132 		if (y == 1 && x < FLT_EPSILON / 128) {  in do_hard_work()
133 			*sqrt_A2my2 = sqrtf(x) * sqrtf((A + y) / 2);  in do_hard_work()
134 		} else if (x >= FLT_EPSILON * fabsf(y - 1)) {  in do_hard_work()
135 			Amy = f(x, y + 1, R) + f(x, y - 1, S);  in do_hard_work()
138 			*sqrt_A2my2 = x * (4 / FLT_EPSILON / FLT_EPSILON) * y /  in do_hard_work()
150 	float x, y, ax, ay, rx, ry, B, sqrt_A2my2, new_y;  in casinhf()  local
154 	x = crealf(z);  in casinhf()
156 	ax = fabsf(x);  in casinhf()
159 	if (isnan(x) || isnan(y)) {  in casinhf()
160 		if (isinf(x))  in casinhf()
161 			return (CMPLXF(x, y + y));  in casinhf()
163 			return (CMPLXF(y, x + x));  in casinhf()
165 			return (CMPLXF(x + x, y));  in casinhf()
166 		return (CMPLXF(nan_mix(x, y), nan_mix(x, y)));  in casinhf()
170 		if (signbit(x) == 0)  in casinhf()
174 		return (CMPLXF(copysignf(crealf(w), x),  in casinhf()
178 	if (x == 0 && y == 0)  in casinhf()
191 	return (CMPLXF(copysignf(rx, x), copysignf(ry, y)));  in casinhf()
205 	float x, y, ax, ay, rx, ry, B, sqrt_A2mx2, new_x;  in cacosf()  local
210 	x = crealf(z);  in cacosf()
212 	sx = signbit(x);  in cacosf()
214 	ax = fabsf(x);  in cacosf()
217 	if (isnan(x) || isnan(y)) {  in cacosf()
218 		if (isinf(x))  in cacosf()
221 			return (CMPLXF(x + x, -y));  in cacosf()
222 		if (x == 0)  in cacosf()
224 		return (CMPLXF(nan_mix(x, y), nan_mix(x, y)));  in cacosf()
236 	if (x == 1 && y == 0)  in cacosf()
242 		return (CMPLXF(pio2_hi - (x - pio2_lo), -y));  in cacosf()
282 	float x, y;  in clog_for_large_values()  local
285 	x = crealf(z);  in clog_for_large_values()
287 	ax = fabsf(x);  in clog_for_large_values()
296 		return (CMPLXF(logf(hypotf(x / m_e, y / m_e)) + 1,  in clog_for_large_values()
297 		    atan2f(y, x)));  in clog_for_large_values()
300 		return (CMPLXF(logf(hypotf(x, y)), atan2f(y, x)));  in clog_for_large_values()
302 	return (CMPLXF(logf(ax * ax + ay * ay) / 2, atan2f(y, x)));  in clog_for_large_values()
306 sum_squares(float x, float y)  in sum_squares()  argument
310 		return (x * x);  in sum_squares()
312 	return (x * x + y * y);  in sum_squares()
316 real_part_reciprocal(float x, float y)  in real_part_reciprocal()  argument
322 	GET_FLOAT_WORD(hx, x);  in real_part_reciprocal()
328 	if (ix - iy >= CUTOFF << 23 || isinf(x))  in real_part_reciprocal()
329 		return (1 / x);  in real_part_reciprocal()
331 		return (x / y / y);  in real_part_reciprocal()
333 		return (x / (x * x + y * y));  in real_part_reciprocal()
335 	x *= scale;  in real_part_reciprocal()
337 	return (x / (x * x + y * y) * scale);  in real_part_reciprocal()
343 	float x, y, ax, ay, rx, ry;  in catanhf()  local
345 	x = crealf(z);  in catanhf()
347 	ax = fabsf(x);  in catanhf()
351 		return (CMPLXF(atanhf(x), y));  in catanhf()
353 	if (x == 0)  in catanhf()
354 		return (CMPLXF(x, atanf(y)));  in catanhf()
356 	if (isnan(x) || isnan(y)) {  in catanhf()
357 		if (isinf(x))  in catanhf()
358 			return (CMPLXF(copysignf(0, x), y + y));  in catanhf()
360 			return (CMPLXF(copysignf(0, x),  in catanhf()
362 		return (CMPLXF(nan_mix(x, y), nan_mix(x, y)));  in catanhf()
366 		return (CMPLXF(real_part_reciprocal(x, y),  in catanhf()
386 	return (CMPLXF(copysignf(rx, x), copysignf(ry, y)));  in catanhf()