Memory-efficient Arnoldi algorithms for linearizations of matrix polynomials in Chebyshev basis

Abstract

Novel memory-efficient Arnoldi algorithms for solving matrix polynomial eigenvalue problems are presented. More specifically, we consider the case of matrix polynomials expressed in the Chebyshev basis, which is often numerically more appropriate than the standard monomial basis for a larger degree $d$. The standard way of solving polynomial eigenvalue problems proceeds by linearization, which increases the problem size by a factor $d$. Consequently, the memory requirements of Krylov subspace methods applied to the linearization grow by this factor. In this paper, we develop two variants of the Arnoldi method that build the Krylov subspace basis implicitly, in a way that only vectors of length equal to the size of the original problem need to be stored. The proposed variants are generalizations of the so called Q-Arnoldi and TOAR methods, which have been developed for the monomial case. We also show how the typical ingredients of a full implementation of the Arnoldi method, including shift-and-invert and restarting, can be incorporated. Numerical experiments are presented for matrix polynomials up to degree $30$ arising from the interpolation of nonlinear eigenvalue problems which stem from boundary element discretizations of PDE eigenvalue problems.