×

QMR:非厄米线性系统的准最小残差法。 (英语) Zbl 0754.65034号

基本的Lanczos双正交方法[参见。C.兰佐斯《国家研究杂志》。伯尔。站立。49,33-53(1952;MR 14.501)]对于线性系统的解(Ax=b),(A)非Hermitian,生成序列(v_1,v_2,dots,v_n})和(w_1,w_2,dotes,w_n}1}=A^Tv_j-\alpha_j w_j-\ gamma_jw_{j-1})其中选择标量系数以满足双正交条件(w^T_kv_l=d_k\delta_{kl}\)。双共轭梯度(BCG)方法是Lanczos方法的变体。请注意,如果\(w^T_{n+1}v_{n+1}=0\),则必须终止上述过程,以防止在下一步中被零除。所谓的BCG前瞻性变体试图克服这一困难[cfB.N.帕莱特,D.R.泰勒,Z.A.刘,数学。计算。44, 105-124 (1985;Zbl 0564.65022号)].
本文提出了拟最小残差(QMR)方法,这是BCG的推广,克服了数值不稳定的趋势。它合并了look-ahead BCG的第(n)次迭代,从(v_1=r_0/\|r_0)开始,其中(r_0。将提供实现细节,以及进一步的属性和错误边界。
最后,给出了QMR和本文提到的其他迭代方法的大量数值实验结果。

MSC公司:

65层10 线性系统的迭代数值方法
65N22型 含偏微分方程边值问题离散方程的数值解
PDF格式BibTeX公司 XML格式引用
全文: 内政部 欧洲DML

参考文献:

[1] Duff,I.S.,Grimes,R.G.,Lewis,J.G.(1989):稀疏矩阵测试问题。ACM事务处理。数学。软15,1-14·Zbl 0667.65040号 ·数字对象标识代码:10.1145/62038.62043
[2] Faber,V.,Manteuffel,T.(1984):共轭梯度法存在的充要条件。SIAM J.数字。分析21,352-362·Zbl 0546.65010号 ·doi:10.1137/0721026
[3] Fischer,B.,Freund,R.W.(1990):关于椭圆上的约束Chebyshev逼近问题。J.近似理论62,297-315·Zbl 0728.41023号 ·doi:10.1016/0021-9045(90)90054-T
[4] Fletcher,R.(1976):不定系统的共轭梯度方法。收录于:G.A.Watson主编,《邓迪数值分析》,1975年,第73-89页。数学课堂讲稿506。施普林格,柏林-海德堡-纽约
[5] Freund,R.W.(1989):具有复对称系数矩阵的线性系统的共轭梯度型方法。技术报告89.54,RIACS,NASA艾姆斯研究中心
[6] Freund,R.W.,Gutknecht,M.H.,Nachtigal,N.M.(1990):非厄米矩阵前瞻Lanczos算法的实现,第部分。I.技术报告90.45,RIACS,NASA艾姆斯研究中心
[7] Freund,R.W.,Nachtigal,N.M.(1990):非厄米矩阵的look-ahead Lanczos算法的实现,第二部分。技术报告90.46,RIACS,NASA艾姆斯研究中心·Zbl 0770.65022号
[8] Golub,G.H.,Van Loan,C.F.(1983):矩阵计算。约翰·霍普金斯大学出版社,巴尔的摩·Zbl 0559.65011号
[9] Gutknecht,M.H.(1990):非对称Lanczos过程和相关算法的完整理论,第二部分。IPS研究报告第90-16号,Z?富有的
[10] Gutknecht,M.H.(1990):非对称Lanczos过程和相关算法的完整理论,第二部分。IPS研究报告第90-16号,Z?富有的
[11] Hestenes,M.R.,Stiefel,E.(1952):求解线性系统的共轭梯度方法。《国家研究杂志》。伯尔。展位49,409-436·Zbl 0048.09901号
[12] Lanczos,C.(1950):求解线性微分和积分算子特征值问题的迭代方法。《国家研究杂志》。伯尔。展位45,255-282
[13] Lanczos,C.(1952):通过最小化迭代求解线性方程组。《国家研究杂志》。伯尔。展位49,33-53
[14] Manteuffel,T.A.(1977):非对称线性系统的切比雪夫迭代。数字。数学28,307-327·Zbl 0361.65024号 ·doi:10.1007/BF01389971
[15] Meijerink,J.A.,van der Vorst,H.A.(1977):系数矩阵为对称M-矩阵的线性系统的迭代解。数学。组件31,148-162·Zbl 0349.65020号
[16] Paige,C.C.,Saunders,M.A.(1975):线性方程组稀疏不定系统的求解。SIAM J.数字。分析12,617-629·Zbl 0319.65025号 ·数字对象标识代码:10.1137/0712047
[17] Parlett,B.N.(1990):简化为三对角形式和最小实现。预印本,伯克利·Zbl 0754.65040号
[18] Parlett,B.N.,Taylor,D.R.,Liu,Z.A.(1985):非对称矩阵的look-ahead Lanczos算法。数学。组件44,105-124·Zbl 0564.65022号
[19] Saad,Y.(1982):求解大型非对称系统的Lanczos双正交算法和其他斜投影方法。SIAM J.数字。分析19,485-506·Zbl 0483.65022号 ·doi:10.1137/0719031
[20] Saad,Y.(1990):SPARSKIT:用于稀疏矩阵计算的基本工具包。技术报告90.20,RIACS,NASA艾姆斯研究中心
[21] Saad,Y.,Schultz,M.H.(1986):GMRES是求解非对称线性系统的广义最小残差算法。SIAM J.科学。统计计算7,856-869·Zbl 0599.65018号 ·doi:10.1137/0907058
[22] Sonneveld,P.(1989):CGS,非对称线性系统的快速Lanczos型解算器。SIAM J.科学。统计计算10,36-52·Zbl 0666.65029号 ·数字对象标识代码:10.1137/0910004
[23] Taylor,D.R.(1982):前瞻性Lanczos算法分析。加州大学伯克利分校博士论文
[24] van der Vorst,H.A.(1990):Bi-CGSTAB:用于非对称线性系统解的Bi-CG的一种快速且平滑收敛的变体。预印本,乌得勒支·Zbl 0761.65023号
[25] Wilkinson,J.H.(1965):代数特征值问题。牛津大学出版社·Zbl 0258.65037号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。