Solving Ax=b, when b is an eigenvector

  • Thread starter Thread starter azdang
  • Start date Start date
  • Tags Tags
    Eigenvector
Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
3 replies · 2K views
azdang
Messages
80
Reaction score
0
Let A be a nonsingular matrix. What can you say about the convergence of GMRES to the solution of Ax=b when b is an eigenvector of A?

I know that if A is a nonsingular matrix with minimal polynomial of degree m, then GMRES solves Ax=b in at most m iterations.

Because b is an eigenvector of A then: (A-[tex]\lambda[/tex]I)b=0.

The only thing I've sort of come up with is that because b is an eigenvector, it cannot equal 0. At first I though that that implies that A-[tex]\lambda[/tex]I=0, but I'm thinking that that is not necessarily true. However, IF it is true that A-[tex]\lambda[/tex]I=0, I completely understand, because this would imply that the minimal polynomial is q(x)=x-[tex]\lambda[/tex]. This would be a minimal polynomial of degree 1, so GMRES would converge in 1 iteration. Can anyone help me out because I don't think this reasoning is correct, but I'm thinking that A-[tex]\lambda[/tex] is not necessarily nonsingular...or is it? Thank you!
 
Physics news on Phys.org
Can anyone offer any insight?
 
I just read the Wikipedia page on GMRES. Perhaps I misinterpret, but it seems like GMRES sould solve [tex]Ax = b[/tex] in just one step if [tex]b[/tex] is an eigenvector, since the exact solution [tex]x = b/\lambda[/tex] is in the span of [tex]b[/tex] (i.e., the first Krylov subspace).
 
Oh wow! That makes so much sense and seems so obvious. Thank you!