reuse the simplified-SDC integrator for the ODE update

[zingale]

this eliminates all of the custom VODE code here

PR summary

PR motivation

PR checklist

• test suite needs to be run on this PR
• this PR will change answers in the test suite to more than roundoff level
• all newly-added functions have docstrings as per the coding conventions
• the CHANGES file has been updated, if appropriate
• if appropriate, this change is described in the docs

[zingale]

this requires: AMReX-Astro/Microphysics#1341

[zingale]

this shows the expected 4th order convergence in reacting_convergence:

 $\rho$                           2.135967e+17   3.850       1.481741e+16   3.959       9.528251e+14
 $\rho u$                          3.41459e+25   3.849       2.369186e+24   3.948       1.535472e+23
 $\rho v$                          3.41459e+25   3.849       2.369186e+24   3.948       1.535472e+23
 $\rho E$                          1.95356e+34   3.884       1.323478e+33   3.938       8.637672e+31
 $\rho e$                         1.678975e+34   3.892       1.130798e+33   3.940       7.365531e+31
 $T$                              1.251912e+20   3.710       9.563941e+18   3.966       6.119475e+17
 $\rho X(\isotm{He}{4})$          2.156035e+17   3.852       1.492992e+16   3.957       9.610957e+14
 $\rho X(\isotm{C}{12})$          8.879796e+15   3.791       6.414778e+14   3.939       4.182066e+13
 $\rho X(\isotm{O}{16})$          5.791356e+12   3.747       4.312966e+11   3.915       2.859503e+10
 $\rho X(\isotm{Fe}{56})$             21359610   3.850            1481740   3.959           95282.51

[zingale]

this is ready now