Article ID Journal Published Year Pages File Type
4638087 Journal of Computational and Applied Mathematics 2016 10 Pages PDF
Abstract

We consider the problem of computing the spectral projections of a regular matrix pair associated to eigenvalues that are inside the unit circle. We discuss two algorithmic variants that both rely on the generalized Schur decomposition with re-ordering. From there, they both obtain a block-diagonalization but in different ways. One does so by solving a generalized Sylvester equation while the other uses a reverse ordering of the eigenvalues in the Schur decomposition. The block diagonalization allows accurate computation of the needed left and right spectral projections. Numerical comparisons show that this technique is a strong competitor to the spectral dichotomy approach.

Related Topics
Physical Sciences and Engineering Mathematics Applied Mathematics
Authors
, ,