Knowledge base dedicated to Linux and applied mathematics.
All the versions of this article: [English] [français]
We will study an iterative method for solving linear systems: the Gauss-Seidel method. The aim is to build a sequence of approximations that converges to the true solution.
The Gauss-Seidel method is an iterative method for solving linear systems such as
where is an invertible matrix.
If is solution of then
Let be the error vector
We put , which gives
The algorithm converges if (null matrix).
Theorem: if and only if the spectral radius of the matrix
Theorem: If A is strictly diagonally dominant,
We decompose in the following way :
In the Gauss-Seidel method we choose and (in the Jacobi method, et ).
For the stop criteria , we can use the residual vector, wich gives for a given precision :