Towards automatic global error control: Computable weak error expansion for the tau-leap method

Jesper Karlsson 1  and Raúl Tempone 1
  • 1 Mathematical and Computer Sciences and Engineering (MCSE), King Abdullah University of Science and Technology (KAUST), Saudi Arabia.

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 Slgorithm. 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.

Purchase article
Get instant unlimited access to the article.
$42.00
Log in
Already have access? Please log in.


or
Log in with your institution

Journal + Issues

This quarterly published journal presents original articles on the theory and applications of Monte Carlo and Quasi-Monte Carlo methods. Launched in 1995 the journal covers all stochastic numerics topics with emphasis on the theory of Monte Carlo methods and new applications in all branches of science and technology.

Search