Lines Matching refs:y
101 do_hard_work(float x, float y, float *rx, int *B_is_usable, float *B, in do_hard_work() argument
107 R = hypotf(x, y + 1); in do_hard_work()
108 S = hypotf(x, y - 1); in do_hard_work()
115 if (y == 1 && x < FLT_EPSILON * FLT_EPSILON / 128) { in do_hard_work()
117 } else if (x >= FLT_EPSILON * fabsf(y - 1)) { in do_hard_work()
118 Am1 = f(x, 1 + y, R) + f(x, 1 - y, S); in do_hard_work()
120 } else if (y < 1) { in do_hard_work()
121 *rx = x / sqrtf((1 - y) * (1 + y)); in do_hard_work()
123 *rx = log1pf((y - 1) + sqrtf((y - 1) * (y + 1))); in do_hard_work()
129 *new_y = y; in do_hard_work()
131 if (y < FOUR_SQRT_MIN) { in do_hard_work()
134 *new_y = y * (2 / FLT_EPSILON); in do_hard_work()
138 *B = y / A; in do_hard_work()
143 if (y == 1 && x < FLT_EPSILON / 128) { in do_hard_work()
144 *sqrt_A2my2 = sqrtf(x) * sqrtf((A + y) / 2); in do_hard_work()
145 } else if (x >= FLT_EPSILON * fabsf(y - 1)) { in do_hard_work()
146 Amy = f(x, y + 1, R) + f(x, y - 1, S); in do_hard_work()
147 *sqrt_A2my2 = sqrtf(Amy * (A + y)); in do_hard_work()
148 } else if (y > 1) { in do_hard_work()
149 *sqrt_A2my2 = x * (4 / FLT_EPSILON / FLT_EPSILON) * y / in do_hard_work()
150 sqrtf((y + 1) * (y - 1)); in do_hard_work()
151 *new_y = y * (4 / FLT_EPSILON / FLT_EPSILON); in do_hard_work()
153 *sqrt_A2my2 = sqrtf((1 - y) * (1 + y)); in do_hard_work()
161 float x, y, ax, ay, rx, ry, B, sqrt_A2my2, new_y; in casinhf() local
166 y = cimagf(z); in casinhf()
168 ay = fabsf(y); in casinhf()
170 if (isnan(x) || isnan(y)) { in casinhf()
172 return (CMPLXF(x, y + y)); in casinhf()
173 if (isinf(y)) in casinhf()
174 return (CMPLXF(y, x + x)); in casinhf()
175 if (y == 0) in casinhf()
176 return (CMPLXF(x + x, y)); in casinhf()
177 return (CMPLXF(x + 0.0L + (y + 0), x + 0.0L + (y + 0))); in casinhf()
186 copysignf(cimagf(w), y))); in casinhf()
189 if (x == 0 && y == 0) in casinhf()
202 return (CMPLXF(copysignf(rx, x), copysignf(ry, y))); in casinhf()
216 float x, y, ax, ay, rx, ry, B, sqrt_A2mx2, new_x; in cacosf() local
222 y = cimagf(z); in cacosf()
224 sy = signbit(y); in cacosf()
226 ay = fabsf(y); in cacosf()
228 if (isnan(x) || isnan(y)) { in cacosf()
230 return (CMPLXF(y + y, -INFINITY)); in cacosf()
231 if (isinf(y)) in cacosf()
232 return (CMPLXF(x + x, -y)); in cacosf()
234 return (CMPLXF(pio2_hi + pio2_lo, y + y)); in cacosf()
235 return (CMPLXF(x + 0.0L + (y + 0), x + 0.0L + (y + 0))); in cacosf()
247 if (x == 1 && y == 0) in cacosf()
248 return (CMPLXF(0, -y)); in cacosf()
253 return (CMPLXF(pio2_hi - (x - pio2_lo), -y)); in cacosf()
293 float x, y; in clog_for_large_values() local
297 y = cimagf(z); in clog_for_large_values()
299 ay = fabsf(y); in clog_for_large_values()
307 return (CMPLXF(logf(hypotf(x / m_e, y / m_e)) + 1, in clog_for_large_values()
308 atan2f(y, x))); in clog_for_large_values()
311 return (CMPLXF(logf(hypotf(x, y)), atan2f(y, x))); in clog_for_large_values()
313 return (CMPLXF(logf(ax * ax + ay * ay) / 2, atan2f(y, x))); in clog_for_large_values()
317 sum_squares(float x, float y) in sum_squares() argument
320 if (y < SQRT_MIN) in sum_squares()
323 return (x * x + y * y); in sum_squares()
327 real_part_reciprocal(float x, float y) in real_part_reciprocal() argument
335 GET_FLOAT_WORD(hy, y); in real_part_reciprocal()
342 return (x / y / y); in real_part_reciprocal()
344 return (x / (x * x + y * y)); in real_part_reciprocal()
347 y *= scale; in real_part_reciprocal()
348 return (x / (x * x + y * y) * scale); in real_part_reciprocal()
354 float x, y, ax, ay, rx, ry; in catanhf() local
357 y = cimagf(z); in catanhf()
359 ay = fabsf(y); in catanhf()
361 if (y == 0 && ax <= 1) in catanhf()
362 return (CMPLXF(atanhf(x), y)); in catanhf()
365 return (CMPLXF(x, atanf(y))); in catanhf()
367 if (isnan(x) || isnan(y)) { in catanhf()
369 return (CMPLXF(copysignf(0, x), y + y)); in catanhf()
370 if (isinf(y)) in catanhf()
372 copysignf(pio2_hi + pio2_lo, y))); in catanhf()
373 return (CMPLXF(x + 0.0L + (y + 0), x + 0.0L + (y + 0))); in catanhf()
377 return (CMPLXF(real_part_reciprocal(x, y), in catanhf()
378 copysignf(pio2_hi + pio2_lo, y))); in catanhf()
397 return (CMPLXF(copysignf(rx, x), copysignf(ry, y))); in catanhf()