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

We just added support for be forwarddiff as an extension under 1.9. Can you see if it works as desired for your use case?