CFD Online Logo CFD Online URL
Home > Wiki > Jacobi method

Jacobi method

From CFD-Wiki

Jump to: navigation, search

We seek the solution to set of linear equations:

 A \cdot X = Q

For the given matrix A and vectors X and Q.
In matrix terms, the definition of the Jacobi method can be expressed as :
x^{(k)}  = D^{ - 1} \left( {L + U} \right)x^{(k - 1)}  + D^{ - 1} Q


Chose an intital guess X^{0} to the solution
for k := 1 step 1 untill convergence do
for i := 1 step until n do
 \sigma = 0
for j := 1 step until n do
if j != i then
end if
end (j-loop)
end (i-loop)
check if convergence is reached
end (k-loop)

Note: The major difference between the Gauss-Seidel method and Jacobi method lies in the fact that for Jacobi method the values of solution of previous iteration (here k) are used, where as in Gauss-Seidel method the latest available values of solution vector X are used.

Return to Numerical Methods

My wiki