Lines Matching refs:ay
418 double ay, h;
427 ay = fabs(z.imag);
428 if (z.real > CM_SQRT_LARGE_DOUBLE || ay > CM_SQRT_LARGE_DOUBLE) {
442 } else if (z.real == 1. && ay < CM_SQRT_DBL_MIN) {
444 if (ay == 0.) {
449 r.real = -log(sqrt(ay)/sqrt(hypot(ay, 2.)));
450 r.imag = copysign(atan2(2., -ay)/2, z.imag);
454 r.real = m_log1p(4.*z.real/((1-z.real)*(1-z.real) + ay*ay))/4.;
455 r.imag = -atan2(-2.*z.imag, (1-z.real)*(1+z.real) - ay*ay)/2.;
637 double ax, ay, am, an, h;
642 ay = fabs(z.imag);
644 if (ax > CM_LARGE_DOUBLE || ay > CM_LARGE_DOUBLE) {
645 r.real = log(hypot(ax/2., ay/2.)) + M_LN2;
646 } else if (ax < DBL_MIN && ay < DBL_MIN) {
647 if (ax > 0. || ay > 0.) {
648 /* catch cases where hypot(ax, ay) is subnormal */
650 ldexp(ay, DBL_MANT_DIG))) - DBL_MANT_DIG*M_LN2;
660 h = hypot(ax, ay);
662 am = ax > ay ? ax : ay; /* max(ax, ay) */
663 an = ax > ay ? ay : ax; /* min(ax, ay) */
819 double ax, ay;
830 ay = fabs(z.imag);
832 if (ax < DBL_MIN && ay < DBL_MIN && (ax > 0. || ay > 0.)) {
833 /* here we catch cases where hypot(ax, ay) is subnormal */
835 s = ldexp(sqrt(ax + hypot(ax, ldexp(ay, CM_SCALE_UP))),
839 s = 2.*sqrt(ax + hypot(ax, ay/8.));
841 d = ay/(2.*s);