Lines Matching defs:ay
277 double x, y, ax, ay, rx, ry, B, sqrt_A2my2, new_y;
284 ay = fabs(y);
304 if (ax > RECIP_EPSILON || ay > RECIP_EPSILON) {
320 if (ax < SQRT_6_EPSILON / 4 && ay < SQRT_6_EPSILON / 4)
323 do_hard_work(ax, ay, &rx, &B_is_usable, &B, &sqrt_A2my2, &new_y);
358 double x, y, ax, ay, rx, ry, B, sqrt_A2mx2, new_x;
368 ay = fabs(y);
388 if (ax > RECIP_EPSILON || ay > RECIP_EPSILON) {
405 if (ax < SQRT_6_EPSILON / 4 && ay < SQRT_6_EPSILON / 4)
408 do_hard_work(ay, ax, &ry, &B_is_usable, &B, &sqrt_A2mx2, &new_x);
458 double ax, ay, t;
463 ay = fabs(y);
464 if (ax < ay) {
466 ax = ay;
467 ay = t;
484 if (ax > QUARTER_SQRT_MAX || ay < SQRT_MIN)
487 return (CMPLX(log(ax * ax + ay * ay) / 2, atan2(y, x)));
569 double x, y, ax, ay, rx, ry;
574 ay = fabs(y);
600 if (ax > RECIP_EPSILON || ay > RECIP_EPSILON)
604 if (ax < SQRT_3_EPSILON / 2 && ay < SQRT_3_EPSILON / 2) {
614 if (ax == 1 && ay < DBL_EPSILON)
615 rx = (m_ln2 - log(ay)) / 2;
617 rx = log1p(4 * ax / sum_squares(ax - 1, ay)) / 4;
620 ry = atan2(2, -ay) / 2;
621 else if (ay < DBL_EPSILON)
622 ry = atan2(2 * ay, (1 - ax) * (1 + ax)) / 2;
624 ry = atan2(2 * ay, (1 - ax) * (1 + ax) - ay * ay) / 2;
