PolyFEM
|
Form of the augmented lagrangian. More...
#include <AugmentedLagrangianForm.hpp>
Public Member Functions | |
AugmentedLagrangianForm () | |
virtual | ~AugmentedLagrangianForm () |
virtual void | update_lagrangian (const Eigen::VectorXd &x, const double k_al)=0 |
virtual double | compute_error (const Eigen::VectorXd &x) const =0 |
void | set_initial_weight (const double k_al) |
double | lagrangian_weight () const |
const StiffnessMatrix & | constraint_matrix () const |
const Eigen::MatrixXd & | constraint_value () const |
const StiffnessMatrix & | constraint_projection_matrix () const |
const Eigen::MatrixXd & | constraint_projection_vector () const |
bool | has_projection () const |
virtual bool | can_project () const |
virtual void | project_gradient (Eigen::VectorXd &grad) const |
virtual void | project_hessian (StiffnessMatrix &hessian) const |
void | set_scale (const double scale) override |
sets the scale for the form | |
![]() | |
virtual | ~Form () |
virtual std::string | name () const =0 |
virtual void | init (const Eigen::VectorXd &x) |
Initialize the form. | |
virtual void | finish () |
virtual double | value (const Eigen::VectorXd &x) const |
Compute the value of the form multiplied with the weigth. | |
Eigen::VectorXd | value_per_element (const Eigen::VectorXd &x) const |
Compute the value of the form multiplied with the weigth. | |
virtual void | first_derivative (const Eigen::VectorXd &x, Eigen::VectorXd &gradv) const |
Compute the first derivative of the value wrt x 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 void | solution_changed (const Eigen::VectorXd &new_x) |
Update cached fields upon a change in the solution. | |
virtual void | update_quantities (const double t, const Eigen::VectorXd &x) |
Update time-dependent fields. | |
virtual void | init_lagging (const Eigen::VectorXd &x) |
Initialize lagged fields TODO: more than one step. | |
virtual void | update_lagging (const Eigen::VectorXd &x, const int iter_num) |
Update lagged fields. | |
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 | |
double | L_weight () const |
double | A_weight () const |
![]() | |
std::string | resolve_output_path (const std::string &path) const |
virtual double | value_unweighted (const Eigen::VectorXd &x) const =0 |
Compute the value of the form. | |
virtual Eigen::VectorXd | value_per_element_unweighted (const Eigen::VectorXd &x) const |
Compute the value of the form multiplied per element. | |
virtual void | first_derivative_unweighted (const Eigen::VectorXd &x, Eigen::VectorXd &gradv) const =0 |
Compute the first derivative of the value wrt x. | |
virtual void | second_derivative_unweighted (const Eigen::VectorXd &x, StiffnessMatrix &hessian) const =0 |
Compute the second derivative of the value wrt x. | |
Protected Attributes | |
double | k_al_ |
penalty parameter | |
Eigen::VectorXd | lagr_mults_ |
vector of lagrange multipliers | |
StiffnessMatrix | A_ |
Constraints matrix. | |
Eigen::MatrixXd | b_ |
Constraints value. | |
StiffnessMatrix | A_proj_ |
Constraints projection matrix. | |
Eigen::MatrixXd | b_proj_ |
Constraints projection value. | |
![]() | |
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_ |
Private Attributes | |
double | k_scale_ = 1 |
Form of the augmented lagrangian.
Definition at line 8 of file AugmentedLagrangianForm.hpp.
|
inline |
Definition at line 12 of file AugmentedLagrangianForm.hpp.
|
inlinevirtual |
Definition at line 14 of file AugmentedLagrangianForm.hpp.
|
inlineprotected |
Definition at line 42 of file AugmentedLagrangianForm.hpp.
References k_al_, and k_scale_.
Referenced by polyfem::solver::BCLagrangianForm::first_derivative_unweighted(), polyfem::solver::MacroStrainLagrangianForm::first_derivative_unweighted(), polyfem::solver::MatrixLagrangianForm::first_derivative_unweighted(), polyfem::solver::PeriodicLagrangianForm::first_derivative_unweighted(), polyfem::solver::BCLagrangianForm::second_derivative_unweighted(), polyfem::solver::MacroStrainLagrangianForm::second_derivative_unweighted(), polyfem::solver::MatrixLagrangianForm::second_derivative_unweighted(), polyfem::solver::PeriodicLagrangianForm::second_derivative_unweighted(), polyfem::solver::BCLagrangianForm::value_unweighted(), polyfem::solver::MacroStrainLagrangianForm::value_unweighted(), polyfem::solver::MatrixLagrangianForm::value_unweighted(), and polyfem::solver::PeriodicLagrangianForm::value_unweighted().
|
inlinevirtual |
Reimplemented in polyfem::solver::BCLagrangianForm, and polyfem::solver::PeriodicLagrangianForm.
Definition at line 32 of file AugmentedLagrangianForm.hpp.
|
pure virtual |
|
inline |
Definition at line 24 of file AugmentedLagrangianForm.hpp.
References A_.
|
inline |
Definition at line 27 of file AugmentedLagrangianForm.hpp.
References A_proj_.
|
inline |
Definition at line 28 of file AugmentedLagrangianForm.hpp.
References b_proj_.
|
inline |
Definition at line 25 of file AugmentedLagrangianForm.hpp.
References b_.
|
inline |
Definition at line 30 of file AugmentedLagrangianForm.hpp.
References A_proj_.
|
inlineprotected |
Definition at line 41 of file AugmentedLagrangianForm.hpp.
References k_scale_.
Referenced by polyfem::solver::BCLagrangianForm::first_derivative_unweighted(), polyfem::solver::MacroStrainLagrangianForm::first_derivative_unweighted(), polyfem::solver::MatrixLagrangianForm::first_derivative_unweighted(), polyfem::solver::PeriodicLagrangianForm::first_derivative_unweighted(), polyfem::solver::BCLagrangianForm::value_unweighted(), polyfem::solver::MacroStrainLagrangianForm::value_unweighted(), polyfem::solver::MatrixLagrangianForm::value_unweighted(), and polyfem::solver::PeriodicLagrangianForm::value_unweighted().
|
inline |
Definition at line 22 of file AugmentedLagrangianForm.hpp.
References k_al_.
|
inlinevirtual |
Reimplemented in polyfem::solver::BCLagrangianForm, and polyfem::solver::PeriodicLagrangianForm.
Definition at line 33 of file AugmentedLagrangianForm.hpp.
|
inlinevirtual |
Reimplemented in polyfem::solver::BCLagrangianForm, and polyfem::solver::PeriodicLagrangianForm.
Definition at line 34 of file AugmentedLagrangianForm.hpp.
|
inline |
Definition at line 20 of file AugmentedLagrangianForm.hpp.
References k_al_.
|
inlineoverridevirtual |
sets the scale for the form
scale |
Reimplemented from polyfem::solver::Form.
Definition at line 38 of file AugmentedLagrangianForm.hpp.
References k_scale_.
|
pure virtual |
|
protected |
Constraints matrix.
Definition at line 48 of file AugmentedLagrangianForm.hpp.
Referenced by polyfem::solver::BCLagrangianForm::compute_error(), polyfem::solver::MatrixLagrangianForm::compute_error(), polyfem::solver::PeriodicLagrangianForm::compute_error(), constraint_matrix(), polyfem::solver::BCLagrangianForm::first_derivative_unweighted(), polyfem::solver::MatrixLagrangianForm::first_derivative_unweighted(), polyfem::solver::PeriodicLagrangianForm::first_derivative_unweighted(), polyfem::solver::BCLagrangianForm::init_masked_lumped_mass(), polyfem::solver::MatrixLagrangianForm::MatrixLagrangianForm(), polyfem::solver::PeriodicLagrangianForm::PeriodicLagrangianForm(), polyfem::solver::BCLagrangianForm::second_derivative_unweighted(), polyfem::solver::PeriodicLagrangianForm::second_derivative_unweighted(), polyfem::solver::BCLagrangianForm::update_lagrangian(), polyfem::solver::MatrixLagrangianForm::update_lagrangian(), polyfem::solver::PeriodicLagrangianForm::update_lagrangian(), polyfem::solver::BCLagrangianForm::value_unweighted(), polyfem::solver::MatrixLagrangianForm::value_unweighted(), and polyfem::solver::PeriodicLagrangianForm::value_unweighted().
|
protected |
Constraints projection matrix.
Definition at line 51 of file AugmentedLagrangianForm.hpp.
Referenced by constraint_projection_matrix(), has_projection(), polyfem::solver::BCLagrangianForm::init_masked_lumped_mass(), polyfem::solver::MatrixLagrangianForm::MatrixLagrangianForm(), and polyfem::solver::PeriodicLagrangianForm::PeriodicLagrangianForm().
|
protected |
Constraints value.
Definition at line 49 of file AugmentedLagrangianForm.hpp.
Referenced by polyfem::solver::BCLagrangianForm::BCLagrangianForm(), polyfem::solver::BCLagrangianForm::compute_error(), polyfem::solver::MatrixLagrangianForm::compute_error(), polyfem::solver::PeriodicLagrangianForm::compute_error(), constraint_value(), polyfem::solver::BCLagrangianForm::first_derivative_unweighted(), polyfem::solver::PeriodicLagrangianForm::first_derivative_unweighted(), polyfem::solver::MatrixLagrangianForm::MatrixLagrangianForm(), polyfem::solver::PeriodicLagrangianForm::PeriodicLagrangianForm(), polyfem::solver::BCLagrangianForm::update_lagrangian(), polyfem::solver::MatrixLagrangianForm::update_lagrangian(), polyfem::solver::PeriodicLagrangianForm::update_lagrangian(), polyfem::solver::BCLagrangianForm::update_target(), polyfem::solver::BCLagrangianForm::value_unweighted(), polyfem::solver::MatrixLagrangianForm::value_unweighted(), and polyfem::solver::PeriodicLagrangianForm::value_unweighted().
|
protected |
Constraints projection value.
Definition at line 52 of file AugmentedLagrangianForm.hpp.
Referenced by polyfem::solver::BCLagrangianForm::BCLagrangianForm(), constraint_projection_vector(), polyfem::solver::MatrixLagrangianForm::MatrixLagrangianForm(), polyfem::solver::PeriodicLagrangianForm::PeriodicLagrangianForm(), and polyfem::solver::BCLagrangianForm::update_target().
|
protected |
penalty parameter
Definition at line 44 of file AugmentedLagrangianForm.hpp.
Referenced by A_weight(), lagrangian_weight(), set_initial_weight(), polyfem::solver::BCLagrangianForm::update_lagrangian(), polyfem::solver::MacroStrainLagrangianForm::update_lagrangian(), polyfem::solver::MatrixLagrangianForm::update_lagrangian(), and polyfem::solver::PeriodicLagrangianForm::update_lagrangian().
|
private |
Definition at line 54 of file AugmentedLagrangianForm.hpp.
Referenced by A_weight(), L_weight(), and set_scale().
|
protected |
vector of lagrange multipliers
Definition at line 46 of file AugmentedLagrangianForm.hpp.
Referenced by polyfem::solver::BCLagrangianForm::first_derivative_unweighted(), polyfem::solver::MacroStrainLagrangianForm::first_derivative_unweighted(), polyfem::solver::MatrixLagrangianForm::first_derivative_unweighted(), polyfem::solver::PeriodicLagrangianForm::first_derivative_unweighted(), polyfem::solver::BCLagrangianForm::init_masked_lumped_mass(), polyfem::solver::MacroStrainLagrangianForm::MacroStrainLagrangianForm(), polyfem::solver::MatrixLagrangianForm::MatrixLagrangianForm(), polyfem::solver::PeriodicLagrangianForm::PeriodicLagrangianForm(), polyfem::solver::BCLagrangianForm::update_lagrangian(), polyfem::solver::MacroStrainLagrangianForm::update_lagrangian(), polyfem::solver::MatrixLagrangianForm::update_lagrangian(), polyfem::solver::PeriodicLagrangianForm::update_lagrangian(), polyfem::solver::BCLagrangianForm::value_unweighted(), polyfem::solver::MacroStrainLagrangianForm::value_unweighted(), polyfem::solver::MatrixLagrangianForm::value_unweighted(), and polyfem::solver::PeriodicLagrangianForm::value_unweighted().