Compressed Row Storage (CRS) Successive Over-Relaxation (SOR) Method

Function uses Successive Over-Relaxation (SOR) method to solve a linear system in CRS format

This function solves a linear system Ax=b using the Successive Over-Relaxation iteration method. The A-matrix is supplied in Compressed Row Storage (CRS) format. The iteration is terminated when either (1) the user-specified maximum number of iterations has been reached or (2) the difference in the norm of successive iterates is less than a user-specified epsilon.

References:
(1) http://mathworld.wolfram.com/SuccessiveOverrelaxationMethod.html
(2) Yang, S., Gobbert, M., "The Optimal Relaxation Parameter for the SOR Method Applied to a Classical Model Problem", University of Maryland, Department of Mathematics and Statistics (2007).
(http://userpages.umbc.edu/~gobber/papers/YangGobbert2007SOR.pdf)