Article ID | Journal | Published Year | Pages | File Type |
---|---|---|---|---|
4638087 | Journal of Computational and Applied Mathematics | 2016 | 10 Pages |
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
Miloud Sadkane, Roger B. Sidje,