PeriodicLagrangianForm.hpp

Go to the documentation of this file.

#pragma once
 
#include "AugmentedLagrangianForm.hpp"
#include <polyfem/assembler/PeriodicBoundary.hpp>
 
namespace polyfem::solver
{
    class PeriodicLagrangianForm : public AugmentedLagrangianForm
    {
    public:
        PeriodicLagrangianForm(const int ndof,
                               const std::shared_ptr<utils::PeriodicBoundary> &periodic_bc);
 
        std::string name() const override
        {
            return "periodic-alagrangian";
        }
 
        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_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;
 
        virtual bool can_project() const override;
        virtual void project_gradient(Eigen::VectorXd &grad) const override;
        virtual void project_diag(Eigen::VectorXd &diag) const override;
        virtual void project_hessian(StiffnessMatrix &hessian) const override;
 
    private:
        const std::shared_ptr<utils::PeriodicBoundary> periodic_bc_;
        const int n_dofs_;
        Eigen::VectorXi constraints_;     
        Eigen::VectorXi not_constraints_; 
    };
} // namespace polyfem::solver