InversionBarrierForm.hpp

Go to the documentation of this file.

 
#pragma once
 
#include "Form.hpp"
 
#include <polyfem/utils/Types.hpp>
 
namespace polyfem::solver
{
    class InversionBarrierForm : public polyfem::solver::Form
    {
    public:
        InversionBarrierForm(
            const Eigen::MatrixXd &rest_positions, const Eigen::MatrixXi &elements, const int dim, const double vhat);
 
        std::string name() const override { return "inversion_barrier"; }
 
    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;
 
    public:
        bool is_step_valid(const Eigen::VectorXd &x0, const Eigen::VectorXd &x1) const override;
 
    protected:
        static double element_volume(const Eigen::MatrixXd &element_vertices);
        static Eigen::VectorXd element_volume_gradient(const Eigen::MatrixXd &element_vertices);
        static Eigen::MatrixXd element_volume_hessian(const Eigen::MatrixXd &element_vertices);
 
    private:
        Eigen::MatrixXd rest_positions_;
        Eigen::MatrixXi elements_;
        int dim_;
        double vhat_;
    };
} // namespace polyfem::solver