Added Gauss-Seidel relaxation routine