ALSolver.hpp

Go to the documentation of this file.

#pragma once
 
#include <polyfem/solver/NLProblem.hpp>
#include <polysolve/nonlinear/Solver.hpp>
#include <polyfem/solver/forms/BCLagrangianForm.hpp>
#include <polyfem/solver/forms/BCPenaltyForm.hpp>
#include <polyfem/Common.hpp>
 
#include <Eigen/Core>
 
#include <functional>
#include <vector>
 
namespace polyfem::solver
{
    class ALSolver
    {
        using NLSolver = polysolve::nonlinear::Solver;
 
    public:
        ALSolver(
            std::shared_ptr<BCLagrangianForm> lagr_form,
            std::shared_ptr<BCPenaltyForm> pen_form,
            const double initial_al_weight,
            const double scaling,
            const double max_al_weight,
            const double eta_tol,
            const std::function<void(const Eigen::VectorXd &)> &update_barrier_stiffness);
 
        void solve_al(std::shared_ptr<NLSolver> nl_solver, NLProblem &nl_problem, Eigen::MatrixXd &sol);
        void solve_reduced(std::shared_ptr<NLSolver> nl_solver, NLProblem &nl_problem, Eigen::MatrixXd &sol);
 
        std::function<void(const double)> post_subsolve = [](const double) {};
 
    protected:
        void set_al_weight(NLProblem &nl_problem, const Eigen::VectorXd &x, const double weight);
 
        std::shared_ptr<BCLagrangianForm> lagr_form;
        std::shared_ptr<BCPenaltyForm> pen_form;
        const double initial_al_weight;
        const double scaling;
        const double max_al_weight;
        const double eta_tol;
 
        // TODO: replace this with a member function
        std::function<void(const Eigen::VectorXd &)> update_barrier_stiffness;
    };
} // namespace polyfem::solver