Short recurrences for computing extended Krylov bases for Hermitian and unitary matrices

Abstract

It is well known that the projection of a matrix A onto a Krylov subspace results in a matrix of Hessenberg form. We show that the projection of the same matrix A onto an extended Krylov subspace, i.e. a succesion of positive and negative powers of A multiplied with a starting vector, is a matrix of so-called extended Hessenberg form which can be characterized uniquely by its QR-factorization. In case A is a Hermitian or unitary matrix, this extended Hessenberg matrix is banded, resulting in short recurrence relations. For the unitary case, coupled two term recurrence relations are derived of which the coefficients capture all information necessary for a sparse factorization of the corresponding extended Hessenberg matrix. This generalizes the approach used by Watkins to retrieve the CMV-form for unitary matrices.nrpages: 23status: publishe