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

06/12/2020
by   Koen Ruymbeek, et al.
0

In this paper we extend the Residual Arnoldi method for calculating an extreme eigenvalue (e.g. largest real part, dominant,...) 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 applies the Residual Arnoldi method (RA) for a grid of parameter points at the same time. 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.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
12/01/2020

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

In this paper we solve m-parameter eigenvalue problems (mEPs), with m an...
research
10/16/2020

A residual concept for Krylov subspace evaluation of the φ matrix function

An efficient Krylov subspace algorithm for computing actions of the φ ma...
research
12/05/2019

An accurate restarting for shift-and-invert Krylov subspaces computing matrix exponential actions of nonsymmetric matrices

An accurate residual–time (AccuRT) restarting for computing matrix expon...
research
06/12/2017

Asymptotic convergence of spectral inverse iterations for stochastic eigenvalue problems

We consider and analyze applying a spectral inverse iteration algorithm ...
research
10/03/2018

Enhanced image approximation using shifted rank-1 reconstruction

Low rank approximation has been extensively studied in the past. It is m...
research
08/03/2020

A parameter-dependent smoother for the multigrid method

The solution of parameter-dependent linear systems, by classical methods...
research
12/25/2018

ART: adaptive residual--time restarting for Krylov subspace matrix exponential evaluations

In this paper a new restarting method for Krylov subspace matrix exponen...

Please sign up or login with your details

Forgot password? Click here to reset