{"title":"求解对称特征值问题的贪婪随机坐标更新迭代方法","authors":"Zhong-Zhi Bai","doi":"10.1016/j.apnum.2025.04.010","DOIUrl":null,"url":null,"abstract":"<div><div>In order to compute the smallest eigenvalue and its corresponding eigenvector of a large-scale, real, and symmetric matrix, we propose a class of greedy randomized coordinate updating iteration methods based on the principle that the indices of larger entries in absolute value of the current residual are selected with a higher probability and, with respect to this index set, the next iterate is updated from the current iterate along with the selected coordinate such that the corresponding entry of the residual is annihilated, resulting in fast convergence rates of the proposed iteration methods. Under appropriate conditions, we prove the convergence of both sequences of the Rayleigh-quotients and the acute angles between the iterates and the eigenvector in terms of the expectation. By numerical experiments, we show that this class of greedy randomized coordinate updating iteration methods are advantageous over the parameterized power method and the coordinate descent method in both iteration counts and computing times. Moreover, with theoretical analysis and computational performance, we confirm that the convergence property of this class of iteration methods can be improved significantly by suitably choosing the arbitrary nonnegative parameter involved in the greedy probability criterion.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"216 ","pages":"Pages 76-97"},"PeriodicalIF":2.2000,"publicationDate":"2025-05-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"On greedy randomized coordinate updating iteration methods for solving symmetric eigenvalue problems\",\"authors\":\"Zhong-Zhi Bai\",\"doi\":\"10.1016/j.apnum.2025.04.010\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"<div><div>In order to compute the smallest eigenvalue and its corresponding eigenvector of a large-scale, real, and symmetric matrix, we propose a class of greedy randomized coordinate updating iteration methods based on the principle that the indices of larger entries in absolute value of the current residual are selected with a higher probability and, with respect to this index set, the next iterate is updated from the current iterate along with the selected coordinate such that the corresponding entry of the residual is annihilated, resulting in fast convergence rates of the proposed iteration methods. Under appropriate conditions, we prove the convergence of both sequences of the Rayleigh-quotients and the acute angles between the iterates and the eigenvector in terms of the expectation. By numerical experiments, we show that this class of greedy randomized coordinate updating iteration methods are advantageous over the parameterized power method and the coordinate descent method in both iteration counts and computing times. Moreover, with theoretical analysis and computational performance, we confirm that the convergence property of this class of iteration methods can be improved significantly by suitably choosing the arbitrary nonnegative parameter involved in the greedy probability criterion.</div></div>\",\"PeriodicalId\":8199,\"journal\":{\"name\":\"Applied Numerical Mathematics\",\"volume\":\"216 \",\"pages\":\"Pages 76-97\"},\"PeriodicalIF\":2.2000,\"publicationDate\":\"2025-05-09\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"Applied Numerical Mathematics\",\"FirstCategoryId\":\"100\",\"ListUrlMain\":\"https://www.sciencedirect.com/science/article/pii/S0168927425000972\",\"RegionNum\":2,\"RegionCategory\":\"数学\",\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q1\",\"JCRName\":\"MATHEMATICS, APPLIED\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Applied Numerical Mathematics","FirstCategoryId":"100","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0168927425000972","RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"MATHEMATICS, APPLIED","Score":null,"Total":0}
On greedy randomized coordinate updating iteration methods for solving symmetric eigenvalue problems
In order to compute the smallest eigenvalue and its corresponding eigenvector of a large-scale, real, and symmetric matrix, we propose a class of greedy randomized coordinate updating iteration methods based on the principle that the indices of larger entries in absolute value of the current residual are selected with a higher probability and, with respect to this index set, the next iterate is updated from the current iterate along with the selected coordinate such that the corresponding entry of the residual is annihilated, resulting in fast convergence rates of the proposed iteration methods. Under appropriate conditions, we prove the convergence of both sequences of the Rayleigh-quotients and the acute angles between the iterates and the eigenvector in terms of the expectation. By numerical experiments, we show that this class of greedy randomized coordinate updating iteration methods are advantageous over the parameterized power method and the coordinate descent method in both iteration counts and computing times. Moreover, with theoretical analysis and computational performance, we confirm that the convergence property of this class of iteration methods can be improved significantly by suitably choosing the arbitrary nonnegative parameter involved in the greedy probability criterion.
期刊介绍:
The purpose of the journal is to provide a forum for the publication of high quality research and tutorial papers in computational mathematics. In addition to the traditional issues and problems in numerical analysis, the journal also publishes papers describing relevant applications in such fields as physics, fluid dynamics, engineering and other branches of applied science with a computational mathematics component. The journal strives to be flexible in the type of papers it publishes and their format. Equally desirable are:
(i) Full papers, which should be complete and relatively self-contained original contributions with an introduction that can be understood by the broad computational mathematics community. Both rigorous and heuristic styles are acceptable. Of particular interest are papers about new areas of research, in which other than strictly mathematical arguments may be important in establishing a basis for further developments.
(ii) Tutorial review papers, covering some of the important issues in Numerical Mathematics, Scientific Computing and their Applications. The journal will occasionally publish contributions which are larger than the usual format for regular papers.
(iii) Short notes, which present specific new results and techniques in a brief communication.