Modified Richardson iteration

{{Short description|Iterative method used to solve a linear system of equations}}

Modified Richardson iteration is an iterative method for solving a system of linear equations. Richardson iteration was proposed by Lewis Fry Richardson in his work dated 1910. 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.$

The Richardson iteration is

$x^{(k+1)} = x^{(k)} + \omega \left( b - A x^{(k)} \right),$

where $\omega$ is a scalar parameter that has to be chosen such that the sequence $x^{(k)}$ converges.

It is easy to see that the method has the correct fixed points, because if it converges, then $x^{(k+1)} \approx x^{(k)}$ and $x^{(k)}$ has to approximate a solution of $A x = b$ .

Convergence

Subtracting the exact solution $x$ , and introducing the notation for the error $e^{(k)} = x^{(k)} - x$ , we get the equality for the errors

$e^{(k+1)} = e^{(k)} - \omega A e^{(k)} = \left(I - \omega A\right) e^{(k)}.$

Thus,

$\left\|e^{(k+1)}\right\| = \left\|\left(I - \omega A\right) e^{(k)}\right\| \leq \left\|I - \omega A\right\| \left\|e^{(k)}\right\|,$

for any vector norm and the corresponding induced matrix norm. Thus, if $\|I-\omega A\| < 1$ , the method converges.

Suppose that $A$ is symmetric positive definite and that $(\lambda_j)_j$ are the eigenvalues of $A$ . The error converges to $0$ if $| 1 - \omega \lambda_j | < 1$ for all eigenvalues $\lambda_j$ . If, e.g., all eigenvalues are positive, this can be guaranteed if $\omega$ is chosen such that $0 < \omega < \omega_\text{max}\,, \ \omega_\text{max}:= 2/\lambda_{\text{max}}(A)$ . The optimal choice, minimizing all $| 1 - \omega \lambda_j |$ , is $\omega_\text{opt} := 2/(\lambda_\text{min}(A) + \lambda_\text{max}(A))$ , which gives the simplest Chebyshev iteration. This optimal choice yields a spectral radius of

$\min_{\omega\in (0,\omega_\text{max}) } \rho \left(I - \omega A\right)
= \rho \left(I - \omega_\text{opt} A\right)
= 1 - \frac{2}{\kappa(A)+1} \,,$

where $\kappa(A)$ is the condition number.

If there are both positive and negative eigenvalues, the method will diverge for any $\omega$ if the initial error $e^{(0)}$ has nonzero components in the corresponding eigenvectors.

Equivalence to [[gradient descent]]

Consider minimizing the function $F(x) = \frac{1}{2} \left\|\tilde{A}x - \tilde{b}\right\|_2^2$ . Since this is a convex function, a sufficient condition for optimality is that the gradient is zero ( $\nabla F(x) = 0$ ) which gives rise to the equation

$\tilde{A}^T \tilde{A} x = \tilde{A}^T\tilde{b}.$

Define $A = \tilde{A}^T\tilde{A}$ and $b = \tilde{A}^T\tilde{b}$ .

Because of the form of A, it is a positive semi-definite matrix, so it has no negative eigenvalues.

A step of gradient descent is

$x^{(k+1)} = x^{(k)} - t \nabla F(x^{(k)}) = x^{(k)} - t \left( Ax^{(k)} - b \right)$

which is equivalent to the Richardson iteration by making $t=\omega$ .

References

{{cite journal

| last = Richardson

| first = L.F.

| year = 1910

| journal = Philosophical Transactions of the Royal Society A

| title = The approximate arithmetical solution by finite differences of physical problems involving differential equations, with an application to the stresses in a masonry dam

| volume = 210

| issue = 459–470

| pages = 307–357

| jstor = 90994

| doi=10.1098/rsta.1911.0009

| doi-access =

}}

{{SpringerEOM

| title = Chebyshev iteration method

| first = Vyacheslav Ivanovich

| last = Lebedev

| author-link = Vyacheslav Ivanovich Lebedev

| accessdate= 2010-05-25

| editor = Michiel Hazewinkel

| editor-link = Michiel Hazewinkel

| isbn = 1-4020-0609-8

}}

Category:Numerical linear algebra

Category:Iterative methods

Modified Richardson iteration

Convergence

Equivalence to [[gradient descent]]

See also

References