MatrixLagrangianForm.hpp

Go to the documentation of this file.

#pragma once
 
#include "AugmentedLagrangianForm.hpp"
 
namespace polyfem::solver
{
    class MatrixLagrangianForm : public AugmentedLagrangianForm
    {
    public:
        MatrixLagrangianForm(const StiffnessMatrix &A,
                             const Eigen::MatrixXd &b,
                             const StiffnessMatrix &A_proj = {},
                             const Eigen::MatrixXd &b_pro = {});
 
        std::string name() const override { return "generic-lagrangian"; }
 
        double value_unweighted(const Eigen::VectorXd &x) const override;
 
        void first_derivative_unweighted(const Eigen::VectorXd &x, Eigen::VectorXd &gradv) const override;
 
        void second_derivative_unweighted(const Eigen::VectorXd &x, StiffnessMatrix &hessian) const override;
 
    public:
        void update_lagrangian(const Eigen::VectorXd &x, const double k_al) override;
        double compute_error(const Eigen::VectorXd &x) const override;
 
    private:
        StiffnessMatrix AtA;
        Eigen::VectorXd Atb;
    };
} // namespace polyfem::solver