Fourier Analysis of the SOR Iteration

Abstract
The SOR iteration for solving linear systems of equations depends upon an overrelaxation factor ω. We show that, for the standard model problem of Poisson's equation on a rectangle, the optimal ω and corresponding convergence rate can be obtained rigorously by Fourier analysis. The trick is to tilt the space-time grid so that the SOR stencil becomes symmetrical. The tilted grid also gives new insight into the relationships between the Gauss-Seidel and Jacobi iterations and between the lexicographic and red-black orderings, and into the modified equation analysis of Garabedian.

This publication has 0 references indexed in Scilit: