The constrained conjugate gradient algorithm allows the specification of a Hessian product function. The optional args parameter is currently not passed to the Hessian product, because the code checks the signature (number of arguments) of the function and adjust the parameters.
When cleaning this up, we should consider moving to scipy's LinearOperator for passing in the Hessian product.
The constrained conjugate gradient algorithm allows the specification of a Hessian product function. The optional
argsparameter is currently not passed to the Hessian product, because the code checks the signature (number of arguments) of the function and adjust the parameters.When cleaning this up, we should consider moving to scipy's
LinearOperatorfor passing in the Hessian product.