A deterministic Kaczmarz algorithm for solving linear systems

Abstract

We propose a new deterministic Kaczmarz algorithm for solving consistent linear systems Ax=b. Basically, the algorithm replaces orthogonal projections with reflections in the original scheme of Stefan Kaczmarz. Building on this, we give a geometric description of solutions of linear systems. Suppose A is m× n, we show that the algorithm generates a series of points distributed with patterns on an (n-1)-sphere centered on a solution. These points lie evenly on 2m lower-dimensional spheres \k0,k1\k=1m, with the property that for any k, the midpoint of the centers of k0,k1 is exactly a solution of Ax=b. With this discovery, we prove that taking the average of O(η(A)(1/)) points on any k0k1 effectively approximates a solution up to relative error , where η(A) characterizes the eigengap of the orthogonal matrix produced by the product of m reflections generated by the rows of A. We also analyze the connection between η(A) and (A), the condition number of A. In the worst case η(A)=O(2(A) m), while for random matrices η(A)=O((A)) on average. Finally, we prove that the algorithm indeed solves the linear system AT W-1A x = AT W-1 b, where W is the lower-triangular matrix such that W+WT = 2AAT. The connection between this linear system and the original one is studied. The numerical tests indicate that this new Kaczmarz algorithm has comparable performance to randomized (block) Kaczmarz algorithms.

0

Turn this paper into a lesson

ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…