Solving chemical master equations by adaptive wavelet compression

JOURNAL OF COMPUTATIONAL PHYSICS

卷:229

Article

Jahnke, Tobias¹  Udrescu, Tudor¹ 

[1] Karlsruhe Inst Technol, Fak Math, Inst Angew & Numer Math, D-76133 Karlsruhe, Germany

关键词: Chemical master equation;    Wavelet compression;    Adaptive step-size selection;    Adaptive Galerkin approximation;    Rothe's method;    Stochastic reaction kinetics;    Gene regulatory networks;   

DOI  :  10.1016/j.jcp.2010.04.015

来源: Elsevier

【 摘 要 】

Solving chemical master equations numerically on a large state space is known to be a difficult problem because the huge number of unknowns is far beyond the capacity of traditional methods. We present an adaptive method which compresses the problem very efficiently by representing the solution in a sparse wavelet basis that is updated in each step. The step-size is chosen adaptively according to estimates of the temporal and spatial approximation errors. Numerical examples demonstrate the reliability of the error estimation and show that the method can solve large problems with bimodal solution profiles. (C) 2010 Elsevier Inc. All rights reserved.

【 授权许可】

Free