The underlying OrdinaryDiffEq.jl package supports solving from a batch of initial conditions. It's very likely I am misunderstanding something but in principle it should be possible to support the same even here. Would love to know if there is a particular reason that the dynamics and readout functions have their type signatures restricted to AbstractVector.