Method solve_limits

src/PolyMath.cpp:190–198

Uses the Brent solver to find the roots of p(x_in,y_in)-z_in @param res Poly2DResidual object to calculate residuals and derivatives @param min double value that represents the minimum value @param max double value that represents the maximum value

/// @param min double value that represents the minimum value
/// @param max double value that represents the maximum value
double Polynomial2D::solve_limits(Poly2DResidual* res, const double& min, const double& max) {
    if (do_debug()) std::cout << format("Called solve_limits with: min=%f and max=%f", min, max) << '\n';
    double macheps = DBL_EPSILON;
    double tol = DBL_EPSILON * 1e3;
    int maxiter = 10;
    double result = Brent(res, min, max, macheps, tol, maxiter);
    if (this->do_debug()) std::cout << "Brent solver message: " << res->errstring << '\n';
    return result;
}

/// Uses the Newton solver to find the roots of p(x_in,y_in)-z_in
/// @param res Poly2DResidual object to calculate residuals and derivatives