The Jacobi Method

next up previous contents index
Next: Convergence of the Up: Stationary Iterative Methods Previous: Stationary Iterative Methods

The Jacobi Method


The Jacobi method is easily derived by examining each of the equations in the linear system in isolation. If in the th equation

we solve for the value of while assuming the other entries of remain fixed, we obtain


This suggests an iterative method defined by


which is the Jacobi method. Note that the order in which the equations are examined is irrelevant, since the Jacobi method treats them independently. For this reason, the Jacobi method is also known as the method of simultaneous displacements, since the updates could in principle be done simultaneously.

Simultaneous displacements, method of: Jacobi method.

In matrix terms, the definition of the Jacobi method in (gif) can be expressed as


where the matrices , and represent the diagonal, the strictly lower-triangular, and the strictly upper-triangular parts of , respectively.

The pseudocode for the Jacobi method is given in Figure gif. Note that an auxiliary storage vector, is used in the algorithm. It is not possible to update the vector in place, since values from are needed throughout the computation of .

Figure: The Jacobi Method

Jack Dongarra
Mon Nov 20 08:52:54 EST 1995