Modified Richardson iteration is an iterative method for solving a system of linear equations. Richardson iteration was proposed by Lewis Fry Richardson in work that was published in 1911. It is similar to the Jacobi and Gauss–Seidel method. We seek the solution to a set of linear equations, expressed in matrix terms as
A x = b . {\displaystyle Ax=b.}
The Richardson iteration is
x ( k + 1 ) = x ( k ) + ω ( b − A x ( k ) ) , {\displaystyle x^{(k+1)}=x^{(k)}+\omega \left(b-Ax^{(k)}\right),}
where ω {\displaystyle \omega } is a scalar parameter that has to be chosen such that the sequence x ( k ) {\displaystyle x^{(k)}} converges. It is easy to see that the method has the correct fixed points, because if it converges, then x ( k + 1 ) ≈ x ( k ) {\displaystyle x^{(k+1)}\approx x^{(k)}} and x ( k ) {\displaystyle x^{(k)}} has to approximate a solution of A x = b {\displaystyle Ax=b} .
Convergence Subtracting the exact solution x {\displaystyle x} , and introducing the notation for the error e ( k ) = x ( k ) − x {\displaystyle e^{(k)}=x^{(k)}-x} , we get the equality for the errors
e ( k + 1 ) = e ( k ) − ω A e ( k ) = ( I − ω A ) e ( k ) . {\displaystyle e^{(k+1)}=e^{(k)}-\omega Ae^{(k)}=\left(I-\omega A\right)e^{(k)}.}
Thus,
‖ e ( k + 1 ) ‖ = ‖ ( I − ω A ) e ( k ) ‖ ≤ ‖ I − ω A ‖ ‖ e ( k ) ‖ , {\displaystyle \left\|e^{(k+1)}\right\|=\left\|\left(I-\omega A\right)e^{(k)}\right\|\leq \left\|I-\omega A\right\|\left\|e^{(k)}\right\|,}
… excerpt ends here. Continue reading the full article.
