2018/02/01 by Jana Lipková, Lipková, Jana, Georgios Arampatzis +7
Business, Management and Accounting · Decision Sciences · Physics and Astronomy · #Advanced Queuing Theory Analysis #Advanced Thermodynamics and Statistical Mechanics #FOS: Biological sciences #FOS: Mathematics #Molecular Networks (q-bio.MN) #Probability (math.PR) #Simulation Techniques and Applications
paper · pdf · doi:10.48550/arxiv.1802.00296
openalex publication_date 2018/02/01 · openalex created_date 2024/04/11 · openalex updated_date 2026/07/28
We propose the S-leaping algorithm for the acceleration of Gillespie's stochastic simulation algorithm that combines the advantages of the two main accelerated methods; the τ-leaping and R-leaping algorithms. These algorithms are known to be efficient under different conditions; the τ-leaping is efficient for non-stiff systems or systems with partial equilibrium, while the R-leaping performs better in stiff system thanks to an efficient sampling procedure. However, even a small change in a system's set up can critically affect the nature of the simulated system and thus reduce the efficiency of an accelerated algorithm. The proposed algorithm combines the efficient time step selection from the τ-leaping with the effective sampling procedure from the R-leaping algorithm. The S-leaping is shown to maintain its efficiency under different conditions and in the case of large and stiff systems or systems with fast dynamics, the S-leaping outperforms both methods. We demonstrate the performance and the accuracy of the S-leaping in comparison with the τ-leaping and R-leaping on a number of benchmark systems involving biological reaction networks.