We should use the iterator interface and nonlinear solve cache from NonlinearSolve.jl in the NonlinearSolveLocalSearch implementation (similar to how SimpleContinuation.jl does). Also, it's probably worth also investigating if a similar setup could be implemented for the LBFGSLocalSearch.