An implementation of SSOR (symmetric successive over relaxation, symmetric Gauss-Seidel) preconditioning that incorporates Eisenstat's trick to reduce the amount of computation needed.
Only implemented for the SeqAIJ matrix format.
Not a true parallel SOR, in parallel this implementation corresponds to block
Jacobi with SOR on each block.