HMDR and FMDR algorithms for the generalized eigenvalue problem

D. C. Dzeng*, Wen-Wei Lin

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

1 Scopus citations

Abstract

A new method which is based on two transformations, called the HMDR and the FMDR transformation, is presented for solving the generalized eigensystem Ax=λBx, where A,B are real square matrices. To overcome difficulties arising in the presence of infinite eigenvalues with high-order Jordan block, we introduce a new procedure which could be used as "preprocessing" to deflate the infinite eigenvalue. This method initially reduces B to a diagonal matrix with nonnegative entries and deflates the infinite eigenvalues if necessary. Maintaining the diagonal form of the right-hand-side matrix, A-λB is first reduced to Hessenberg diagonal form and then in an iterative part to quasitriangular diagonal form. In computational work it shows higher efficiency than the QZ algorithm, and in stability it is analogous to the MDR algorithm. Our method is closely related to the QR algorithm in essence.

Original languageEnglish
Pages (from-to)169-187
Number of pages19
JournalLinear Algebra and Its Applications
Volume112
Issue numberC
DOIs
StatePublished - 1 Jan 1989

Fingerprint Dive into the research topics of 'HMDR and FMDR algorithms for the generalized eigenvalue problem'. Together they form a unique fingerprint.

Cite this