Maps to pressure boundary. More...

#include <PressureBoundaryVariableToSimulation.hpp>

Inheritance diagram for polyfem::solver::PressureBoundaryVariableToSimulation:

Collaboration diagram for polyfem::solver::PressureBoundaryVariableToSimulation:

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

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

Public Member Functions
	PressureBoundaryVariableToSimulation (StatePtrs states, DiffCachePtrs diff_caches, CompositeParametrization parametrizations, Eigen::VectorXi active_boundary_ids, Eigen::VectorXi active_time_slices)
	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

Private Attributes
bool	is_transient_

int	time_steps_

StatePtrs	states_

DiffCachePtrs	diff_caches_

CompositeParametrization	parametrization_

Eigen::VectorXi	active_boundary_ids_

Eigen::VectorXi	active_time_slices_

Detailed Description

Maps to pressure boundary.

Maps one boundary solution to all boundaries in one boundary selection. Basically just the pressure version of dirichet boundary var2sim.

Expect optimiztion variable with layout: [ solution of boundary selection 1 at t0] [solution of boundary selection 2 at t0] ... [ solution of boundary selection 1 at t1] [solution of boundary selection 2 at t1] ...

Definition at line 21 of file PressureBoundaryVariableToSimulation.hpp.

Member Typedef Documentation

◆ DiffCachePtrs

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

Definition at line 25 of file PressureBoundaryVariableToSimulation.hpp.

◆ StatePtrs

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

Definition at line 24 of file PressureBoundaryVariableToSimulation.hpp.

Constructor & Destructor Documentation

◆ PressureBoundaryVariableToSimulation()

polyfem::solver::PressureBoundaryVariableToSimulation::PressureBoundaryVariableToSimulation	(	StatePtrs	states,
		DiffCachePtrs	diff_caches,
		CompositeParametrization	parametrizations,
		Eigen::VectorXi	active_boundary_ids,
		Eigen::VectorXi	active_time_slices
	)

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_boundary_ids	Active pressure boundary ids. Empty implies all active.
[in]	active_time_slices_	Active time slices. Empty implies all active. Ignored for static problem.

Exceptions

std::runtime_error Throw if input is invalid.

Definition at line 23 of file PressureBoundaryVariableToSimulation.cpp.

References active_boundary_ids_, active_time_slices_, diff_caches_, polyfem::solver::is_active_pressure_boundary_ids_valid(), polyfem::solver::is_active_time_slices_valid(), is_transient_, polyfem::utils::json_as_array(), polyfem::log_and_throw_adjoint_error(), states_, and time_steps_.

Here is the call graph for this function:

Member Function Documentation

◆ affect_state()

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

overridevirtual

Return true if current var2sim maps to target state.

Implements polyfem::solver::VariableToSimulation.

Definition at line 100 of file PressureBoundaryVariableToSimulation.cpp.

References states_.

◆ apply_parametrization_jacobian()

Eigen::VectorXd polyfem::solver::PressureBoundaryVariableToSimulation::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 263 of file PressureBoundaryVariableToSimulation.cpp.

References polyfem::log_and_throw_adjoint_error(), and name().

Here is the call graph for this function:

◆ compute_adjoint_term()

Eigen::VectorXd polyfem::solver::PressureBoundaryVariableToSimulation::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.