Convergency Gauss Jacobi Method
Convergency Gauss Jacobi Method
unknown x:
where:
Convergency
The standard convergence condition (for any iterative method) is when the spectral radius of the iteration matrix is less than 1:
The method is guaranteed to converge if the matrix A is strictly or irreducibly diagonally dominant. Strict row diagonal dominance means that for each row, the absolute value of the diagonal term is greater than the sum of absolute values of other terms:
The Jacobi method sometimes converges even if these conditions are not satisfied.
where:
Then A can be decomposed into a lower triangular component , and a strictly upper triangular component U:
The GaussSeidel method is an iterative technique that solves the left hand side of this expression for x, using previous value for x on the right hand side. Analytically, this may be written as:
However, by taking advantage of the triangular form of , the elements of x(k+1) can be computed sequentially using forward substitution:
Convergence The convergence properties of the GaussSeidel method are dependent on the matrix A. Namely, the procedure is known to converge if either: A is symmetric positive-definite, or A is strictly or irreducibly diagonally dominant. The GaussSeidel method sometimes converges even if these conditions are not satisfied.