Complex Shift and Invert Strategies for Real Matrices

摘要

When using an iterative method for solving a generalized nonsymmetric eigenvalue problem of the form Fu = lambda Mu, where f and M are real matrices, it is foten desirable to work with the shifted and inverted operator B = (K -A-M) to the minus first power M in order to enhance the eigenvalue separation and improve efficiency. Unfortunately, the shift sigma is generally complex and so is the matrix B. The question then is whether it is possible to avoid complex arithmetic while preserving the advantages of bandedness of the pair (F,M). For the classical problem where M = I and F is banded, complex arithmetic can be avoided by using double shifts, i.e., by working with the real matrix BB whose bandwidth is double that of F. This satisfactory solution extends to the case where M is diagonal as well. In the generalized case the answer to the above question is negative, in the sense that complex arithmetic can be avoided only at the expense of loosing the advantage of bandedness. One solution is to factor the shifted matrix F - sigma M in complex arithmetic but employ real arithmetic subsequently in the iterative procedure. This paper examines several approaches and discusses their respective merits under different circumstances. (Author)