L2ProjectionForm.hpp

Go to the documentation of this file.

 
#pragma once
 
#include "Form.hpp"
 
#include <polyfem/utils/Types.hpp>
 
namespace polyfem::solver
{
 
    class L2ProjectionForm : public Form
    {
    public:
        L2ProjectionForm(
            const StiffnessMatrix &M,
            const StiffnessMatrix &A,
            const Eigen::VectorXd &y);
 
        std::string name() const override { return "L2_projection"; }
 
    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;
 
        StiffnessMatrix M_;
        Eigen::VectorXd rhs_;
    };
 
} // namespace polyfem::solver