Jacobian-Free Krylov Versions for TR/LM/GN

[avik-pal]

• [x] Trust Region
• [x] Levenberg Marquardt
  • [x] LM can't support Krylov versions due to the $D^TD$ updates
  • [x] LinearSolve aliases u. Should we allocate 2 caches? Inefficient for factorizations where we want to reuse factorization. How do we prevent the aliasing
• [x] Gauss Newton
  • [x] Non-Square JacVec https://github.com/JuliaDiff/SparseDiffTools.jl/pull/274
• [x] Allow external control of the AD
• [x] Testing
• [x] Default to Zygote for Line Search if Zygote is loaded

Fixes #167

[codecov[bot]]

Codecov Report

Attention: 9 lines in your changes are missing coverage. Please review.

Comparison is base (0026bc1) 49.06% compared to head (bcfcc16) 94.04%.

Files	Patch %	Lines
src/jacobian.jl	91.52%	5 Missing :warning:
src/linesearch.jl	72.72%	3 Missing :warning:
src/NonlinearSolve.jl	0.00%	1 Missing :warning:

Additional details and impacted files

```diff @@ Coverage Diff @@ ## master #282 +/- ## =========================================== + Coverage 49.06% 94.04% +44.97% =========================================== Files 19 20 +1 Lines 1824 1896 +72 =========================================== + Hits 895 1783 +888 + Misses 929 113 -816 ```

:umbrella: View full report in Codecov by Sentry.
:loudspeaker: Have feedback on the report? Share it here.

[avik-pal]

@ChrisRackauckas this is good to go on my end