Regularized randomized iterative algorithms for factorized linear systems

From MaRDI portal
Publication:6128669



Abstract: Randomized iterative algorithms for solving a factorized linear system, mathbfAmathbfBmathbfx=mathbfb with mathbfAinmathbbRmimesell, mathbfBinmathbbRellimesn, and mathbfbinmathbbRm, have recently been proposed. They take advantage of the factorized form and avoid forming the matrix mathbfC=mathbfAmathbfB explicitly. However, they can only find the minimum norm (least squares) solution. In contrast, the regularized randomized Kaczmarz (RRK) algorithm can find solutions with certain structures from consistent linear systems. In this work, by combining the randomized Kaczmarz algorithm or the randomized Gauss--Seidel algorithm with the RRK algorithm, we propose two novel regularized randomized iterative algorithms to find (least squares) solutions with certain structures of mathbfAmathbfBmathbfx=mathbfb. We prove linear convergence of the new algorithms. Computed examples are given to illustrate that the new algorithms can find sparse (least squares) solutions of mathbfAmathbfBmathbfx=mathbfb and can be better than the existing randomized iterative algorithms for the corresponding full linear system mathbfCmathbfx=mathbfb with mathbfC=mathbfAmathbfB.












This page was built for publication: Regularized randomized iterative algorithms for factorized linear systems

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6128669)