411     double ay, h;  in cmath_atanh_impl()  local
420     ay = fabs(z.imag);  in cmath_atanh_impl()
421     if (z.real > CM_SQRT_LARGE_DOUBLE || ay > CM_SQRT_LARGE_DOUBLE) {  in cmath_atanh_impl()
435     } else if (z.real == 1. && ay < CM_SQRT_DBL_MIN) {  in cmath_atanh_impl()
437         if (ay == 0.) {  in cmath_atanh_impl()
442             r.real = -log(sqrt(ay)/sqrt(hypot(ay, 2.)));  in cmath_atanh_impl()
443             r.imag = copysign(atan2(2., -ay)/2, z.imag);  in cmath_atanh_impl()
447         r.real = m_log1p(4.*z.real/((1-z.real)*(1-z.real) + ay*ay))/4.;  in cmath_atanh_impl()
448         r.imag = -atan2(-2.*z.imag, (1-z.real)*(1+z.real) - ay*ay)/2.;  in cmath_atanh_impl()
630     double ax, ay, am, an, h;  in c_log()  local
635     ay = fabs(z.imag);  in c_log()
637     if (ax > CM_LARGE_DOUBLE || ay > CM_LARGE_DOUBLE) {  in c_log()
638         r.real = log(hypot(ax/2., ay/2.)) + M_LN2;  in c_log()
639     } else if (ax < DBL_MIN && ay < DBL_MIN) {  in c_log()
640         if (ax > 0. || ay > 0.) {  in c_log()
643                      ldexp(ay, DBL_MANT_DIG))) - DBL_MANT_DIG*M_LN2;  in c_log()
653         h = hypot(ax, ay);  in c_log()
655             am = ax > ay ? ax : ay;  /* max(ax, ay) */  in c_log()
656             an = ax > ay ? ay : ax;  /* min(ax, ay) */  in c_log()
812     double ax, ay;  in cmath_sqrt_impl()  local
823     ay = fabs(z.imag);  in cmath_sqrt_impl()
825     if (ax < DBL_MIN && ay < DBL_MIN && (ax > 0. || ay > 0.)) {  in cmath_sqrt_impl()
828         s = ldexp(sqrt(ax + hypot(ax, ldexp(ay, CM_SCALE_UP))),  in cmath_sqrt_impl()
832         s = 2.*sqrt(ax + hypot(ax, ay/8.));  in cmath_sqrt_impl()
834     d = ay/(2.*s);  in cmath_sqrt_impl()