English

Subspace method for multiparameter-eigenvalue problems based on tensor-train representations

Numerical Analysis 2020-12-03 v1 Numerical Analysis

Abstract

In this paper we solve mm-parameter eigenvalue problems (mmEPs), with mm any natural number by representing the problem using Tensor-Trains (TT) and designing a method based on this format. mmEPs typically arise when separation of variables is applied to separable boundary value problems. Often, methods for solving mmEP are restricted to m=3m = 3, due to the fact that, to the best of our knowledge, no available solvers exist for m>3m>3 and reasonable size of the involved matrices. In this paper, we prove that computing the eigenvalues of a mmEP can be recast into computing the eigenvalues of TT-operators. We adapted the algorithm in \cite{Dolgov2014a} for symmetric eigenvalue problems in TT-format to an algorithm for solving generic mmEPs. This leads to a subspace method whose subspace dimension does not depend on mm, in contrast to other subspace methods for mmEPS. This allows us to tackle mmEPs with m>3m > 3 and reasonable size of the matrices. We provide theoretical results and report numerical experiments. The MATLAB code is publicly available.

Keywords

Cite

@article{arxiv.2012.00815,
  title  = {Subspace method for multiparameter-eigenvalue problems based on tensor-train representations},
  author = {Koen Ruymbeek and Karl Meerbergen and Wim Michiels},
  journal= {arXiv preprint arXiv:2012.00815},
  year   = {2020}
}

Comments

17 pages, 14 figures