Michael Zingale

some compilation failures: http://groot.astro.sunysb.edu/Microphysics/test-suite/gfortran/2021-10-07-003/index.html

note that that is what VODE is -- it is BDF. It is just a different implementation.

this seems to work, although it comes with the overhead of a backup `burn_t`

CPU tests pass: http://groot.astro.sunysb.edu/Microphysics/test-suite/gfortran/2022-09-07-008/index.html

note: this is for the analytic Jacobian

this can integrate aprox13 it only uses rtol_spec and atol_spec for error estimation.

Here's a comparison using burn_cell between VODE and Backward Euler: 

yes, I just need to make that not depend on VODE. I'll do that in a separate PR.

I think I might have a first cut in a branch somewhere...

I'll add a test today