QTT-isogeometric solver in two dimensions

The goal of this paper is to develop a numerical algorithm that solves a two-dimensional elliptic partial differential equation in a polygonal domain using tensor methods and ideas from isogeometric analysis. The proposed algorithm is based on the Finite Element (FE) approximation with Quantized Tensor Train decomposition (QTT) used for matrix representation and solution approximation. In this paper we propose a special discretisation scheme that allows to construct the global stiffness matrix in the QTT-format. The algorithm has $O(\log n)$ complexity, where $n=2^d$ is the number of nodes per quadrangle side. A new operation called z-kron is introduced for QTT-format. It makes it possible to build a matrix in z-order if the matrix can be expressed in terms of Kronecker products and sums. An algorithm for building a QTT coefficient matrix for FEM in z-order "on the fly", as opposed to the transformation of a calculated matrix into QTT, is presented. This algorithm has $O(\log n)$ complexity for $n$ as above.

[1]  Claude Brezinski,et al.  Extrapolation algorithms and Pade´ approximations: a historical survey , 1996 .

[2]  Anders Logg,et al.  Automated Solution of Differential Equations by the Finite Element Method: The FEniCS Book , 2012 .

[3]  S. V. Dolgov,et al.  ALTERNATING MINIMAL ENERGY METHODS FOR LINEAR SYSTEMS IN HIGHER DIMENSIONS∗ , 2014 .

[4]  Ivan V. Oseledets,et al.  Approximation of 2d˟2d Matrices Using Tensor Decomposition , 2010, SIAM J. Matrix Anal. Appl..

[5]  Bruno Lévy,et al.  Quad‐Mesh Generation and Processing: A Survey , 2013, Comput. Graph. Forum.

[6]  Boris N. Khoromskij,et al.  Simultaneous state-time approximation of the chemical master equation using tensor product formats , 2015, Numer. Linear Algebra Appl..

[7]  Ted D. Blacker,et al.  Paving: A new approach to automated quadrilateral mesh generation , 1991 .

[8]  B. Khoromskij.  O(dlog N)-Quantics Approximation of N-d Tensors in High-Dimensional Numerical Modeling , 2011 .

[9]  Christoph Schwab,et al.  Low-rank tensor structure of linear diffusion operators in the TT and QTT formats☆ , 2013 .

[10]  Ivan Oseledets,et al.  Tensor-Train Decomposition , 2011, SIAM J. Sci. Comput..

[11]  E. Liu,et al.  Fundamental Methods of Extrapolation Fundamental Methods of Numerical Extrapolation With Applications , 2006 .

[12]  Eugene E. Tyrtyshnikov,et al.  Breaking the Curse of Dimensionality, Or How to Use SVD in Many Dimensions , 2009, SIAM J. Sci. Comput..

[13]  Markus Bachmayr,et al.  Stability of Low-Rank Tensor Representations and Structured Multilevel Preconditioning for Elliptic PDEs , 2018, Foundations of Computational Mathematics.

[14]  VLADIMIR A. KAZEEV,et al.  Low-Rank Explicit QTT Representation of the Laplace Operator and Its Inverse , 2012, SIAM J. Matrix Anal. Appl..

[15]  G. Strang,et al.  An Analysis of the Finite Element Method , 1974 .

[16]  Vladimir A. Kazeev,et al.  Quantized tensor-structured finite elements for second-order elliptic PDEs in two dimensions , 2018, Numerische Mathematik.

[17]  Ivan V. Oseledets,et al.  Solution of Linear Systems and Matrix Inversion in the TT-Format , 2012, SIAM J. Sci. Comput..

[18]  Jennifer Seberry,et al.  The Strong Kronecker Product , 1994, J. Comb. Theory, Ser. A.