Longze Tan , Xueping Guo , Mingyu Deng , Jingrun Chen
{"title":"On the adaptive deterministic block Kaczmarz method with momentum for solving large-scale consistent linear systems","authors":"Longze Tan , Xueping Guo , Mingyu Deng , Jingrun Chen","doi":"10.1016/j.cam.2024.116328","DOIUrl":null,"url":null,"abstract":"<div><div>The Kaczmarz method is a traditional and widely used iterative algorithm for solving large-scale consistent linear systems, while its improved block Kaczmarz-type methods have received much attention and research in recent years due to their excellent numerical performance. Hence, in this paper, we present a deterministic block Kaczmarz method with momentum, which is based on Polyak’s heavy ball method and a row selection criterion for a set of block-controlled indices defined by the Euclidean norm of the residual vector. The proposed method does not need to compute the pseudo-inverses of a row submatrix at each iteration and it adaptively selects and updates the set of block control indices, thus this is different from the block Kaczmarz-type methods that are based on projection and pre-partitioning of row indices. The theoretical analysis of the proposed method shows that it converges linearly to the unique least-norm solutions of the consistent linear systems. Numerical experiments demonstrate that the deterministic block Kaczmarz method with momentum method is more efficient than the existing block Kaczmarz-type methods.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"457 ","pages":"Article 116328"},"PeriodicalIF":2.1000,"publicationDate":"2024-10-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Computational and Applied Mathematics","FirstCategoryId":"100","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0377042724005764","RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"MATHEMATICS, APPLIED","Score":null,"Total":0}
引用次数: 0
Abstract
The Kaczmarz method is a traditional and widely used iterative algorithm for solving large-scale consistent linear systems, while its improved block Kaczmarz-type methods have received much attention and research in recent years due to their excellent numerical performance. Hence, in this paper, we present a deterministic block Kaczmarz method with momentum, which is based on Polyak’s heavy ball method and a row selection criterion for a set of block-controlled indices defined by the Euclidean norm of the residual vector. The proposed method does not need to compute the pseudo-inverses of a row submatrix at each iteration and it adaptively selects and updates the set of block control indices, thus this is different from the block Kaczmarz-type methods that are based on projection and pre-partitioning of row indices. The theoretical analysis of the proposed method shows that it converges linearly to the unique least-norm solutions of the consistent linear systems. Numerical experiments demonstrate that the deterministic block Kaczmarz method with momentum method is more efficient than the existing block Kaczmarz-type methods.
期刊介绍:
The Journal of Computational and Applied Mathematics publishes original papers of high scientific value in all areas of computational and applied mathematics. The main interest of the Journal is in papers that describe and analyze new computational techniques for solving scientific or engineering problems. Also the improved analysis, including the effectiveness and applicability, of existing methods and algorithms is of importance. The computational efficiency (e.g. the convergence, stability, accuracy, ...) should be proved and illustrated by nontrivial numerical examples. Papers describing only variants of existing methods, without adding significant new computational properties are not of interest.
The audience consists of: applied mathematicians, numerical analysts, computational scientists and engineers.