Efficient step size selection for the tau-leaping simulation method
Efficient step size selection for the tau-leaping simulation method
复制标题
DOI:
10.1063/1.2159468
复制
发表时间:
2006-01-28
影响因子:
4.4
通讯作者:
Petzold, LR
中科院分区:
文献类型:
--
作者:
Cao, Y;Gillespie, DT;Petzold, LR
The tau-leaping method of simulating the stochastic time evolution of a well-stirred chemically reacting system uses a Poisson approximation to take time steps that leap over many reaction events. Theory implies that tau leaping should be accurate so long as no propensity function changes its value "significantly" during any time step tau. Presented here is an improved procedure for estimating the largest value for tau that is consistent with this condition. This new tau-selection procedure is more accurate, easier to code, and faster to execute than the currently used procedure. The speedup in execution will be especially pronounced in systems that have many reaction channels. (c) 2006 American Institute of Physics.