User's guide
fast_random_choice()
Developer's guide
About
Fast random choice from a list of weights.
This is at least for small arrays much faster than numpy.random.choice. For the Gillespie, overall this brings for 3 reactions a speedup factor of 2.