Saving (serializing) internal state of Boost.Numeric.Odeint steppers with full precision

Question

I have an application where I need to interrupt ODE evolution after a certain time (or, number of steps), save the state of the system + the ODE engine with binary precision, so