310     double ay, h;  in c_atanh()  local
319     ay = fabs(z.imag);  in c_atanh()
320     if (z.real > CM_SQRT_LARGE_DOUBLE || ay > CM_SQRT_LARGE_DOUBLE) {  in c_atanh()
334     } else if (z.real == 1. && ay < CM_SQRT_DBL_MIN) {  in c_atanh()
336         if (ay == 0.) {  in c_atanh()
341             r.real = -log(sqrt(ay)/sqrt(hypot(ay, 2.)));  in c_atanh()
342             r.imag = copysign(atan2(2., -ay)/2, z.imag);  in c_atanh()
346         r.real = m_log1p(4.*z.real/((1-z.real)*(1-z.real) + ay*ay))/4.;  in c_atanh()
347         r.imag = -atan2(-2.*z.imag, (1-z.real)*(1+z.real) - ay*ay)/2.;  in c_atanh()
529     double ax, ay, am, an, h;  in c_log()  local
534     ay = fabs(z.imag);  in c_log()
536     if (ax > CM_LARGE_DOUBLE || ay > CM_LARGE_DOUBLE) {  in c_log()
537         r.real = log(hypot(ax/2., ay/2.)) + M_LN2;  in c_log()
538     } else if (ax < DBL_MIN && ay < DBL_MIN) {  in c_log()
539         if (ax > 0. || ay > 0.) {  in c_log()
542                      ldexp(ay, DBL_MANT_DIG))) - DBL_MANT_DIG*M_LN2;  in c_log()
552         h = hypot(ax, ay);  in c_log()
554             am = ax > ay ? ax : ay;  /* max(ax, ay) */  in c_log()
555             an = ax > ay ? ay : ax;  /* min(ax, ay) */  in c_log()
698     double ax, ay;  in c_sqrt()  local
709     ay = fabs(z.imag);  in c_sqrt()
711     if (ax < DBL_MIN && ay < DBL_MIN && (ax > 0. || ay > 0.)) {  in c_sqrt()
714         s = ldexp(sqrt(ax + hypot(ax, ldexp(ay, CM_SCALE_UP))),  in c_sqrt()
718         s = 2.*sqrt(ax + hypot(ax, ay/8.));  in c_sqrt()
720     d = ay/(2.*s);  in c_sqrt()