Completeness for use in OrdinaryDiffEq.jl

[ChrisRackauckas]

This package is almost complete to be used in OrdinaryDiffEq.jl. For other discussions on this topic, see SciML/OrdinaryDiffEq.jl#1570

However, there are some special things done in its nonlinear solvers which we would need to allow for here:

1. First step is completing the iterator interface Complete the iterator interface #241
2. We need to have a version of a Quasi-Newton method which forces Jacobian reuse and terminates if convergence starts to fail. I think this can be done just with a SingleFactorizeNewton and some extended TerminationCondition stuff.
3. For this SingleFactorizeNewton, we'd need to allow it in the iterator form to have a mode where changing u does not cause a new Jacobian to be taken. Thus in the cache interface, we'd need a reevalute_jacobian! that can be handled outside the solver for such cases.

Once those are together, then I think we're good to try it in OrdinaryDiffEq?

[ChrisRackauckas]

Related #206 @frankschae

[oscardssmith]

@ChrisRackauckas I believe this can now be closed, right?