GenericLagrangianForm.hpp

Go to the documentation of this file.

#pragma once
 
#include "AugmentedLagrangianForm.hpp"
 
namespace polyfem::solver
{
    class GenericLagrangianForm : public AugmentedLagrangianForm
    {
    public:
        GenericLagrangianForm(const std::vector<int> &constraint_nodes,
                              const StiffnessMatrix &A,
                              const Eigen::VectorXd &b);
 
        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:
        const StiffnessMatrix A;
        const Eigen::VectorXd b;
    };
} // namespace polyfem::solver