MacroStrainLagrangianForm.hpp

Go to the documentation of this file.

#pragma once
 
#include "AugmentedLagrangianForm.hpp"
 
namespace polyfem::assembler
{
    class MacroStrainValue;
}
 
namespace polyfem::solver
{
    class MacroStrainLagrangianForm : public AugmentedLagrangianForm
    {
    public:
        MacroStrainLagrangianForm(const assembler::MacroStrainValue &macro_strain_constraint);
 
        std::string name() const override { return "strain-Lagrangian"; }
 
        void update_quantities(const double t, const Eigen::VectorXd &x) override;
        void update_lagrangian(const Eigen::VectorXd &x, const double k_al) override;
        double compute_error(const Eigen::VectorXd &x) const override;
 
    protected:
        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;
 
    private:
        Eigen::VectorXd values;
 
        Eigen::VectorXd lagr_mults_; 
        const assembler::MacroStrainValue &macro_strain_constraint_;
    };
} // namespace polyfem::solver