PolyFEM
|
#include <SpatialIntegralForms.hpp>
Public Member Functions | |
SpatialIntegralForm (const VariableToSimulationGroup &variable_to_simulations, const State &state, const json &args) | |
std::string | name () const override |
const State & | get_state () |
void | set_integral_type (const SpatialIntegralType type) |
Eigen::VectorXd | compute_adjoint_rhs_step (const int time_step, const Eigen::VectorXd &x, const State &state) const override |
double | value_unweighted_step (const int time_step, const Eigen::VectorXd &x) const override |
void | compute_partial_gradient_step (const int time_step, const Eigen::VectorXd &x, Eigen::VectorXd &gradv) const override |
Public Member Functions inherited from polyfem::solver::StaticForm | |
virtual | ~StaticForm ()=default |
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 void | solution_changed_step (const int time_step, const Eigen::VectorXd &new_x) |
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. | |
virtual void | set_apply_DBC (const Eigen::VectorXd &x, bool apply_DBC) final override |
Set if the Dirichlet boundary conditions should be enforced. | |
Public Member Functions inherited from polyfem::solver::Form | |
virtual | ~Form () |
virtual void | init (const Eigen::VectorXd &x) |
Initialize the form. | |
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) |
Protected Member Functions | |
virtual IntegrableFunctional | get_integral_functional () const =0 |
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_ |
SpatialIntegralType | spatial_integral_type_ |
std::set< int > | ids_ |
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 } |
Definition at line 7 of file SpatialIntegralForms.hpp.
|
inline |
Definition at line 10 of file SpatialIntegralForms.hpp.
|
overridevirtual |
Implements polyfem::solver::StaticForm.
Definition at line 82 of file SpatialIntegralForms.cpp.
References polyfem::State::diff_cached, polyfem::solver::AdjointTools::dJ_du_step(), get_integral_functional(), ids_, polyfem::State::ndof(), polyfem::solver::DiffCache::size(), spatial_integral_type_, state_, polyfem::solver::DiffCache::u(), and polyfem::solver::Form::weight().
Referenced by compute_partial_gradient_step().
|
overridevirtual |
Implements polyfem::solver::StaticForm.
Reimplemented in polyfem::solver::StressNormForm, and polyfem::solver::StressForm.
Definition at line 60 of file SpatialIntegralForms.cpp.
References polyfem::solver::VariableToSimulationGroup::apply_parametrization_jacobian(), compute_adjoint_rhs_step(), polyfem::solver::AdjointTools::compute_shape_derivative_functional_term(), polyfem::State::diff_cached, polyfem::solver::DiffCache::disp_grad(), polyfem::utils::flatten(), polyfem::State::mesh, polyfem::solver::SolveData::nl_problem, polyfem::solver::PeriodicShape, polyfem::State::primitive_to_node(), polyfem::solver::NLHomoProblem::reduced_to_full_shape_derivative(), polyfem::solver::Shape, polyfem::solver::DiffCache::size(), polyfem::State::solve_data, state_, polyfem::solver::DiffCache::u(), polyfem::utils::unflatten(), polyfem::solver::AdjointForm::variable_to_simulations_, polyfem::solver::Form::weight(), and x.
Referenced by polyfem::solver::ElasticEnergyForm::compute_partial_gradient_step(), polyfem::solver::StressNormForm::compute_partial_gradient_step(), polyfem::solver::ComplianceForm::compute_partial_gradient_step(), and polyfem::solver::StressForm::compute_partial_gradient_step().
|
protectedpure virtual |
Implemented in polyfem::solver::ElasticEnergyForm, polyfem::solver::StressNormForm, polyfem::solver::ComplianceForm, polyfem::solver::PositionForm, polyfem::solver::AccelerationForm, polyfem::solver::KineticForm, polyfem::solver::StressForm, polyfem::solver::VolumeForm, polyfem::solver::TargetForm, polyfem::solver::SDFTargetForm, and polyfem::solver::MeshTargetForm.
Referenced by compute_adjoint_rhs_step(), and value_unweighted_step().
|
inline |
Definition at line 16 of file SpatialIntegralForms.hpp.
References state_.
|
inlineoverridevirtual |
Reimplemented from polyfem::solver::StaticForm.
Reimplemented in polyfem::solver::StressNormForm, polyfem::solver::StressForm, and polyfem::solver::TargetForm.
Definition at line 14 of file SpatialIntegralForms.hpp.
Referenced by polyfem::solver::AccelerationForm::get_integral_functional(), and polyfem::solver::KineticForm::get_integral_functional().
|
inline |
Definition at line 18 of file SpatialIntegralForms.hpp.
References spatial_integral_type_.
Referenced by polyfem::solver::AccelerationForm::AccelerationForm(), polyfem::solver::ComplianceForm::ComplianceForm(), polyfem::solver::ElasticEnergyForm::ElasticEnergyForm(), polyfem::solver::KineticForm::KineticForm(), polyfem::solver::MeshTargetForm::MeshTargetForm(), polyfem::solver::PositionForm::PositionForm(), polyfem::solver::SDFTargetForm::SDFTargetForm(), polyfem::solver::StressForm::StressForm(), polyfem::solver::StressNormForm::StressNormForm(), polyfem::solver::TargetForm::TargetForm(), and polyfem::solver::VolumeForm::VolumeForm().
|
overridevirtual |
Implements polyfem::solver::StaticForm.
Definition at line 54 of file SpatialIntegralForms.cpp.
References polyfem::State::diff_cached, get_integral_functional(), ids_, polyfem::solver::AdjointTools::integrate_objective(), polyfem::solver::DiffCache::size(), spatial_integral_type_, state_, and polyfem::solver::DiffCache::u().
|
protected |
Definition at line 29 of file SpatialIntegralForms.hpp.
Referenced by polyfem::solver::AccelerationForm::AccelerationForm(), polyfem::solver::ComplianceForm::ComplianceForm(), compute_adjoint_rhs_step(), polyfem::solver::ElasticEnergyForm::ElasticEnergyForm(), polyfem::solver::KineticForm::KineticForm(), polyfem::solver::MeshTargetForm::MeshTargetForm(), polyfem::solver::PositionForm::PositionForm(), polyfem::solver::SDFTargetForm::SDFTargetForm(), polyfem::solver::StressForm::StressForm(), polyfem::solver::StressNormForm::StressNormForm(), polyfem::solver::TargetForm::TargetForm(), value_unweighted_step(), and polyfem::solver::VolumeForm::VolumeForm().
|
protected |
Definition at line 28 of file SpatialIntegralForms.hpp.
Referenced by compute_adjoint_rhs_step(), set_integral_type(), and value_unweighted_step().
|
protected |
Definition at line 27 of file SpatialIntegralForms.hpp.
Referenced by compute_adjoint_rhs_step(), compute_partial_gradient_step(), polyfem::solver::ElasticEnergyForm::compute_partial_gradient_step(), polyfem::solver::StressNormForm::compute_partial_gradient_step(), polyfem::solver::ComplianceForm::compute_partial_gradient_step(), polyfem::solver::StressForm::compute_partial_gradient_step(), polyfem::solver::ElasticEnergyForm::get_integral_functional(), polyfem::solver::StressNormForm::get_integral_functional(), polyfem::solver::ComplianceForm::get_integral_functional(), polyfem::solver::AccelerationForm::get_integral_functional(), polyfem::solver::KineticForm::get_integral_functional(), polyfem::solver::StressForm::get_integral_functional(), polyfem::solver::TargetForm::get_integral_functional(), get_state(), polyfem::solver::SDFTargetForm::set_bspline_target(), polyfem::solver::SDFTargetForm::set_bspline_target(), polyfem::solver::TargetForm::set_reference(), polyfem::solver::MeshTargetForm::set_surface_mesh_target(), polyfem::solver::SDFTargetForm::solution_changed_step(), polyfem::solver::MeshTargetForm::solution_changed_step(), and value_unweighted_step().