English

Lambda admissible subspaces of self adjoint matrices

Numerical Analysis 2026-02-13 v1 Numerical Analysis

Abstract

Given a self-adjoint matrix AA and an index hh such that λh(A)\lambda_h(A) lies in a cluster of eigenvalues of AA, we introduce the novel class of Λ\Lambda-admissible subspaces of AA of dimension hh. First, we show that the low-rank approximation of the form PTAPTP_{\mathcal{T}} A P_{\mathcal{T}}, for a subspace T\mathcal{T} that is close to any Λ\Lambda-admissible subspace of AA, 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 Λ\Lambda-admissible subspaces. We obtain upper bounds for the distance between subspaces obtained by the Rayleigh-Ritz method applied to AA and the class of Λ\Lambda-admissible subspaces. We also find upper bounds for the condition number of the (set-valued) map computing the class of Λ\Lambda-admissible subspaces of AA. Finally, we include numerical examples that show the advantage of considering this new class of subspaces in the clustered eigenvalue setting.

Cite

@article{arxiv.2602.11976,
  title  = {Lambda admissible subspaces of self adjoint matrices},
  author = {Francisco Arrieta Zuccalli and Pedro Massey},
  journal= {arXiv preprint arXiv:2602.11976},
  year   = {2026}
}

Comments

32 pages, 5 figures

R2 v1 2026-07-01T10:33:42.717Z