There is an example program in your class directory: examples/pgi/laplace.hpf which solves this problem by Jacobi iteration. You may extend this program to a red/black SOR solver. (Compiler note: some of us have had difficulty with spreads and modules - try to avoid these constructs for now.)
***************************************************************************
By comparing the numbers of iteration, using red/black SOR technique is better than just using point-Jacobi method and it just take almost as half amount of iteration and less CPU time as using point-Jacobi method.
By observing the result, we know that the optimum value of parameter w for this problem is between 1.0 to 1.5, and then by comparing the converge behavior of eps value, we can estimate the optimum value of the over-relaxation parameter w probably equal to 1.3.
Ho-Yen Chang