I would like to use **scipy.integrate.ode** (or scipy.integrate.odeint) instances in **multiple threads** (one for each CPU core) in order to solve multiple IVPs at a time. However the documentation says: "*This integrator is not re-entrant. You cannot have two ode instances using the “vode” integrator at the same time.*"

(Also odeint causes internal errors if instantiated multiple times although the documentation does not say so.)

Any idea what can be done?