A Maximal Residual Two Subspace Projection Algorithm for Solving Least-Squares Problems

Abstract

We study a greedy coordinate descent method to solve large linear least-squares problems expanding on a randomized coordinate descent method presented by Leventhal and Lewis (Math. Oper. Res. 35: 641–654, 2010). For an overdetermined system, they proved its exponential convergence, regardless of its consistency. In our work, we study a greedy selection rule for the coordinate descent method which we refer to as the two-dimensional maximal residual Gauss-Seidel (D2MRGS) method. In this method, we select two coordinates in every iteration and treat the current approximation in those directions. Convergence is analyzed for the stated method and numerical experiments are provided to demonstrate its efficiency.

Keywords

Coordinate descent / Iterative techniques / Least-squares problem / Petrov-Galerkin condition / Projection methods / Randomized methods

References

[1] BaiZ-Z, WangL. On convergence rates of Kaczmarz-type methods with different selection rules of working rows. Appl. Numer. Math., 2023, 186: 289-319.

[2] BaiZ-Z, WangL. On multi-step randomized extended Kaczmarz method for solving large sparse inconsistent linear systems. Appl. Numer. Math., 2023, 192: 197-213.

[3] BaiZ-Z, WuW-T. On greedy randomized Kaczmarz method for solving large sparse linear systems. SIAM J. Sci. Comput., 2018, 40: A592-A606.

[4] BaiZ-Z, WuW-T. On greedy randomized coordinate descent methods for solving large linear least-squares problems. Numer. Linear Algebra Appl., 2019, 26. e2237

[5] BaiZ-Z, WuW-T. On greedy randomized augmented Kaczmarz method for solving large sparse inconsistent linear systems. SIAM J. Sci. Comput., 2021, 43: A3892-A3911.

[6] BaiZ-Z, WuW-T. Randomized Kaczmarz iteration methods: algorithmic extensions and convergence theory. Jpn. J. Ind. Appl. Math., 2023, 40: 1421-1443.

[7] DuK. Tight upper bounds for the convergence of the randomized extended Kaczmarz and Gauss-Seidel algorithms. Numer. Linear Algebra Appl., 2019, 26: 1-14.

[8] JinL-L, LiH-B. Greedy double subspaces coordinate descent method for solving linear least-squares problems. J. Comput Sci., 2023, 70. 102029

[9] KaczmarzS. Angenäherte auösung von systemen linearer gleichungen. Bull. Int. Acad. Polom. Sci. Lett. A, 1937, 35: 355-357

[10] KaczmarzS. Approximate solution of systems of linear equations. Int. J. Control, 1993, 57: 1269-1271.

[11] LeventhalD, LewisAS. Randomized methods for linear constraints: convergence rates and conditioning. Math. Oper. Res., 2010, 35: 641-654.

[12] LiaoY-M, LuT-X, YinF. A two-step randomized Gauss-Seidel method for solving large-scale linear least squares problems. Electron. Res. Arch., 2022, 30: 755-779.

[13] LiuY, JiangX-L, GuC-Q. On maximum residual block and two-step Gauss-Seidel algorithms for linear least-squares problems. Calcolo, 2021, 58: 13.

[14] LvZ, BaoW, LiW, WangF, WuG. On extended Kaczmarz methods with random sampling and maximum-distance for solving large inconsistent linear systems. Results Appl. Math., 2022, 13. 100240

[15] MaA, NeedellD, RamdasA. Convergence properties of the randomized extended Gauss-Seidel and Kaczmarz methods. SIAM J. Matrix Anal. Appl., 2015, 36: 1590-1604.

[16] MustafaA, SahaM. A generalized projection iterative method for solving non-singular linear systems. Math. Found. Comput., 2022, 5: 343-350.

[17] Mustafa, A., Saha, M.: A two dimensional randomized extended Gauss-Seidel algorithm for solving least squares problems. Numer. Algorithm 96, 665–686 (2023)

[18] MustafaA, SahaM. Maximal residual extended Kaczmarz and Gauss-Seidel methods-convergence properties and applications. J. Comput. Appl. Math., 2024, 448. 115920

[19] NeedellD, WardR. Two-subspace projection method for coherent overdetermined systems. J. Fourier Anal. Appl., 2013, 19: 256-269.

[20] NiuY-Q, ZhengB. A new randomized Gauss-Seidel method for solving linear least-squares problems. Appl. Math. Lett., 2021, 116. 107057

[21] Popa, C.: Convergence rates for Kaczmarz-type algorithms. Numer. Algorithm 79, 1–17 (2018)

[22] SaadYIterative Methods for Sparse Linear Systems, 20032PhiladelphiaSIAM.

[23] StrohmerT, VershyninR. A randomized Kaczmarz algorithm with exponential convergence. J. Fourier Anal. Appl., 2009, 15: 262-278.

[24] WatkinsDSFundamentals of Matrix Computations, 20042New YorkWiley

[25] Wu, W.: Paving the randomized Gauss-Seidel method. BSc Thesis, Scripps College, Claremont, California (2017)

[26] Wu, W.-T.: On two-subspace randomized extended Kaczmarz method for solving large linear least-squares problems. Numer. Algorithm 89, 1–31 (2022)

[27] Zouzias, A., Freris, N.M.: Randomized extended Kaczmarz for solving least squares. SIAM J. Matrix Anal. Appl. 34, 773–793 (2013)

RIGHTS & PERMISSIONS

Shanghai University