Lambda admissible subspaces of self adjoint matrices

Abstract

Given a self-adjoint matrix $A$ and an index $h$ such that $\hat{I}{»}_h(A)$ lies in a cluster of eigenvalues of $A$, we introduce the novel class of $\hat{I}›$-admissible subspaces of $A$ of dimension $h$. First, we show that the low-rank approximation of the form $P_{\mathcal{T}}A{P}_{\mathcal{T}}$, for a subspace $\mathcal{T}$ that is close to any $\hat{I}›$-admissible subspace of $A$, has nice properties. Then, we prove that some well-known iterative algorithms (such as the Subspace Iteration Method, or the Krylov subspace method) produce subspaces that become arbitrarily close to $\hat{I}›$-admissible subspaces. We obtain upper bounds for the distance between subspaces obtained by the Rayleigh-Ritz method applied to $A$ and the class of $\hat{I}›$-admissible subspaces. We also find upper bounds for the condition number of the (set-valued) map computing the class of $\hat{I}›$-admissible subspaces of $A$. Finally, we include numerical examples that show the advantage of considering this new class of subspaces in the clustered eigenvalue setting.

32 pages, 5 figures