Implement `NonlinearCG` variant which uses memoryless BFGS update

[patrick-kidger]

This would allow us to use all our line searches and descents with the nonlinear CG approximate Hessian. See Conjugate Gradient Methods with Inexact Line Search by Shanno.