MCSOR

Multicolour Gauss-Seidel/SOR.

Notes

This implements a true parallel Gauss-Seidel method (as opposed to PETSc’s parallel SOR which is actually block Jacobi with Gauss-Seidel in the blocks.

Implemented for MATAIJ and MATLRC matrices (with MATAIJ as the base matrix type).

Developer notes

Should this be a PC?