Lambda admissible subspaces of self adjoint matrices

Abstract

Given a self-adjoint matrix A and an index h such that λh(A) lies in a cluster of eigenvalues of A, we introduce the novel class of -admissible subspaces of A of dimension h. First, we show that the low-rank approximation of the form PT A PT, for a subspace T that is close to any -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 -admissible subspaces. We obtain upper bounds for the distance between subspaces obtained by the Rayleigh-Ritz method applied to A and the class of -admissible subspaces. We also find upper bounds for the condition number of the (set-valued) map computing the class of -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.

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…