reintroduce the idea of multiple integrators

[zingale]

With the integrators now templated, it should be straightforward to create a wrapper that can call multiple different integrators depending on the thermodynamic conditions. For example, we could try RKC if it seems appropriate and VODE otherwise.