Solving linear systems by methods based on a probabilistic interpretation

Abstract

In this paper it is demonstrated how the probabilistic concept of a stopping time in a random process may be used to generate an iterative method for solving a system of linear equations. Actually all known iterative approximation methods for solving linear equations are generated by various choices of a stopping time e. g. the point and block Jacobi methods, the point and block Gauss-Seidel Methods and overrelaxation methods are covered.

The probabilistic approach offers—in a natural way—the possibility of adapting the solution technique to the special structure of the problem. Moreover, posterior bounds for the solution are constructed, which lead to faster convergence of the approximations than with usual prior bounds.

Zusammenfassung

Es wird gezeigt, wie das probabilistische Konzept der Stoppzeit eines Zufallsprozesses für eine iterative Methode zur Lösung eines Systems linearer Gleichungen herangezogen werden kann. Alle bekannten iterativen Näherungsmethoden zur Lösung linearer Gleichungssysteme, wie z. B. die Blockmethoden nach Jakobi und Gauss-Seidel und Überrelaxationsmethoden, entsprechen verschieden gewählten Stoppzeiten.

Das probabilistische Verfahren bietet die Möglichkeit, Lösungstechniken für die spezielle Struktur des Problems zu adaptieren. Darüber hinaus führen die vorgeschlagenen Methoden zu einer schnelleren Konvergenz.