#include <BarrierForms.hpp>

Inheritance diagram for polyfem::solver::SmoothContactForceForm:

Collaboration diagram for polyfem::solver::SmoothContactForceForm:

Public Member Functions
	SmoothContactForceForm (const VariableToSimulationGroup &variable_to_simulations, const State &state, const json &args)

	~SmoothContactForceForm ()=default

double	value_unweighted_step (const int time_step, const Eigen::VectorXd &x) const override

Eigen::VectorXd	compute_adjoint_rhs_step (const int time_step, const Eigen::VectorXd &x, const State &state) const override

void	compute_partial_gradient_step (const int time_step, const Eigen::VectorXd &x, Eigen::VectorXd &gradv) const override

void	solution_changed_step (const int time_step, const Eigen::VectorXd &x) override

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

virtual std::string	name () const override

Eigen::MatrixXd	compute_adjoint_rhs (const Eigen::VectorXd &x, const State &state) const final override

void	compute_partial_gradient (const Eigen::VectorXd &x, Eigen::VectorXd &gradv) const final override

double	value_unweighted (const Eigen::VectorXd &x) const final override
	Compute the value of the form.

virtual Eigen::VectorXd	compute_adjoint_rhs_step_prev (const int time_step, const Eigen::VectorXd &x, const State &state) const

virtual void	solution_changed (const Eigen::VectorXd &new_x) final override
	Update cached fields upon a change in the solution.

virtual bool	depends_on_step_prev () const final

	AdjointForm (const VariableToSimulationGroup &variable_to_simulations)

Public Member Functions inherited from polyfem::solver::AdjointForm
	AdjointForm (const VariableToSimulationGroup &variable_to_simulations)

virtual	~AdjointForm ()

void	enable_energy_print (const std::string &print_energy_keyword)

double	value (const Eigen::VectorXd &x) const override
	Compute the value of the form multiplied with the weigth.

const VariableToSimulationGroup &	get_variable_to_simulations () const

virtual Eigen::MatrixXd	compute_reduced_adjoint_rhs (const Eigen::VectorXd &x, const State &state) const

virtual void	first_derivative (const Eigen::VectorXd &x, Eigen::VectorXd &gradv) const final override
	Compute the first derivative of the value wrt x multiplied with the weigth.

virtual void	update_quantities (const double t, const Eigen::VectorXd &x) final override
	Update time-dependent fields.

virtual void	init_lagging (const Eigen::VectorXd &x) final override
	Initialize lagged fields TODO: more than one step.

virtual void	update_lagging (const Eigen::VectorXd &x, const int iter_num) final override
	Update lagged fields.

Public Member Functions inherited from polyfem::solver::Form
virtual	~Form ()

virtual void	init (const Eigen::VectorXd &x)
	Initialize the form.

virtual void	finish ()

Eigen::VectorXd	value_per_element (const Eigen::VectorXd &x) const
	Compute the value of the form multiplied with the weigth.

void	second_derivative (const Eigen::VectorXd &x, StiffnessMatrix &hessian) const
	Compute the second derivative of the value wrt x multiplied with the weigth.

virtual bool	is_step_valid (const Eigen::VectorXd &x0, const Eigen::VectorXd &x1) const
	Determine if a step from solution x0 to solution x1 is allowed.

virtual double	max_step_size (const Eigen::VectorXd &x0, const Eigen::VectorXd &x1) const
	Determine the maximum step size allowable between the current and next solution.

virtual void	line_search_begin (const Eigen::VectorXd &x0, const Eigen::VectorXd &x1)
	Initialize variables used during the line search.

virtual void	line_search_end ()
	Clear variables used during the line search.

virtual void	post_step (const polysolve::nonlinear::PostStepData &data)
	Update fields after a step in the optimization.

virtual int	max_lagging_iterations () const
	Get the maximum number of lagging iteration allowable.

virtual bool	uses_lagging () const
	Does this form require lagging?

void	set_project_to_psd (bool val)
	Set project to psd.

bool	is_project_to_psd () const
	Get if the form's second derivative is projected to psd.

void	enable ()
	Enable the form.

void	disable ()
	Disable the form.

void	set_enabled (const bool enabled)
	Set if the form is enabled.

bool	enabled () const
	Determine if the form is enabled.

virtual double	weight () const
	Get the form's multiplicative constant weight.

void	set_weight (const double weight)
	Set the form's multiplicative constant weight.

virtual bool	is_step_collision_free (const Eigen::VectorXd &x0, const Eigen::VectorXd &x1) const
	Checks if the step is collision free.

void	set_output_dir (const std::string &output_dir)

