Theoretical and Computable Optimal Subspace Expansions for Matrix Eigenvalue Problems
Abstract
Consider the optimal subspace expansion problem for the matrix eigenvalue problem Ax=λ x: Which vector w in the current subspace V, after multiplied by A, provides an optimal subspace expansion for approximating a desired eigenvector x in the sense that x has the smallest angle with the expanded subspace Vw=V+ span\Aw\, i.e., wopt=w∈V(Vw,x)? This problem is important as many iterative methods construct nested subspaces that successively expand V to Vw. An expression of wopt by Ye (Linear Algebra Appl., 428 (2008), pp. 911--918) for A general, but it could not be exploited to construct a computable (nearly) optimally expanded subspace. He turns to deriving a maximization characterization of (Vw,x) for a given w∈ V when A is Hermitian. We generalize Ye's maximization characterization to the general case and find its maximizer. Our main contributions consist of explicit expressions of wopt, (I-PV)Awopt and the optimally expanded subspace Vwopt for A general, where PV is the orthogonal projector onto V. These results are fully exploited to obtain computable optimally expanded subspaces within the framework of the standard, harmonic, refined, and refined harmonic Rayleigh--Ritz methods. We show how to efficiently implement the proposed subspace expansion approaches. Numerical experiments demonstrate the effectiveness of our computable optimal expansions.