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
At the moment the Hessian is computed as: H = J.transpose() * J; (1)
this can be optimized using: _H.selfadjointView<Eigen::Lower>().rankUpdate(_J.transpose()); (2)
however, we must pay attention to sparsity. Basically, the solver needs to be initialized with the sparsity from (1) but then in the solve we need to compute the Hessian as in (2), otherwise errors rise and it only works when the Hessian is diagonal.
The text was updated successfully, but these errors were encountered:
At the moment the Hessian is computed as:
H = J.transpose() * J;
(1)this can be optimized using:
_H.selfadjointView<Eigen::Lower>().rankUpdate(_J.transpose());
(2)however, we must pay attention to sparsity. Basically, the solver needs to be initialized with the sparsity from (1) but then in the
solve
we need to compute the Hessian as in (2), otherwise errors rise and it only works when the Hessian is diagonal.The text was updated successfully, but these errors were encountered: