Preconditioned Low-Rank Methods for High-Dimensional Elliptic PDE Eigenvalue Problems

Abstract We consider elliptic PDE eigenvalue problems on a tensorized domain, discretized such that the resulting matrix eigenvalue problem Ax=λx exhibits Kronecker product structure. In particular, we are concerned with the case of high dimensions, where standard approaches to the solution of matrix eigenvalue problems fail due to the exponentially growing degrees of freedom. Recent work shows that this curse of dimensionality can in many cases be addressed by approximating the desired solution vector x in a low-rank tensor format. In this paper, we use the hierarchical Tucker decomposition to develop a low-rank variant of LOBPCG, a classical preconditioned eigenvalue solver. We also show how the ALS and MALS (DMRG) methods known from computational quantum physics can be adapted to the hierarchical Tucker decomposition. Finally, a combination of ALS and MALS with LOBPCG and with our low-rank variant is proposed. A number of numerical experiments indicate that such combinations represent the methods of choice.

[1]  Daniel Kressner,et al.  Low-Rank Tensor Krylov Subspace Methods for Parametrized Linear Systems , 2011, SIAM J. Matrix Anal. Appl..

[2]  Gene H. Golub,et al.  Matrix computations , 1983 .

[3]  Vladimir A. Kazeev,et al.  On explicit QTT representation of Laplace operator and its inverse , 2010 .

[4]  Andrew V. Knyazev,et al.  Toward the Optimal Preconditioned Eigensolver: Locally Optimal Block Preconditioned Conjugate Gradient Method , 2001, SIAM J. Sci. Comput..

[5]  Dietrich Braess,et al.  F ¨ Ur Mathematik in Den Naturwissenschaften Leipzig Approximation of 1/x by Exponential Sums in [1, ∞) , 2022 .

[6]  Erwan Faou,et al.  Computing Semiclassical Quantum Dynamics with Hagedorn Wavepackets , 2009, SIAM J. Sci. Comput..

[7]  Lars Grasedyck,et al.  F ¨ Ur Mathematik in Den Naturwissenschaften Leipzig a Projection Method to Solve Linear Systems in Tensor Format a Projection Method to Solve Linear Systems in Tensor Format , 2022 .

[8]  U. Schollwoeck.  The density-matrix renormalization group in the age of matrix product states , 2010, 1008.3477.

[9]  F. Verstraete,et al.  Complete-graph tensor network states: a new fermionic wave function ansatz for molecules , 2010, 1004.5303.

[10]  Karl Meerbergen,et al.  Inverse Iteration for Purely Imaginary Eigenvalues with Application to the Detection of Hopf Bifurcations in Large-Scale Problems , 2010, SIAM J. Matrix Anal. Appl..

[11]  B. Khoromskij,et al.  Tensor-structured Galerkin approximation of parametric and stochastic elliptic PDEs , 2010 .

[12]  Yuji Nakatsukasa,et al.  Perturbation of Partitioned Hermitian Definite Generalized Eigenvalue Problems , 2011, SIAM J. Matrix Anal. Appl..

[13]  Hans-Dieter Meyer,et al.  A numerical study on the performance of the multiconfiguration time-dependent Hartree method for density operators , 2000 .

[14]  F. Verstraete,et al.  Matrix product operator representations , 2008, 0804.3976.

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

[16]  Jack Dongarra,et al.  Templates for the Solution of Algebraic Eigenvalue Problems , 2000, Software, environments, tools.

[17]  Lars Grasedyck,et al.  Hierarchical Singular Value Decomposition of Tensors , 2010, SIAM J. Matrix Anal. Appl..

[18]  Karl Meerbergen,et al.  Shift-and-invert iteration for purely imaginary eigenvalues with application to the detection of Hopf Bifurcations in large scale problems , 2008 .

[19]  Cleve B. Moler,et al.  Nineteen Dubious Ways to Compute the Exponential of a Matrix, Twenty-Five Years Later , 1978, SIAM Rev..

[20]  W. Hackbusch,et al.  A New Scheme for the Tensor Representation , 2009 .

[21]  Reinhold Schneider,et al.  On manifolds of tensors of fixed TT-rank , 2012, Numerische Mathematik.

[22]  T. Schulte-Herbrüggen,et al.  Computations in quantum tensor networks , 2012, 1212.5005.

[23]  Lars Grasedyck,et al.  Existence and Computation of Low Kronecker-Rank Approximations for Large Linear Systems of Tensor Product Structure , 2004, Computing.

[24]  Boris N. Khoromskij,et al.  Use of tensor formats in elliptic eigenvalue problems , 2012, Numer. Linear Algebra Appl..

[25]  Daniel Kressner,et al.  Krylov Subspace Methods for Linear Systems with Tensor Product Structure , 2010, SIAM J. Matrix Anal. Appl..

[26]  U. Manthe,et al.  The multi-configurational time-dependent Hartree approach , 1990 .

[27]  R. Schneider,et al.  The Alternating Linear Scheme for Tensor Optimisation in the TT Format , 2022 .