Lines Matching refs:ay
169 long double x, y, ax, ay, rx, ry, B, sqrt_A2my2, new_y;
176 ay = fabsl(y);
188 if (ax > RECIP_EPSILON || ay > RECIP_EPSILON) {
202 if (ax < SQRT_6_EPSILON / 4 && ay < SQRT_6_EPSILON / 4)
205 do_hard_work(ax, ay, &rx, &B_is_usable, &B, &sqrt_A2my2, &new_y);
225 long double x, y, ax, ay, rx, ry, B, sqrt_A2mx2, new_x;
235 ay = fabsl(y);
247 if (ax > RECIP_EPSILON || ay > RECIP_EPSILON) {
261 if (ax < SQRT_6_EPSILON / 4 && ay < SQRT_6_EPSILON / 4)
264 do_hard_work(ay, ax, &ry, &B_is_usable, &B, &sqrt_A2mx2, &new_x);
303 long double ax, ay, t;
308 ay = fabsl(y);
309 if (ax < ay) {
311 ax = ay;
312 ay = t;
319 if (ax > QUARTER_SQRT_MAX || ay < SQRT_MIN)
322 return (CMPLXL(logl(ax * ax + ay * ay) / 2, atan2l(y, x)));
364 long double x, y, ax, ay, rx, ry;
369 ay = fabsl(y);
386 if (ax > RECIP_EPSILON || ay > RECIP_EPSILON)
390 if (ax < SQRT_3_EPSILON / 2 && ay < SQRT_3_EPSILON / 2) {
395 if (ax == 1 && ay < LDBL_EPSILON)
396 rx = (m_ln2 - logl(ay)) / 2;
398 rx = log1pl(4 * ax / sum_squares(ax - 1, ay)) / 4;
401 ry = atan2l(2, -ay) / 2;
402 else if (ay < LDBL_EPSILON)
403 ry = atan2l(2 * ay, (1 - ax) * (1 + ax)) / 2;
405 ry = atan2l(2 * ay, (1 - ax) * (1 + ax) - ay * ay) / 2;