Homotopy continuation method for the numerical solutions of generalised symmetric eigenvalue problems

Abstract We consider a generalised symmetric eigenvalue problem Ax = λMx, where A and M are real n by n symmetric matrices such that M is positive semidefinite. The purpose of this paper is to develop an algorithm based on the homotopy methods in [9, 11] to compute all eigenpairs, or a specified number of eigenvalues, in any part of the spectrum of the eigenvalue problem Ax = λMx. We obtain a special Kronecker structure of the pencil A − λM, and give an algorithm to compute the number of eigenvalues in a prescribed interval. With this information, we can locate the lost eigenpair by using the homotopy algorithm when multiple arrivals occur. The homotopy maintains the structures of the matrices A and M (if any), and the homotopy curves are n disjoint smooth curves. This method can be used to find all/some isolated eigenpairs for large sparse A and M on SIMD machines.

[1]  J. Milnor.  Topology from the differentiable viewpoint , 1965 .

[2]  G. Fix,et al.  An Algorithm for the Ill-Conditioned Generalized Eigenvalue Problem , 1972 .

[3]  G. Stewart,et al.  An Algorithm for Generalized Matrix Eigenvalue Problems. , 1973 .

[4]  J. Bunch,et al.  Some stable methods for calculating inertia and solving symmetric linear systems , 1977 .

[5]  J. H. Wilkinson,et al.  Inverse Iteration, Ill-Conditioned Equations and Newton’s Method , 1979 .

[6]  B. Parlett.  The Symmetric Eigenvalue Problem , 1981 .

[7]  Leigh Tesfatsion,et al.  Individual tracking of an eigenvalue and eigenvector of a parameterized matrix , 1981 .

[8]  Leigh Tesfatsion,et al.  Variational equations for the eigenvalues and eigenvectors of nonsymmetric matrices , 1981 .

[9]  Peter B. Geltner.  General Rayleigh Quotient Iteration , 1981 .

[10]  Angelika Bunse-Gerstner,et al.  An algorithm for the symmetric generalized eigenvalue problem , 1984 .

[11]  M. Chu.  A simple application of the homotopy method to symmetric eigenvalue problems , 1984 .

[12]  Tien-Yien Li,et al.  Homotopy method for generalized eigenvalue problems Ax= ΛBx☆ , 1987 .

[13]  Tien-Yien Li,et al.  Homotopy method for general l-matrix problems , 1988 .

[14]  Wen-Wei Lin,et al.  An application of the homotopy method to the generalised symmetric eigenvalue problem , 1988, The Journal of the Australian Mathematical Society. Series B. Applied Mathematics.

[15]  Tien-Yien Li,et al.  Homotopy algorithm for symmetric eigenvalue problems , 1989 .