2000 年 2000 巻 p. 20000020
This paper presents a PCG method using the Sylvester law of inertia for the eigensolution of large, sparse and symmetric matrices. The proposed method, retaining the advantages of the conjugate gradient method, permits to count the number of sign changes for given matrices by the Sylvester law of inertia, and is able to overcome the numerical difficulty caused in the case where the solution converges to the true eigenvalue. This method is particularly useful to find only small numbers of lower eigenpairs in the large sparse system. The accuracy and stability of this method are confirmed by using several numerical examples. The numerical results give a good agreement even in the systems with multiple eigenvalues.