Tensor-Krylov method for computing eigenvalues of parameter-dependent matrices

• Published in: Journal of computational and applied mathematics Vol. 408; p. 113869
• Main Authors: Ruymbeek, Koen, Meerbergen, Karl, Michiels, Wim
• Format: Journal Article
• Language: English
• Published: Elsevier B.V 01-07-2022
• Subjects: Eigenvalues; Parameter-dependent matrices; Tensor-Krylov; Eigenvalues; 65F15; Tensor-Krylov; Parameter-dependent matrices
• ISSN: 0377-0427; 1879-1778
• DOI: 10.1016/j.cam.2021.113869

In this paper we extend the Residual Arnoldi method for calculating an extreme eigenvalue (e.g. largest real part, dominant, etc.) to the case where the matrices depend on parameters. The difference between this Arnoldi method and the classical Arnoldi algorithm is that in the former the residual is added to the subspace. We develop a Tensor-Krylov method that can be interpreted as an application of the Residual Arnoldi algorithm to a set of matrices obtained by evaluating a parameter-dependent matrix in parameter values on a grid. The subspace contains an approximate Krylov space for all these points. Instead of adding the residuals for all parameter values to the subspace we create a low-rank approximation of the matrix consisting of these residuals and add only the column space to the subspace. In order to keep the computations efficient, it is needed to limit the dimension of the subspace and to restart once the subspace has reached the prescribed maximal dimension. The novelty of this approach is twofold. Firstly, we observed that a large error in the low-rank approximations is allowed without slowing down the convergence, which implies that we can do more iterations before restarting. Secondly, we pay particular attention to the way the subspace is restarted, since classical restarting techniques give a too large subspace in our case. We motivate why it is good enough to just keep the approximation of the searched eigenvector. At the end of the paper we extend this algorithm to shift-and-invert Residual Arnoldi method to calculate the eigenvalue close to a shift σ for a specific parameter dependency. We provide theoretical results and report numerical experiments. The Matlab code is publicly available.