PressureForm.hpp

Go to the documentation of this file.

#pragma once
 
#include "Form.hpp"
 
#include <polyfem/assembler/PressureAssembler.hpp>
#include <polyfem/mesh/LocalBoundary.hpp>
#include <polyfem/utils/ElasticityUtils.hpp>
 
#include <polyfem/utils/Types.hpp>
 
namespace polyfem::solver
{
    class PressureForm : public Form
    {
    public:
        PressureForm(const int ndof,
                     const std::vector<mesh::LocalBoundary> &local_pressure_boundary,
                     const std::unordered_map<int, std::vector<mesh::LocalBoundary>> &local_pressure_cavity,
                     const std::vector<int> &dirichlet_nodes,
                     const QuadratureOrders &n_boundary_samples,
                     const assembler::PressureAssembler &pressure_assembler,
                     const bool is_time_dependent);
 
        std::string name() const override { return "pressure"; }
 
    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:
        void update_quantities(const double t, const Eigen::VectorXd &x) override;
 
        void force_shape_derivative(
            const int n_verts,
            const double t,
            const Eigen::MatrixXd &x,
            const Eigen::MatrixXd &adjoint,
            Eigen::VectorXd &term);
 
        double force_pressure_derivative(
            const int n_verts,
            const double t,
            const int pressure_boundary_id,
            const Eigen::MatrixXd &x,
            const Eigen::MatrixXd &adjoint);
 
    private:
        double t_;       
        const int ndof_; 
 
        const std::vector<mesh::LocalBoundary> &local_pressure_boundary_;
        const std::unordered_map<int, std::vector<mesh::LocalBoundary>> &local_pressure_cavity_;
        const std::vector<int> &dirichlet_nodes_;
        const QuadratureOrders n_boundary_samples_;
 
        const assembler::PressureAssembler &pressure_assembler_; 
    };
} // namespace polyfem::solver