Abstract:The Kaczmarz algorithm is a simple iterative scheme for solving consistent linear systems. At each step, the method projects the current iterate onto the solution space of a single constraint. Hence, it requires low cost per iteration and storage, and it has a linear rate of convergence. Distributed implementations of Kaczmarz have recently become the de facto architectural choice for large-scale linear systems. Therefore, in this paper we develop a family of randomized block Kaczmarz algorithms that uses at each step a subset of the constraints and extrapolated stepsizes, and can be deployed on distributed computing units. Our approach is based on several new ideas and tools, including stochastic selection rules for the blocks of rows, stochastic conditioning of linear systems, and novel strategies for designing extrapolated stepsizes. We prove that randomized block Kaczmarz algorithms converge linearly in expectation, with a rate depending on the geometric properties of the matrix and its submatrices and on the size of the blocks. Our convergence analysis reveals that the algorithm is most effective when it is given a good sampling of the rows into well-conditioned blocks. Besides providing a general framework for the design and analysis of randomized block Kaczmarz methods, our results resolve an open problem in the literature related to the theoretical understanding of observed practical efficiency of extrapolated block Kaczmarz methods. We also propose an accelerated block Kaczmarz scheme, that is, acceleration in the sense of Chebyshev semi-iterative methods, where the stepsize is chosen based on the roots of Chebyshev polynomials, and we derive convergence rates depending on the square root of the geometric properties of the matrix. Finally, numerical examples illustrate the benefits of the new algorithms.

Randomized Kaczmarz with tail averaging

An accelerated randomized Kaczmarz algorithm

Faster randomized block sparse Kaczmarz by averaging

A Kaczmarz Method with Simple Random Sampling for Solving Large Linear Systems

Randomized Kaczmarz Algorithm for Inconsistent Linear Systems: An Exact MSE Analysis

Randomized Kaczmarz algorithm with averaging and block projection

Faster Randomized Block Kaczmarz Algorithms

Randomized Kaczmarz method with adaptive stepsizes for inconsistent linear systems

Randomized Extended Kaczmarz for Solving Least-Squares

On the Randomized Kaczmarz Algorithm

Randomized block Kaczmarz method with projection for solving least squares

The randomized Kaczmarz algorithm with the probability distribution depending on the angle

Parallelization Strategies for the Randomized Kaczmarz Algorithm on Large-Scale Dense Systems

An Asynchronous Parallel Randomized Kaczmarz Algorithm

A Note on the Randomized Kaczmarz Algorithm for Solving Doubly Noisy Linear Systems

Paved with good intentions: Analysis of a randomized block Kaczmarz method

A Note on Randomized Kaczmarz Algorithm for Solving Doubly-Noisy Linear Systems

On the alternating randomized block Kaczmarz method

Rates of convergence of randomized Kaczmarz algorithms in Hilbert spaces

Rows vs. Columns: Randomized Kaczmarz or Gauss-Seidel for Ridge Regression

A Randomized Distributed Kaczmarz Algorithm and Anomaly Detection