Search

MiS Preprint Repository

We have decided to discontinue the publication of preprints on our preprint server as of 1 March 2024. The publication culture within mathematics has changed so much due to the rise of repositories such as ArXiV (www.arxiv.org) that we are encouraging all institute members to make their preprints available there. An institute's repository in its previous form is, therefore, unnecessary. The preprints published to date will remain available here, but we will not add any new preprints here.

MiS Preprint
80/2011

Fast solution of multi-dimensional parabolic problems in the TT/QTT-format with initial application to the Fokker-Planck equation

Sergey Dolgov, Boris N. Khoromskij and Ivan V. Oseledets

Abstract

In this paper we propose two schemes of using the QTT tensor approximations for solution of multi-dimensional parabolic problems. First, we present a simple one-step implicit time integration scheme and modify it using the matrix multiplication and a linear ALS-type solver in the TT format. As the second approach, we propose the global space-time formulation, resulting in a large block linear system, encapsulating all time steps, and solve it at once in the QTT format.

We prove the QTT-rank estimate for certain classes of multivariate potentials and respective solutions in (x,t) variables. We observe the log-linear complexity of storage and the solution algorithm in both spatial and time grid sizes, and at most cubic scaling in the QTT ranks of the discretized operator matrix and solution.

The method is applied to the Fokker-Planck equation arising from the beads-springs models of polymeric liquids. For the dumbbell model numerical experiments are shown to demonstrate logarithmic behavior of computational time versus number of grid points in space and time, as well as accuracy. However, in numerical tests for the case of multispring Hookean potential we observe, that the rank properties of more general models might make the straightforward application of the tensor product approximations inefficient, requiring modifications in model descriptions and tensor discretizations.

Received:
Dec 8, 2011
Published:
Dec 8, 2011
MSC Codes:
35K20, 65F50, 15A69, 65D15, 33F05, 65F10, 35Q84, 82D60
Keywords:
parabolic problems, QTT-format, DMRG, higher dimensions, tensor methods, fokker-planck equation, dumbbell model

Related publications

inJournal
2012 Repository Open Access
Sergey Dolgov, Boris N. Khoromskij and Ivan V. Oseledets

Fast solution of parabolic problems in the tensor train/quantized tensor train format with initial application to the Fokker-Planck equation

In: SIAM journal on scientific computing, 34 (2012) 6, A3016-A3038