Skip to content
Licensed Unlicensed Requires Authentication Published by De Gruyter October 19, 2011

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

  • Jesper Karlsson EMAIL logo and Raúl Tempone
From the journal


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.

Received: 2010-05-25
Revised: 2011-07-26
Published Online: 2011-10-19
Published in Print: 2011-September

© de Gruyter 2011

Downloaded on 27.3.2023 from
Scroll Up Arrow