https://hgpu.org/?p=5089
Optimized Parallel Implementation of Gillespie's First Reaction Method on Graphics Processing Units