Stabilization of Stochastic Iterative Methods for Singular and Nearly Singular Linear Systems

Abstract

We consider linear systems of equations, Ax = b, with an emphasis on the case where A is singular. Under certain conditions, necessary as well as sufficient, linear deterministic iterative methods generate sequences {x[subscript k]} that converge to a solution as long as there exists at least one solution. This convergence property can be impaired when these methods are implemented with stochastic simulation, as is often done in important classes of large-scale problems. We introduce additional conditions and novel algorithmic stabilization schemes under which {x[subscript k]} converges to a solution when A is singular and may also be used with substantial benefit when A is nearly singular.