You can use the integrator interface with reinit
.
https://diffeq.sciml.ai/stable/basics/integrator/
Or you can use the ensemble interface to auto-parallelize it.
https://diffeq.sciml.ai/stable/features/ensemble/
This includes automated GPU parallelism with DiffEqGPU.jl