virtual void	set_scale (const double scale)
	sets the scale for the form

Protected Member Functions
void	build_collision_mesh ()

ipc::SmoothCollisions	get_smooth_collision_set (const Eigen::MatrixXd &displaced_surface)

Protected Member Functions inherited from polyfem::solver::AdjointForm
virtual void	first_derivative_unweighted (const Eigen::VectorXd &x, Eigen::VectorXd &gradv) const final override
	Compute the first derivative of the value wrt x.

virtual void	second_derivative_unweighted (const Eigen::VectorXd &x, StiffnessMatrix &hessian) const final override
	Compute the second derivative of the value wrt x.

Protected Member Functions inherited from polyfem::solver::Form
std::string	resolve_output_path (const std::string &path) const

virtual Eigen::VectorXd	value_per_element_unweighted (const Eigen::VectorXd &x) const
	Compute the value of the form multiplied per element.

Protected Attributes
const State &	state_

std::set< int >	boundary_ids_

std::map< int, std::set< int > >	boundary_ids_to_dof_

ipc::CollisionMesh	collision_mesh_

ipc::SmoothCollisions	collisions_

const ipc::ParameterType	params_

const double	dmin_ = 0

ipc::SmoothContactPotential	potential_

Protected Attributes inherited from polyfem::solver::StaticForm
bool	depends_on_step_prev_ = false

Protected Attributes inherited from polyfem::solver::AdjointForm
const VariableToSimulationGroup	variable_to_simulations_

PrintStage	print_energy_ = PrintStage::Inactive

std::string	print_energy_keyword_

Protected Attributes inherited from polyfem::solver::Form
bool	project_to_psd_ = false
	If true, the form's second derivative is projected to be positive semidefinite.

double	weight_ = 1
	weight of the form (e.g., AL penalty weight or Δt²)

bool	enabled_ = true
	If true, the form is enabled.

std::string	output_dir_

Additional Inherited Members
Protected Types inherited from polyfem::solver::AdjointForm
enum class	PrintStage { Inactive , AlreadyPrinted , ToPrint }

Detailed Description

Definition at line 105 of file BarrierForms.hpp.

Constructor & Destructor Documentation

◆ SmoothContactForceForm()

polyfem::solver::SmoothContactForceForm::SmoothContactForceForm	(	const VariableToSimulationGroup &	variable_to_simulations,
		const State &	state,
		const json &	args
	)

Definition at line 375 of file BarrierForms.cpp.

References boundary_ids_, and build_collision_mesh().

Here is the call graph for this function:

◆ ~SmoothContactForceForm()

polyfem::solver::SmoothContactForceForm::~SmoothContactForceForm ( )

default

Member Function Documentation

◆ build_collision_mesh()

void polyfem::solver::SmoothContactForceForm::build_collision_mesh ( )

protected

Definition at line 390 of file BarrierForms.cpp.

References polyfem::State::bases, polyfem::State::collision_mesh, collision_mesh_, polyfem::State::mesh, polyfem::State::n_bases, and state_.

Referenced by SmoothContactForceForm(), and solution_changed_step().

Here is the caller graph for this function:

◆ compute_adjoint_rhs_step()

Eigen::VectorXd polyfem::solver::SmoothContactForceForm::compute_adjoint_rhs_step	(	const int	time_step,
		const Eigen::VectorXd &	x,
		const State &	state
	)		const

overridevirtual

Implements polyfem::solver::StaticForm.

Definition at line 444 of file BarrierForms.cpp.

References collision_mesh_, collisions_, polyfem::solver::SolveData::contact_form, polyfem::State::diff_cached, potential_, polyfem::State::solve_data, state_, polyfem::solver::DiffCache::u(), polyfem::utils::unflatten(), and polyfem::solver::Form::weight().

Here is the call graph for this function:

◆ compute_partial_gradient_step()

void polyfem::solver::SmoothContactForceForm::compute_partial_gradient_step	(	const int	time_step,
		const Eigen::VectorXd &	x,
		Eigen::VectorXd &	gradv
	)		const

overridevirtual

Implements polyfem::solver::StaticForm.

Definition at line 462 of file BarrierForms.cpp.

References polyfem::solver::VariableToSimulationGroup::apply_parametrization_jacobian(), polyfem::State::basis_nodes_to_gbasis_nodes, collision_mesh_, collisions_, polyfem::solver::SolveData::contact_form, polyfem::State::diff_cached, polyfem::solver::AdjointTools::map_node_to_primitive_order(), potential_, polyfem::solver::Shape, polyfem::State::solve_data, state_, polyfem::solver::DiffCache::u(), polyfem::utils::unflatten(), polyfem::solver::AdjointForm::variable_to_simulations_, polyfem::solver::Form::weight(), and x.

