Hello, I guess the computation of the solution of the equation involving the hessian is computationally the most expensive step in the optimization process. Is there a way of passing a function to calculate the inverse of the hessian instead of the hessian to the optimizer? Cheers, -Ajo.