Tensor-product approach to global time-space-parametric discretization of chemical master equation
Sergey Dolgov and Boris N. Khoromskij
Contact the author: Please use for correspondence this email.
Submission date: 12. Nov. 2012 (revised version: November 2012)
MSC-Numbers: 65F50, 15A69, 65F10, 82C31, 80A30, 34B08
Keywords and phrases: multilinear algebra, tensor products, chemical master equation, parameter dependent problems
Download full preprint: PDF (978 kB)
We study the application of the novel tensor formats (TT, QTT, QTT-Tucker) to the solution of d-dimensional chemical master equations, applied mostly to gene regulating networks (signaling cascades, toggle switches, phage-λ). For some important cases, e.g. signaling cascade models, we prove good separability properties of the system operator. The time is treated as an additional variable, with the Quantized tensor representations (QTT, QTT-Tucker) employed, leading to the log-complexity in the system size. This global space-time (d + 1)-dimensional system, approximated in the QTT or QTT-Tucker formats, is solved in the block-diagonal form by the ALS-type iterations. Another issue considered is the quantification of uncertainty, which means that some model parameters are not known exactly, but only their ranges can be estimated. It occurs frequently in real-life systems. In this case, we introduce the unknown parameters as auxiliary variables discretized on the corresponding grids, and solve the global space-parametric system at once in the tensor formats.