{"title":"对Rn×n中HSS迭代方法的一个wachpress - habetler扩展","authors":"Thomas Smotzer, John Buoni","doi":"10.1016/j.apnum.2025.04.002","DOIUrl":null,"url":null,"abstract":"<div><div>In the study of implicit iterations for the two dimensional heat and Helmhotz equations, one constructs a splitting of the form <span><math><mi>A</mi><mo>=</mo><mi>U</mi><mo>+</mo><mi>V</mi></math></span> where <em>U</em> and <em>V</em> are the difference approximations parallel to the <em>x</em> and <em>y</em> axis, respectively. In the commutative case for <em>U</em> and <em>V</em>, several investigations have taken place. For the noncommutative case, a symmetric positive definite matrix <em>F</em> is found, such that <span><math><mi>U</mi><mi>F</mi><mi>V</mi><mo>=</mo><mi>V</mi><mi>F</mi><mi>U</mi></math></span> and then, the investigations use this to address the non-commutative nature of <em>U</em> with <em>V</em>. The purpose of this paper is to study the same problem type for the commutativity of <em>H</em> and <em>S</em> in the <span><math><mi>H</mi><mi>S</mi><mi>S</mi></math></span> splitting of <span><math><mi>A</mi><mo>=</mo><mi>H</mi><mo>+</mo><mi>S</mi></math></span>, where <em>H</em> and <em>S</em> are the symmetric and skew-symmetric parts of <em>A</em>, respectively. Although throughout the literature in <span><math><msup><mrow><mi>C</mi></mrow><mrow><mi>n</mi><mo>×</mo><mi>n</mi></mrow></msup></math></span>, the <em>H</em> of <span><math><mi>H</mi><mi>S</mi><mi>S</mi></math></span> stands for hermitian, we use it in the symmetric matrix case. One then applies this result to the <span><math><mi>H</mi><mi>S</mi><mi>S</mi></math></span> iteration method. Since it is well known that the commutativity of <em>U</em> and <em>V</em> plays an important role in the analysis of <em>ADI</em> methods; especially for the solution of the Helmhotz Equation, one hopes that this commutativity will improve performance of the <span><math><mi>H</mi><mi>S</mi><mi>S</mi></math></span> method. This extension is similar to that of Wachspress and Habetler's variation of the Peaceman-Rachford method. Along the way, suitable conditions are found for <em>A</em>, which yield a symmetric non-zero matrix <span><math><mi>P</mi><mo>=</mo><msqrt><mrow><mi>F</mi></mrow></msqrt></math></span> for which <span><math><mi>N</mi><mspace></mspace><mo>=</mo><mspace></mspace><mi>P</mi><mi>A</mi><mi>P</mi></math></span> is a normal matrix.</div></div>","PeriodicalId":8199,"journal":{"name":"Applied Numerical Mathematics","volume":"215 ","pages":"Pages 49-58"},"PeriodicalIF":2.2000,"publicationDate":"2025-04-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"A Wachspress-Habetler extension to the HSS iteration method in Rn×n\",\"authors\":\"Thomas Smotzer, John Buoni\",\"doi\":\"10.1016/j.apnum.2025.04.002\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"<div><div>In the study of implicit iterations for the two dimensional heat and Helmhotz equations, one constructs a splitting of the form <span><math><mi>A</mi><mo>=</mo><mi>U</mi><mo>+</mo><mi>V</mi></math></span> where <em>U</em> and <em>V</em> are the difference approximations parallel to the <em>x</em> and <em>y</em> axis, respectively. In the commutative case for <em>U</em> and <em>V</em>, several investigations have taken place. For the noncommutative case, a symmetric positive definite matrix <em>F</em> is found, such that <span><math><mi>U</mi><mi>F</mi><mi>V</mi><mo>=</mo><mi>V</mi><mi>F</mi><mi>U</mi></math></span> and then, the investigations use this to address the non-commutative nature of <em>U</em> with <em>V</em>. The purpose of this paper is to study the same problem type for the commutativity of <em>H</em> and <em>S</em> in the <span><math><mi>H</mi><mi>S</mi><mi>S</mi></math></span> splitting of <span><math><mi>A</mi><mo>=</mo><mi>H</mi><mo>+</mo><mi>S</mi></math></span>, where <em>H</em> and <em>S</em> are the symmetric and skew-symmetric parts of <em>A</em>, respectively. Although throughout the literature in <span><math><msup><mrow><mi>C</mi></mrow><mrow><mi>n</mi><mo>×</mo><mi>n</mi></mrow></msup></math></span>, the <em>H</em> of <span><math><mi>H</mi><mi>S</mi><mi>S</mi></math></span> stands for hermitian, we use it in the symmetric matrix case. One then applies this result to the <span><math><mi>H</mi><mi>S</mi><mi>S</mi></math></span> iteration method. Since it is well known that the commutativity of <em>U</em> and <em>V</em> plays an important role in the analysis of <em>ADI</em> methods; especially for the solution of the Helmhotz Equation, one hopes that this commutativity will improve performance of the <span><math><mi>H</mi><mi>S</mi><mi>S</mi></math></span> method. This extension is similar to that of Wachspress and Habetler's variation of the Peaceman-Rachford method. Along the way, suitable conditions are found for <em>A</em>, which yield a symmetric non-zero matrix <span><math><mi>P</mi><mo>=</mo><msqrt><mrow><mi>F</mi></mrow></msqrt></math></span> for which <span><math><mi>N</mi><mspace></mspace><mo>=</mo><mspace></mspace><mi>P</mi><mi>A</mi><mi>P</mi></math></span> is a normal matrix.</div></div>\",\"PeriodicalId\":8199,\"journal\":{\"name\":\"Applied Numerical Mathematics\",\"volume\":\"215 \",\"pages\":\"Pages 49-58\"},\"PeriodicalIF\":2.2000,\"publicationDate\":\"2025-04-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/S0168927425000819\",\"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/S0168927425000819","RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"MATHEMATICS, APPLIED","Score":null,"Total":0}
A Wachspress-Habetler extension to the HSS iteration method in Rn×n
In the study of implicit iterations for the two dimensional heat and Helmhotz equations, one constructs a splitting of the form where U and V are the difference approximations parallel to the x and y axis, respectively. In the commutative case for U and V, several investigations have taken place. For the noncommutative case, a symmetric positive definite matrix F is found, such that and then, the investigations use this to address the non-commutative nature of U with V. The purpose of this paper is to study the same problem type for the commutativity of H and S in the splitting of , where H and S are the symmetric and skew-symmetric parts of A, respectively. Although throughout the literature in , the H of stands for hermitian, we use it in the symmetric matrix case. One then applies this result to the iteration method. Since it is well known that the commutativity of U and V plays an important role in the analysis of ADI methods; especially for the solution of the Helmhotz Equation, one hopes that this commutativity will improve performance of the method. This extension is similar to that of Wachspress and Habetler's variation of the Peaceman-Rachford method. Along the way, suitable conditions are found for A, which yield a symmetric non-zero matrix for which is a normal matrix.
期刊介绍:
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.