97       T yn = sqrt(z);  in ellint_rd_imp()  local
99       T y0 = yn;  in ellint_rd_imp()
103       while(fabs(xn - yn) >= 2.7 * tools::root_epsilon<T>() * fabs(xn))  in ellint_rd_imp()
105          T t = sqrt(xn * yn);  in ellint_rd_imp()
106          xn = (xn + yn) / 2;  in ellint_rd_imp()
107          yn = t;  in ellint_rd_imp()
109          sum += sum_pow * boost::math::pow<2>(xn - yn);  in ellint_rd_imp()
111       T RF = constants::pi<T>() / (xn + yn);  in ellint_rd_imp()
130    T yn = y;  in ellint_rd_imp()  local
145       ry = sqrt(yn);  in ellint_rd_imp()
151       yn = (yn + lambda) / 4;  in ellint_rd_imp()