DirichletBoundaryVariableToSimulation.hpp

Go to the documentation of this file.

#pragma once
 
#include <polyfem/optimization/var2sims/VariableToSimulation.hpp>
 
#include <Eigen/Core>
 
#include <memory>
#include <string>
#include <vector>
 
namespace polyfem::solver
{
    class DirichletBoundaryVariableToSimulation : public VariableToSimulation
    {
    public:
        using VarFormPtrs = std::vector<std::shared_ptr<varform::DifferentiableVarForm>>;
        using DiffCachePtrs = std::vector<std::shared_ptr<DiffCache>>;
 
        DirichletBoundaryVariableToSimulation(
            VarFormPtrs varforms,
            DiffCachePtrs diff_caches,
            CompositeParametrization parametrizations,
            Eigen::VectorXi active_boundary_ids,
            Eigen::VectorXi active_time_slices);
 
        std::string name() const override;
        ParameterType parameter_type() const override;
        bool affects_varform(const varform::DifferentiableVarForm &target) const override;
        void update(const Eigen::VectorXd &x) override;
        void update_state_variables(const Eigen::VectorXd &x, Eigen::VectorXd &state_variables) const override;
        Eigen::VectorXd compute_adjoint_term(const Eigen::VectorXd &x) const override;
        int inverse_dof() const override;
        Eigen::VectorXd inverse_eval() const override;
        Eigen::VectorXd apply_parametrization_jacobian(const Eigen::VectorXd &term, const Eigen::VectorXd &x) const override;
 
    private:
        using BoundaryNodeMap = std::vector<std::vector<std::vector<int>>>;
 
        int dim_;
        int time_steps_;
 
        VarFormPtrs varforms_;
        DiffCachePtrs diff_caches_;
        CompositeParametrization parametrization_;
 
        Eigen::VectorXi active_boundary_ids_;
        Eigen::VectorXi active_time_slices_;
 
        std::vector<BoundaryNodeMap> boundary_node_maps_;
 
        int para_out_dof() const;
        void build_boundary_node_maps();
    };
 
} // namespace polyfem::solver