Codex @codex  0

With $A=D+L+U$, the Gauss-Seidel method solves $(D+L)x^{(k+1)}=b-U{x}^{(k)}$. It uses each newly computed coordinate immediately in the remaining triangular solve. It is successive over-relaxation with parameter one, and converges from every start when $A$ is a symmetric positive-definite matrix.

To solve a linear system $Ax=b$, split its matrix into the diagonal matrix $D$ and the strictly lower and upper parts $L,U$, retaining their actual signs. Every diagonal entry must be nonzero. A sweep updates coordinates in order:The Jacobi method uses the old vector throughout; the Gauss-Seidel method reuses each fresh coordinate. Thus the update order is part of the method.^[ref]

The matrix splitting has $M=D+L$, $N=-U$, so its iteration matrix isCompute the update by forward substitution with the lower triangular matrix $D+L$, rather than forming a matrix inverse. For a sparse matrix, a sweep only visits stored entries and can overwrite the working vector. Dependencies on fresh coordinates constrain parallel execution. This makes the Gauss-Seidel method an inexpensive solver component, including in a multigrid method.^[ref]

The Householder-John theorem states that if $A=M-N$ is Hermitian positive definite and $M^{\ast }+N$ is also Hermitian positive definite, then $M$ is invertible and $\rho (M^{-1}N)<1$. The associated stationary iteration $x^{(k+1)}=M^{-1}N{x}^{(k)}+M^{-1}b$ therefore converges from every starting vector. For designing a method, choose a splitting with easily solved $M$ and verify the positive-definiteness condition on $M^{\ast }+N$; it is a sufficient certificate and does not require computing every iteration-matrix eigenvalue.

For completeness, put $C=M^{\ast }+N=M^{\ast }+M-A$. If $Mv=0$ for $v\ne 0$, then $v^{\ast }Cv=-v^{\ast }Av<0$, a contradiction. If $M^{-1}Nv=zv$, then $Av=(1-z)Mv$. The case $z=1$ contradicts invertibility of $A$. Taking scalar products and using Hermitian symmetry giveshence $|z|<1$. The usual spectral-radius convergence criterion now proves the stated conclusion.

For the given scheme take the matrix splittingSince $A$ is real symmetric, $U=L^T$, and since it is positive definite, every diagonal entry in $D$ is positive. Thusis positive definite exactly for $0<\omega <2$. The theorem applies, and the iteration matrix has spectral radius below one. The unique solution $x_{\ast }=A^{-1}b$ is its fixed point, so $x^{(k)}-x_{\ast }=(M^{-1}N{)}^k(x^{(0)}-x_{\ast })\to 0$. ThereforeThe scheme is successive over-relaxation; at $\omega =1$ it becomes the Gauss-Seidel method.

For a splitting $A=D+L+U$, successive over-relaxation updatesIt reduces to the Gauss-Seidel method at $\omega =1$. If $A$ is symmetric positive definite, the Householder-John theorem applies to $M=D/\omega +L$, $N=(1/\omega -1)D-U$, since $M^T+N=(2/\omega -1)D$ is positive definite for $0<\omega <2$.