Maps to mesh vertex position. More...

#include <ShapeVariableToSimulation.hpp>

Inheritance diagram for polyfem::solver::ShapeVariableToSimulation:

Collaboration diagram for polyfem::solver::ShapeVariableToSimulation:

Public Types
using	StatePtrs = std::vector< std::shared_ptr< legacy::State > >

using	DiffCachePtrs = std::vector< std::shared_ptr< DiffCache > >

Public Member Functions
	ShapeVariableToSimulation (StatePtrs states, DiffCachePtrs diff_caches, CompositeParametrization parametrizations, Eigen::VectorXi active_dimensions, Eigen::VectorXi active_geom_nodes)
	Construct ShapeVariableToSimulation.

std::string	name () const override

ParameterType	parameter_type () const override

bool	affect_state (const legacy::State &target) const override
	Return true if current var2sim maps to target state.

void	update (const Eigen::VectorXd &x) override
	Update forward simulation states from optimization variables.

void	update_state_variables (const Eigen::VectorXd &x, Eigen::VectorXd &state_variables) const override
	Update state variables from optimization variables.

Eigen::VectorXd	compute_adjoint_term (const Eigen::VectorXd &x) const override
	Compute adjoint contribution of objective gradient.

int	inverse_dof () const override
	Compute optimization variables dof.

Eigen::VectorXd	inverse_eval () const override
	Compute optimization variables from forward simulation legacy::State.

Eigen::VectorXd	apply_parametrization_jacobian (const Eigen::VectorXd &term, const Eigen::VectorXd &x) const override
	Apply parametrization jacobian to compute the gradient w.r.t.

Public Member Functions inherited from polyfem::solver::VariableToSimulation
virtual	~VariableToSimulation ()=default

Private Member Functions
int	para_out_dof () const
	Return variable dof after parametrization mapping.

Private Attributes
int	dim_

int	vertex_num_

StatePtrs	states_

DiffCachePtrs	diff_caches_

CompositeParametrization	parametrization_

Eigen::VectorXi	active_dimensions_

Eigen::VectorXi	active_geom_nodes_

Detailed Description

Maps to mesh vertex position.

Expect optimization variable layout: [vertex 0 position] [vertex 1 position] ...

Definition at line 15 of file ShapeVariableToSimulation.hpp.

Member Typedef Documentation

◆ DiffCachePtrs

using polyfem::solver::ShapeVariableToSimulation::DiffCachePtrs = std::vector<std::shared_ptr<DiffCache> >

Definition at line 19 of file ShapeVariableToSimulation.hpp.

◆ StatePtrs

using polyfem::solver::ShapeVariableToSimulation::StatePtrs = std::vector<std::shared_ptr<legacy::State> >

Definition at line 18 of file ShapeVariableToSimulation.hpp.

Constructor & Destructor Documentation

◆ ShapeVariableToSimulation()

polyfem::solver::ShapeVariableToSimulation::ShapeVariableToSimulation	(	StatePtrs	states,
		DiffCachePtrs	diff_caches,
		CompositeParametrization	parametrizations,
		Eigen::VectorXi	active_dimensions,
		Eigen::VectorXi	active_geom_nodes
	)

Construct ShapeVariableToSimulation.

Parameters

[in]	states	Shared ptr to all forward sim states.
[in]	diff_caches	Shared ptr to all diff caches.
[in]	parametrizations	Parametrizations.
[in]	active_dimensions	Vector of active dimensions per mesh node. Empty implies all active.
[in]	active_geom_nodes	Vector of active nodes. Empty implies all active.

Exceptions

std::runtime_error Throw if input is invalid.

Definition at line 17 of file ShapeVariableToSimulation.cpp.

References active_dimensions_, active_geom_nodes_, diff_caches_, dim_, polyfem::solver::is_active_dims_valid(), polyfem::solver::is_active_geom_nodes_valid(), polyfem::log_and_throw_adjoint_error(), states_, and vertex_num_.

Here is the call graph for this function:

Member Function Documentation

◆ affect_state()

bool polyfem::solver::ShapeVariableToSimulation::affect_state ( const legacy::State & target ) const

overridevirtual

Return true if current var2sim maps to target state.

Implements polyfem::solver::VariableToSimulation.

Definition at line 65 of file ShapeVariableToSimulation.cpp.

References states_.

◆ apply_parametrization_jacobian()

Eigen::VectorXd polyfem::solver::ShapeVariableToSimulation::apply_parametrization_jacobian	(	const Eigen::VectorXd &	term,
		const Eigen::VectorXd &	x
	)		const

overridevirtual

Apply parametrization jacobian to compute the gradient w.r.t.

to optimization variables.

Parameters

term	Gradient w.r.t. to full inherent dof (Ex. all vertices for shape var2sim)
x	Optimization variables.

Implements polyfem::solver::VariableToSimulation.

Definition at line 201 of file ShapeVariableToSimulation.cpp.

References active_dimensions_, active_geom_nodes_, polyfem::solver::CompositeParametrization::apply_jacobian(), dim_, para_out_dof(), parametrization_, vertex_num_, and x.

Here is the call graph for this function:

◆ compute_adjoint_term()

Eigen::VectorXd polyfem::solver::ShapeVariableToSimulation::compute_adjoint_term ( const Eigen::VectorXd & x ) const

overridevirtual

Compute adjoint contribution of objective gradient.

See arXiv:2205.13643.

Let objective function be J and optimization variable be x. This method computes all adjoint related terms in dJ/dx. For static case (E.q. 11), it's the second term. For dynamic case (E.q. 15), it's everything except ∂qJ.

Parameters

x[in] Optimization variables.