Simulators should support infinite reaction rates.

[kaizu]

import numpy

with reaction_rules():
    A + B > C | 0.01
    C > D | numpy.inf

The second reaction must occur immediately after the first reaction.

[kaizu]

It is hard to support an infinite rate in ode. egfrd seems supporting it already, but not correctly working with multiple reactions with infinite rates.

with reaction_rules():
    A + B > C | 0.01
    C > D | numpy.inf
    C > E | numpy.inf

gillespie and meso work well even in the case now.