Lines Matching defs:ay
161 float x, y, ax, ay, rx, ry, B, sqrt_A2my2, new_y;
168 ay = fabsf(y);
180 if (ax > RECIP_EPSILON || ay > RECIP_EPSILON) {
194 if (ax < SQRT_6_EPSILON / 4 && ay < SQRT_6_EPSILON / 4)
197 do_hard_work(ax, ay, &rx, &B_is_usable, &B, &sqrt_A2my2, &new_y);
216 float x, y, ax, ay, rx, ry, B, sqrt_A2mx2, new_x;
226 ay = fabsf(y);
238 if (ax > RECIP_EPSILON || ay > RECIP_EPSILON) {
252 if (ax < SQRT_6_EPSILON / 4 && ay < SQRT_6_EPSILON / 4)
255 do_hard_work(ay, ax, &ry, &B_is_usable, &B, &sqrt_A2mx2, &new_x);
294 float ax, ay, t;
299 ay = fabsf(y);
300 if (ax < ay) {
302 ax = ay;
303 ay = t;
310 if (ax > QUARTER_SQRT_MAX || ay < SQRT_MIN)
313 return (CMPLXF(logf(ax * ax + ay * ay) / 2, atan2f(y, x)));
354 float x, y, ax, ay, rx, ry;
359 ay = fabsf(y);
376 if (ax > RECIP_EPSILON || ay > RECIP_EPSILON)
380 if (ax < SQRT_3_EPSILON / 2 && ay < SQRT_3_EPSILON / 2) {
385 if (ax == 1 && ay < FLT_EPSILON)
386 rx = (m_ln2 - logf(ay)) / 2;
388 rx = log1pf(4 * ax / sum_squares(ax - 1, ay)) / 4;
391 ry = atan2f(2, -ay) / 2;
392 else if (ay < FLT_EPSILON)
393 ry = atan2f(2 * ay, (1 - ax) * (1 + ax)) / 2;
395 ry = atan2f(2 * ay, (1 - ax) * (1 + ax) - ay * ay) / 2;