Mathematics – Numerical Analysis
Scientific paper
2006-10-16
SIAM Journal on Matrix Analysis and Applications, Vol.30, No.2, pp. 548-559, 2008.
Mathematics
Numerical Analysis
12 pages. Accepted to SIAM Journal on Matrix Analysis and Applications (SIMAX)
Scientific paper
10.1137/070684628
The Rayleigh-Ritz method is widely used for eigenvalue approximation. Given a matrix $X$ with columns that form an orthonormal basis for a subspace $\X$, and a Hermitian matrix $A$, the eigenvalues of $X^HAX$ are called Ritz values of $A$ with respect to $\X$. If the subspace $\X$ is $A$-invariant then the Ritz values are some of the eigenvalues of $A$. If the $A$-invariant subspace $\X$ is perturbed to give rise to another subspace $\Y$, then the vector of absolute values of changes in Ritz values of $A$ represents the absolute eigenvalue approximation error using $\Y$. We bound the error in terms of principal angles between $\X$ and $\Y$. We capitalize on ideas from a recent paper [DOI: 10.1137/060649070] by A. Knyazev and M. Argentati, where the vector of absolute values of differences between Ritz values for subspaces $\X$ and $\Y$ was weakly (sub-)majorized by a constant times the sine of the vector of principal angles between $\X$ and $\Y$, the constant being the spread of the spectrum of $A$. In that result no assumption was made on either subspace being $A$-invariant. It was conjectured there that if one of the trial subspaces is $A$-invariant then an analogous weak majorization bound should only involve terms of the order of sine squared. Here we confirm this conjecture. Specifically we prove that the absolute eigenvalue error is weakly majorized by a constant times the sine squared of the vector of principal angles between the subspaces $\X$ and $\Y$, where the constant is proportional to the spread of the spectrum of $A$. For many practical cases we show that the proportionality factor is simply one, and that this bound is sharp. For the general case we can only prove the result with a slightly larger constant, which we believe is artificial.
Argentati Merico E.
Knyazev Andrew V.
Paige Christopher C.
Panayotov I.
No associations
LandOfFree
Bounds on changes in Ritz values for a perturbed invariant subspace of a Hermitian matrix does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with Bounds on changes in Ritz values for a perturbed invariant subspace of a Hermitian matrix, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Bounds on changes in Ritz values for a perturbed invariant subspace of a Hermitian matrix will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-390629