The algebraic eigenvalue problem
The algebraic eigenvalue problem
LAPACK's user's guide
Computing selected eigenvalues of sparse unsymmetric matrices using subspace iteration
ACM Transactions on Mathematical Software (TOMS)
ACM Transactions on Mathematical Software (TOMS)
Algorithm 570: LOPSI: A Simultaneous Iteration Method for Real Matrices [F2]
ACM Transactions on Mathematical Software (TOMS)
Journal of Computational Physics
Hi-index | 0.01 |
SRRT is a Fortran program to calculate an approximate orthonomral basis fr a dominant invariant subspace of a real matrix A by the method of simultaneous iteration. Specifically, given an integer m, SRRIT computes a matrix Q with m orthonormal columns and real quasi-triangular matrix T or order m such that the equation AQ = QT is satisfied up to a tolerance specified by the user. The eigenvalues of T are approximations to the m eigenvalues of largest absolute magnitude of A and the columns of Q span the invariant subspace corresponding to those eigenvalues. SRRIT references A only through a user-provided subroutine to form the product AQ; hence it is suitable for large sparse problems.