You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Sometimes the rounding error might cause a linear solver that requires strict positive definiteness to crash. However the Eigen::SimplicialLDLT we use by default works just fine.
A simple trick that one can do is to add an epsilon * identitiyMatrix to each of the projected element Hessian, so that the rounding error won't make the global Hessian indefinite.
Similarly, one can also blend the global Hessian with an epsilon scaled lumped mass to help with this issue.
The text was updated successfully, but these errors were encountered:
Sometimes the rounding error might cause a linear solver that requires strict positive definiteness to crash. However the Eigen::SimplicialLDLT we use by default works just fine.
A simple trick that one can do is to add an epsilon * identitiyMatrix to each of the projected element Hessian, so that the rounding error won't make the global Hessian indefinite.
Similarly, one can also blend the global Hessian with an epsilon scaled lumped mass to help with this issue.
The text was updated successfully, but these errors were encountered: