Figure 4
From: Simulation methods with extended stability for stiff biochemical Kinetics
Kullback-Leibler divergence for the Schlögl reaction. Kullback-Leibler divergence between the exact stationary distribution of X in the Schlögl reaction (estimated by 106 samples solved by SSA) and the approximate stationary distributions obtained with the Poisson τ-leap (black), Optimal RK s = 3; ϵ = 0.5 (grey lines) and Optimal RK s = 5, ϵ = 0.5 (white). Bars are shown only for the stable method and τ settings. Asterisks denote methods that have a rate of failure above 10-3.