×

DQGMRES:基于不完全正交化的直接准最小残差算法。 (英语) Zbl 0906.65033号

提出了求解大型稀疏线性代数方程组的截断GMRES方法。作者提出了基于不完全Arnoldi正交化过程的DQGMRES算法,并计算了一系列具有拟最小残差性质的近似解。对新算法进行了理论研究,并在一些数值示例上进行了广泛测试。
在第一部分中,作者描述了不完全Arnoldi过程,其中基向量仅与少数先前向量正交。这种不完全正交化过程用于定义GMRES方法的一种变体,具有与其他截断Krylov方法相当的短期递归性和计算复杂性。证明了给定的DQGMRES算法在且仅当找到精确解(即发生所谓的快乐崩溃)时失效。结果表明,DQGMRES的残差范数与GMRES方法的最优残差范数之间的关系是基向量矩阵的条件数的一个因子。因此,保持由不完全Arnoldi过程计算的基向量的线性无关性对DQGMRES方法的收敛速度非常重要。DQGMRES可以进行预处理,与FGMRES方法类似,它是灵活的,即每个迭代步骤都可以改变预处理器。
论文的第二部分是数值实验。通过基于有限元的软件包FIDAP,对从Navier-Stokes方程解中提取的一组非对称矩阵进行了数值测试。考虑了求解非对称系统的几种迭代共轭梯度型方法:稳定双共轭梯度法(BiCGSTAB)、无转置拟最小残差法(TFQMR)、重启广义最小残差方法(GMRES((m))、不完全正交化方法的直接实现(DIOM((k))、,ORTHOMIN方法的截断版本和新提出的DQGMRES((k))方法。数值实验表明,DQGMRES((k))与其他截断Krylov子空间方法一样有效。采用了三种预处理技术。据报道,通过预处理,DQGMRES((k))的性能优于GMRES方法(2k),但通常比其他截断方法更有效。作者还考虑了内部-外部方案,并报告了选择DQGMRES作为外部迭代方法将导致与在外部迭代过程中使用重新启动的GMRES方法获得的结果大致相同的结果。

MSC公司:

65层10 线性系统的迭代数值方法
65层25 数值线性代数中的正交化
65层35 矩阵范数、条件、缩放的数值计算
65层50 稀疏矩阵的计算方法
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] 夸脱州阿莫尔迪。申请。静音。第17页第9页–(1951年)
[2] 阿克塞尔松,SLAMJ。材料分析。12 (1991)
[3] 布朗,SLAM J.数字。分析。第23页610–(1986)
[4] Eisenstat,SLAM J.数字。分析。第20页,第345页–(1983年)
[5] 法伯,SIAM J.Numer。分析。第21页,第352页–(1984年)
[6] Freund,《计算与应用数学杂志》,第34页,第135页–(1992年)
[7] 拟核多项式和拟最小剩余迭代的收敛结果。《近似理论的数值方法》和,编辑,Birkhauser,巴塞尔,1992年,第77-95页·Zbl 0814.65035号 ·doi:10.1007/978-3-0348-8619-2-5
[8] 纽纳·弗伦德。静音。第60页,315页–(1991年)
[9] SIAMJ弗伦德。科学。计算。第14页,第470页–(1993年)
[10] 在IGMRES(g)上,大型非对称线性系统的不完全广义最小残差方法。《技术报告94-047》,比勒费尔德大学数学系,桑德福斯琼斯贝里奇3431994年。上次修订日期:1995年3月。
[11] Lanczos算法的一个look-ahead变体及其在非厄米线性系统准最小残差法中的应用。麻省理工学院应用数学博士论文,马萨诸塞州剑桥,1991年。
[12] 萨阿德,SIAM J.Sci。统计师。计算。第7页,856页–(1986年)
[13] 和。WGMRES是一种基于不完全或正交化的准最小残差算法。技术报告UMSI-93/131,明尼苏达州超级计算研究所,明尼阿波利斯,1993年。
[14] 和。并行矩阵库(P-SPARSLIB)迭代求解模块的设计。《IUACS会议记录》,格鲁吉亚,1994年,编辑,1995年。曾任明尼苏达大学计算机科学系TR 94-59。
[15] Saad,数学。计算。第37页105–(1981)
[16] 萨阿德,SIAM J.Sci。统计师。计算。第5页203–(1984)
[17] SPARSKIT,一个用于稀疏矩阵计算的基本工具包。《技术报告90-20》,美国国家航空航天局艾姆斯研究中心高级计算机科学研究所,加利福尼亚州莫菲特菲尔德,1990年。
[18] 目前可从ftp//ftp.cs.umn.edu/dept/sparse/下载软件。
[19] 萨阿德,SIAM J.Sci。计算。第14页,461页–(1993年)
[20] 萨阿德,数值线性代数及其应用1第387页–(1994)
[21] 正交,一种求解稀疏线性方程组的迭代方法。第四届油藏模拟研讨会论文集,第149-159页。美国石油工程师学会,1976年。
[22] 范德沃斯特,《数值线性代数及其应用》1 pp 369–(1994)
[23] van der Vorst,SIAM科学杂志。统计师。计算。第13页,631页–(1992年)
[24] Walker,《数值线性代数及其应用》1 pp 571–(1994)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。