Towards Automatic Global Error Control: Computable Weak Error Expansion for the Tau-Leap Method
Abstract
This work develops novel error expansions with computable leading order terms for the global weak error in the tau-leap discretization of pure jump processes arising in kinetic Monte Carlo models. Accurate computable a posteriori error approximations are the basis for adaptive algorithms; a fundamental tool for numerical simulation of both deterministic and stochastic dynamical systems. These pure jump processes are simulated either by the tau-leap method, or by exact simulation, also referred to as dynamic Monte Carlo, the Gillespie algorithm or the Stochastic simulation algorithm. Two types of estimates are presented: an a priori estimate for the relative error that gives a comparison between the work for the two methods depending on the propensity regime, and an a posteriori estimate with computable leading order term.
- Publication:
-
arXiv e-prints
- Pub Date:
- April 2010
- DOI:
- 10.48550/arXiv.1004.2948
- arXiv:
- arXiv:1004.2948
- Bibcode:
- 2010arXiv1004.2948K
- Keywords:
-
- Mathematics - Numerical Analysis;
- 60H10;
- 60H35;
- 60J75;
- 65Y20;
- 65L50
- E-Print:
- Monte Carlo Methods and Applications 2011 17:3, 233-278