RICHARDSON
Usage:
-ksp_type richardson
Options:
-ksp_richardson_scale
damping factor on the correction (defaults to 1.0) where xn+1 = xn + scale*B(b - A xn)
This method often (usually) will not converge unless scale is very small
Uses the matrix splitting A = I – (I – A).
Characterized by the iteration rk+1 = Pk+1(A) r0, where Pk+1(A) = (I – A)k+1.