Lines Matching refs:y
152 do_hard_work(double x, double y, double *rx, int *B_is_usable, double *B, in do_hard_work() argument
158 R = hypot(x, y + 1); /* |z+I| */ in do_hard_work()
159 S = hypot(x, y - 1); /* |z-I| */ in do_hard_work()
176 if (y == 1 && x < DBL_EPSILON * DBL_EPSILON / 128) { in do_hard_work()
182 } else if (x >= DBL_EPSILON * fabs(y - 1)) { in do_hard_work()
187 Am1 = f(x, 1 + y, R) + f(x, 1 - y, S); in do_hard_work()
189 } else if (y < 1) { in do_hard_work()
194 *rx = x / sqrt((1 - y) * (1 + y)); in do_hard_work()
199 *rx = log1p((y - 1) + sqrt((y - 1) * (y + 1))); in do_hard_work()
205 *new_y = y; in do_hard_work()
207 if (y < FOUR_SQRT_MIN) { in do_hard_work()
215 *new_y = y * (2 / DBL_EPSILON); in do_hard_work()
220 *B = y / A; in do_hard_work()
229 if (y == 1 && x < DBL_EPSILON / 128) { in do_hard_work()
234 *sqrt_A2my2 = sqrt(x) * sqrt((A + y) / 2); in do_hard_work()
235 } else if (x >= DBL_EPSILON * fabs(y - 1)) { in do_hard_work()
242 Amy = f(x, y + 1, R) + f(x, y - 1, S); in do_hard_work()
243 *sqrt_A2my2 = sqrt(Amy * (A + y)); in do_hard_work()
244 } else if (y > 1) { in do_hard_work()
252 *sqrt_A2my2 = x * (4 / DBL_EPSILON / DBL_EPSILON) * y / in do_hard_work()
253 sqrt((y + 1) * (y - 1)); in do_hard_work()
254 *new_y = y * (4 / DBL_EPSILON / DBL_EPSILON); in do_hard_work()
260 *sqrt_A2my2 = sqrt((1 - y) * (1 + y)); in do_hard_work()
276 double x, y, ax, ay, rx, ry, B, sqrt_A2my2, new_y; in casinh() local
281 y = cimag(z); in casinh()
283 ay = fabs(y); in casinh()
285 if (isnan(x) || isnan(y)) { in casinh()
288 return (CMPLX(x, y + y)); in casinh()
290 if (isinf(y)) in casinh()
291 return (CMPLX(y, x + x)); in casinh()
293 if (y == 0) in casinh()
294 return (CMPLX(x + x, y)); in casinh()
300 return (CMPLX(nan_mix(x, y), nan_mix(x, y))); in casinh()
309 return (CMPLX(copysign(creal(w), x), copysign(cimag(w), y))); in casinh()
313 if (x == 0 && y == 0) in casinh()
327 return (CMPLX(copysign(rx, x), copysign(ry, y))); in casinh()
357 double x, y, ax, ay, rx, ry, B, sqrt_A2mx2, new_x; in cacos() local
363 y = cimag(z); in cacos()
365 sy = signbit(y); in cacos()
367 ay = fabs(y); in cacos()
369 if (isnan(x) || isnan(y)) { in cacos()
372 return (CMPLX(y + y, -INFINITY)); in cacos()
374 if (isinf(y)) in cacos()
375 return (CMPLX(x + x, -y)); in cacos()
378 return (CMPLX(pio2_hi + pio2_lo, y + y)); in cacos()
384 return (CMPLX(nan_mix(x, y), nan_mix(x, y))); in cacos()
398 if (x == 1 && y == 0) in cacos()
399 return (CMPLX(0, -y)); in cacos()
405 return (CMPLX(pio2_hi - (x - pio2_lo), -y)); in cacos()
456 double x, y; in clog_for_large_values() local
460 y = cimag(z); in clog_for_large_values()
462 ay = fabs(y); in clog_for_large_values()
482 return (CMPLX(log(hypot(x / m_e, y / m_e)) + 1, atan2(y, x))); in clog_for_large_values()
489 return (CMPLX(log(hypot(x, y)), atan2(y, x))); in clog_for_large_values()
491 return (CMPLX(log(ax * ax + ay * ay) / 2, atan2(y, x))); in clog_for_large_values()
508 sum_squares(double x, double y) in sum_squares() argument
512 if (y < SQRT_MIN) in sum_squares()
515 return (x * x + y * y); in sum_squares()
528 real_part_reciprocal(double x, double y) in real_part_reciprocal() argument
540 GET_HIGH_WORD(hy, y); in real_part_reciprocal()
548 return (x / y / y); /* should avoid double div, but hard */ in real_part_reciprocal()
550 return (x / (x * x + y * y)); in real_part_reciprocal()
554 y *= scale; in real_part_reciprocal()
555 return (x / (x * x + y * y) * scale); in real_part_reciprocal()
573 double x, y, ax, ay, rx, ry; in catanh() local
576 y = cimag(z); in catanh()
578 ay = fabs(y); in catanh()
581 if (y == 0 && ax <= 1) in catanh()
582 return (CMPLX(atanh(x), y)); in catanh()
586 return (CMPLX(x, atan(y))); in catanh()
588 if (isnan(x) || isnan(y)) { in catanh()
591 return (CMPLX(copysign(0, x), y + y)); in catanh()
593 if (isinf(y)) in catanh()
595 copysign(pio2_hi + pio2_lo, y))); in catanh()
601 return (CMPLX(nan_mix(x, y), nan_mix(x, y))); in catanh()
605 return (CMPLX(real_part_reciprocal(x, y), in catanh()
606 copysign(pio2_hi + pio2_lo, y))); in catanh()
630 return (CMPLX(copysign(rx, x), copysign(ry, y))); in catanh()