Simulation-level parallelism

SSAs are "embarrassingly" parallel on the simulation level, meaning that N independent simulations can be run at the same time. I don't know how high N is in practice for the users of JumpProcesses but if it's in the thousands, and each simulation is small (e.g. <100 reactions), perhaps each simulation can be run by its own core of a GPU. Memory management would be hard and the networks would need to be small enough to fit in memory, but the potential speedup can be 100-1000x. Are there any fundamental obstacles to this? Perhaps making a prototype with Direct can be a GSoC project for someone already familiar with GPU programming.

The issue names no files, tests, or entry points. Start by determining whether JumpProcesses.jl supports running many small Direct simulations concurrently on a GPU, then assess the memory constraints and fundamental obstacles described in the issue. Done would be a justified feasibility result or a scoped prototype proposal.