Hi,
I'm looking for a function in R that could solve an augmented Lagrangian function (with equality and inequality constraints) that takes the Hessian matrix to compute Newton steps.
So far, I have tried using constrOptim.nl
from the alabama
package, but it does not allow to use second order derivates for faster implementation. I have also tried using auglag
from the same package, but the result is the same.
There is a post on StackOverflow which refers to the trust
package, but this one does not allow for constrained optimization.
Thanks in advance !