6
$\begingroup$

I've been studying many iterative methods, like Jacobi, fixed-point, Newton's and the conjugate gradient methods.

Currently, I'm studying the CG method, but it's not the first time where the eigenvectors (and eigenvalues) of the usual matrix involved are important in determining how the iterative method behaves, like what's the convergence of the method.

In particular, I was reading the book "A first course in numerical methods" (by Greif, Ascher) and at a certain point (pp. 186-7) there's:

If the eigenvalues of A are located in only a few narrow clusters, then the CG method requires only a few iterations to converge. Convergence is slower, though, if the eigenvalues are widely spread.

which makes me wonder why, why if the eigenvalues of $A$ are located in an only a few narrow clusters the CG converges fast?

This is not the first time that I see similar observations where the eigenvalues and eigenvectors of the matrix involved in the iterative method somehow determine the behavior or performance of the iterative method.

My questions, apart from the specific case above, are:

  1. What are the relations between the eigenvectors and eigenvalues of the matrices involved in an iterative method and the behaviour of the iterative method (if this can be generalized)?

  2. What are the general behaviors of iterative methods that we can predict given the eigenvectors or eigenvalues?

If this can't be generalized, either I could ask specific questions, or you could like list a few examples where eigenvectors and eigenvalues influence different iterative methods. Of course, an exhaustive list would be nice for everyone.

Note: I know what are eigenvectors and eigenvalues, even though when I think about them I can't really make sense of why they are useful. You could argue that they are useful like in the example above, but I can't see the relation very well and clearly.

$\endgroup$
0
$\begingroup$

Iterative method can be divided in two big groups:

  • Stationary methods that they use an fixed point iteration idea
  • Krylov methods that they are projection methods. (in this question you are interested in these)

See this @Christian_Clason's answer for a more details.

What are the relations between the eigenvectors and eigenvalues of the matrices involved in an iterative method and the behaviour of the iterative method (if this can be generalized)?

What are general behaviours of iterative methods that we can predict given the eigenvectors or eigenvalues?

For the CG the relation with eigenvalues is given by the error estimate that involves the condition number $k= \lambda_{max} / \lambda_{min}$: $$ || x_\star - x_m||_A \leq 2 \left[ \frac{\sqrt{k} - 1}{\sqrt{k} + 1} \right]^m || x_\star - x_0||_A $$ See pag 187 of your book, or [1] pag 215 equation (6.128) here you can find also the proof.

Note:

  • the use of condition number for the estimate you can consider a "characteristic" of the Krylov methods.
  • for different algorithm you can find the opportune estimate, see [1]

Considering the case:

If the eigenvalues of A are located in only a few narrow clusters

then also $\lambda_{max}, \lambda_{min}$ are near and you have got a better $k$, the perfect case is when $k=1$. In the realistic case is when $k$ as near as possible to 1. Better is $k$ less iterations you need for the convergence.

Note: This is the scope of the precondition preconditioning technique, replace the original linear system with an equivalent whose spectral properties are better (see pag 187-88 of your book, or dedicate chapters of [1]). I.e. that the eigenvalues of the precondition matrix $P^{-1}A$ are in a cluster to 1.


[1] Saad, Yousef. Iterative methods for sparse linear systems. Siam, 2003.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.