Randomized double and triple Kaczmarz for solving extended normal equations

Kui Du,Xiao-Hui Sun
DOI: https://doi.org/10.1007/s10092-021-00406-9
2021-03-19
CALCOLO
Abstract:The randomized Kaczmarz algorithm has received considerable attention recently because of its simplicity, speed, and the ability to approximately solve large-scale linear systems of equations. In this paper we propose randomized double and triple Kaczmarz algorithms to solve extended normal equations of the form <span class="mathjax-tex">\(\mathbf{A}^\top \mathbf{Ax}=\mathbf{A}^\top \mathbf{b}-\mathbf{c}\)</span>. The proposed algorithms avoid forming <span class="mathjax-tex">\(\mathbf{A}^\top \mathbf{A}\)</span> explicitly and work for <i>arbitrary</i> <span class="mathjax-tex">\(\mathbf{A}\in \mathbb {R}^{m\times n}\)</span> (full rank or rank-deficient, <span class="mathjax-tex">\(m\ge n\)</span> or <span class="mathjax-tex">\(m<n\)< span="">). <i>Tight</i> upper bounds showing exponential convergence in the mean square sense of the proposed algorithms are presented and numerical experiments are given to illustrate the theoretical results.</n\)<></span>
mathematics, applied
What problem does this paper attempt to address?