{"title":"一类线性系统的Orthomin(k)收敛速率的尖锐估计","authors":"Andrei Drăgănescu , Florin Spinu","doi":"10.1016/j.cam.2025.116699","DOIUrl":null,"url":null,"abstract":"<div><div>In this work we show that the convergence rate of Orthomin(<span><math><mi>k</mi></math></span>) applied to systems of the form <span><math><mrow><mrow><mo>(</mo><mi>I</mi><mo>+</mo><mi>ρ</mi><mi>U</mi><mo>)</mo></mrow><mi>x</mi><mo>=</mo><mi>b</mi></mrow></math></span>, where <span><math><mi>U</mi></math></span> is a unitary operator and <span><math><mrow><mn>0</mn><mo><</mo><mi>ρ</mi><mo><</mo><mn>1</mn></mrow></math></span>, is less than or equal to <span><math><mi>ρ</mi></math></span>. Moreover, we give examples of operators <span><math><mi>U</mi></math></span> and <span><math><mrow><mi>ρ</mi><mo>></mo><mn>0</mn></mrow></math></span> for which the asymptotic convergence rate of Orthomin(<span><math><mi>k</mi></math></span>) is exactly <span><math><mi>ρ</mi></math></span>, thus showing that the estimate is sharp. While the systems under scrutiny may not be of great interest in themselves, their existence shows that, in general, Orthomin(<span><math><mi>k</mi></math></span>) does not converge faster than Orthomin(1). Furthermore, we give examples of systems for which Orthomin(<span><math><mi>k</mi></math></span>) has the same asymptotic convergence rate as Orthomin(2) for <span><math><mrow><mi>k</mi><mo>≥</mo><mn>2</mn></mrow></math></span>, but smaller than that of Orthomin(1). The latter systems are related to the numerical solution of certain partial differential equations.</div></div>","PeriodicalId":50226,"journal":{"name":"Journal of Computational and Applied Mathematics","volume":"470 ","pages":"Article 116699"},"PeriodicalIF":2.1000,"publicationDate":"2025-04-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Sharp estimates for the convergence rate of Orthomin(k) for a class of linear systems\",\"authors\":\"Andrei Drăgănescu , Florin Spinu\",\"doi\":\"10.1016/j.cam.2025.116699\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"<div><div>In this work we show that the convergence rate of Orthomin(<span><math><mi>k</mi></math></span>) applied to systems of the form <span><math><mrow><mrow><mo>(</mo><mi>I</mi><mo>+</mo><mi>ρ</mi><mi>U</mi><mo>)</mo></mrow><mi>x</mi><mo>=</mo><mi>b</mi></mrow></math></span>, where <span><math><mi>U</mi></math></span> is a unitary operator and <span><math><mrow><mn>0</mn><mo><</mo><mi>ρ</mi><mo><</mo><mn>1</mn></mrow></math></span>, is less than or equal to <span><math><mi>ρ</mi></math></span>. Moreover, we give examples of operators <span><math><mi>U</mi></math></span> and <span><math><mrow><mi>ρ</mi><mo>></mo><mn>0</mn></mrow></math></span> for which the asymptotic convergence rate of Orthomin(<span><math><mi>k</mi></math></span>) is exactly <span><math><mi>ρ</mi></math></span>, thus showing that the estimate is sharp. While the systems under scrutiny may not be of great interest in themselves, their existence shows that, in general, Orthomin(<span><math><mi>k</mi></math></span>) does not converge faster than Orthomin(1). Furthermore, we give examples of systems for which Orthomin(<span><math><mi>k</mi></math></span>) has the same asymptotic convergence rate as Orthomin(2) for <span><math><mrow><mi>k</mi><mo>≥</mo><mn>2</mn></mrow></math></span>, but smaller than that of Orthomin(1). The latter systems are related to the numerical solution of certain partial differential equations.</div></div>\",\"PeriodicalId\":50226,\"journal\":{\"name\":\"Journal of Computational and Applied Mathematics\",\"volume\":\"470 \",\"pages\":\"Article 116699\"},\"PeriodicalIF\":2.1000,\"publicationDate\":\"2025-04-22\",\"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/S0377042725002134\",\"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":"Journal of Computational and Applied Mathematics","FirstCategoryId":"100","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0377042725002134","RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"MATHEMATICS, APPLIED","Score":null,"Total":0}
Sharp estimates for the convergence rate of Orthomin(k) for a class of linear systems
In this work we show that the convergence rate of Orthomin() applied to systems of the form , where is a unitary operator and , is less than or equal to . Moreover, we give examples of operators and for which the asymptotic convergence rate of Orthomin() is exactly , thus showing that the estimate is sharp. While the systems under scrutiny may not be of great interest in themselves, their existence shows that, in general, Orthomin() does not converge faster than Orthomin(1). Furthermore, we give examples of systems for which Orthomin() has the same asymptotic convergence rate as Orthomin(2) for , but smaller than that of Orthomin(1). The latter systems are related to the numerical solution of certain partial differential equations.
期刊介绍:
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.