Here is the call graph for this function:

◆ get_smooth_collision_set()

ipc::SmoothCollisions polyfem::solver::SmoothContactForceForm::get_smooth_collision_set ( const Eigen::MatrixXd & displaced_surface )

protected

Definition at line 420 of file BarrierForms.cpp.

References collision_mesh_, polyfem::solver::SolveData::contact_form, polyfem::State::solve_data, and state_.

Referenced by solution_changed_step().

Here is the caller graph for this function:

◆ solution_changed_step()

void polyfem::solver::SmoothContactForceForm::solution_changed_step	(	const int	time_step,
		const Eigen::VectorXd &	x
	)

overridevirtual

Reimplemented from polyfem::solver::StaticForm.

Definition at line 487 of file BarrierForms.cpp.

References build_collision_mesh(), collision_mesh_, collisions_, polyfem::State::diff_cached, get_smooth_collision_set(), state_, polyfem::solver::DiffCache::u(), and polyfem::utils::unflatten().

Here is the call graph for this function:

◆ value_unweighted_step()

double polyfem::solver::SmoothContactForceForm::value_unweighted_step	(	const int	time_step,
		const Eigen::VectorXd &	x
	)		const

overridevirtual

Implements polyfem::solver::StaticForm.

Definition at line 428 of file BarrierForms.cpp.

References collision_mesh_, collisions_, polyfem::solver::SolveData::contact_form, polyfem::State::diff_cached, potential_, polyfem::State::solve_data, state_, polyfem::solver::DiffCache::u(), and polyfem::utils::unflatten().

Here is the call graph for this function:

Member Data Documentation

◆ boundary_ids_

std::set<int> polyfem::solver::SmoothContactForceForm::boundary_ids_

protected

Definition at line 124 of file BarrierForms.hpp.

Referenced by SmoothContactForceForm().

◆ boundary_ids_to_dof_

std::map<int, std::set<int> > polyfem::solver::SmoothContactForceForm::boundary_ids_to_dof_

protected

Definition at line 125 of file BarrierForms.hpp.

◆ collision_mesh_

ipc::CollisionMesh polyfem::solver::SmoothContactForceForm::collision_mesh_

protected

Definition at line 127 of file BarrierForms.hpp.

Referenced by build_collision_mesh(), compute_adjoint_rhs_step(), compute_partial_gradient_step(), get_smooth_collision_set(), solution_changed_step(), and value_unweighted_step().

◆ collisions_

ipc::SmoothCollisions polyfem::solver::SmoothContactForceForm::collisions_

protected

Definition at line 128 of file BarrierForms.hpp.

Referenced by compute_adjoint_rhs_step(), compute_partial_gradient_step(), solution_changed_step(), and value_unweighted_step().

◆ dmin_

const double polyfem::solver::SmoothContactForceForm::dmin_ = 0

protected

Definition at line 130 of file BarrierForms.hpp.

◆ params_

const ipc::ParameterType polyfem::solver::SmoothContactForceForm::params_

protected

Definition at line 129 of file BarrierForms.hpp.

◆ potential_

ipc::SmoothContactPotential polyfem::solver::SmoothContactForceForm::potential_

protected

Definition at line 132 of file BarrierForms.hpp.

Referenced by compute_adjoint_rhs_step(), compute_partial_gradient_step(), and value_unweighted_step().

◆ state_

const State& polyfem::solver::SmoothContactForceForm::state_

protected

Definition at line 123 of file BarrierForms.hpp.

Referenced by build_collision_mesh(), compute_adjoint_rhs_step(), compute_partial_gradient_step(), get_smooth_collision_set(), solution_changed_step(), and value_unweighted_step().

The documentation for this class was generated from the following files:

src/polyfem/solver/forms/adjoint_forms/BarrierForms.hpp
src/polyfem/solver/forms/adjoint_forms/BarrierForms.cpp

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ SmoothContactForceForm()

◆ ~SmoothContactForceForm()

Member Function Documentation

◆ build_collision_mesh()

◆ compute_adjoint_rhs_step()

◆ compute_partial_gradient_step()

◆ get_smooth_collision_set()

◆ solution_changed_step()

◆ value_unweighted_step()

Member Data Documentation

◆ boundary_ids_

◆ boundary_ids_to_dof_

◆ collision_mesh_

◆ collisions_

◆ dmin_

◆ params_

◆ potential_

◆ state_