A fast exact simulation method for a class of Markov jump processes
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, Massachusetts 10003 (United States)
- School of Mathematics, Georgia Institute of Technology, Atlanta, Georgia 30332 (United States)
A new method of the stochastic simulation algorithm (SSA), named the Hashing-Leaping method (HLM), for exact simulations of a class of Markov jump processes, is presented in this paper. The HLM has a conditional constant computational cost per event, which is independent of the number of exponential clocks in the Markov process. The main idea of the HLM is to repeatedly implement a hash-table-like bucket sort algorithm for all times of occurrence covered by a time step with length τ. This paper serves as an introduction to this new SSA method. We introduce the method, demonstrate its implementation, analyze its properties, and compare its performance with three other commonly used SSA methods in four examples. Our performance tests and CPU operation statistics show certain advantages of the HLM for large scale problems.
- OSTI ID:
- 22493210
- Journal Information:
- Journal of Chemical Physics, Vol. 143, Issue 18; Other Information: (c) 2015 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9606
- Country of Publication:
- United States
- Language:
- English
Similar Records
An adaptive multi-level simulation algorithm for stochastic biological systems
Annealed Importance Sampling Reversible Jump MCMC algorithms