Custom rule for differentiating through Newton solver using ForwardDiff: works for gradient, fails for hessian

Very cool, thanks! We’re currently reshaping the interface but as soon as that is done I’ll take a look