Changing problem and vector size with a callback during a OrdinaryDiffEq time integration retains the original solution vector
I am using OrdinaryDiffEq to solve a set of ordinary differential equations. During the time integration I would like to extend the solution vector and change the problem. For that I use a callback of type DiscreteCallback and reinitialize the integrator with the new problem: