DifferentialEquations is getting support for integrating ODE parameters ensembles on the GPU, in parallel. Another example for an ‘embarassingly parallel’ problem is stochastic simulation with multiple copies of random numbers, to be averaged over later.
Is it feasible to do a Gillespie simulation on the GPU in Julia? What has to be done to get independent parallel streams of random numbers on the GPU? Does/will DiffEqGPU support this?