×

求解具有多个右手边的非对称线性系统的Lanczos型方法-矩阵和多项式解释。 (英语) Zbl 0954.65024号

本文从理论和数值两方面研究了一种新的求解具有多个右手边的线性(对称和非对称)系统的Lanczos型方法。该方法的目的是避免使用Lanczos方法的通常实现中会出现的多个Hankel矩阵,因为它们是由与每个初始残基相关的Krylov空间确定的。相反,这里考虑的Krylov空间是由两个任意的初始向量生成的,因此每次迭代都有一个Hankel矩阵来求解与每个右手边相关的每个系统。其结果之一是,与Lanczos方法不同,对每个解进行逐次逼近的矩阵多项式通常不是正交的。
本文首先详细研究了与Lanczos方法的基本实现相对应的两种实现,即Orthodir和Orthomin,但由于收敛性能差而放弃了它们,尽管它们在内存需求和每次迭代的计算成本方面具有良好的实现特性。Orthomin实现的无转置变体,以及基于H.A.van der Vorst的稳定双共轭梯度法(BiCGSTAB)[SIAM J.Sci.Stat.Compute.13,No.2,631-644(1992;Zbl 0761.65023号)]就收敛性进行了比较,在六个例子中,三个对称的和三个非对称的。
在对称情况下,在一个示例中,前者的行为明显好于BiCGSTAB实现,而在另外两种方法中的行为类似。对于非对称情况,情况就不那么简单了:在一个条件数较低的示例中,Orthomin实现的性能比BiCGSTAB好得多,但在其他两个示例中,后者明显优于前者,其中一个示例来自基准集,两个示例的条件数都较大。这并不奇怪,因为BiCGSTAB是为了处理数值不稳定的非对称问题而设计的。
然而,必须指出的是,论文中也提到的一般缺点是,残留物的总体表现较差,这表明需要进一步改进该方法。

MSC公司:

65层10 线性系统的迭代数值方法
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Beckermann,B.,形式正交多项式的稳定计算,Numer。算法,11,1-23(1996)·Zbl 0848.65007号
[2] Brezinski,C。;Redivo Zaglia,M.,《正交多项式计算中的故障》,数值。算法,1207-221(1991)·Zbl 0752.65011号
[3] Brezinski,C.,《生物正交性及其在数值分析中的应用》(1992),马塞尔·德克尔:马塞尔·戴克尔纽约·兹比尔0757.41001
[4] C.Brezinski,M.Redivo Zaglia,非对称线性系统的无转置Lanczos型算法,Numer。算法,待出现。;C.Brezinski,M.Redivo Zaglia,非对称线性系统的无转置Lanczos型算法,Numer。算法,出现·Zbl 0907.65031号
[5] Brezinski,C.,CGM:线性系统的一类Lanczos型解算器(1991年11月),里尔科技大学数值分析与优化实验室,注释ANO 253
[6] T.F.Chan,L.De Pillis,H.A.van der Vorst,非对称线性系统Lanczos型方法的无转位公式,Numer。算法,待出现。;T.F.Chan,L.De Pillis,H.A.van der Vorst,非对称线性系统的Lanczos型方法的无转座公式,Numer。算法,出现·兹比尔0907.65030
[7] Chan,T.F。;Wan,W.L.,《求解多右手边线性系统的投影方法分析》,SIAM J.Sci。计算。,18, 1698-1721 (1997) ·Zbl 0888.65033号
[8] Draux,A.,Polynómes Orthogonaux Formels。应用,(计算机科学讲义,第974卷(1983),施普林格:施普林格柏林)·Zbl 0504.42001
[9] Higham,N.J.,Matlab测试矩阵工具箱(3.0版),(第276号数值分析报告(1995年9月),曼彻斯特大学数学系)
[10] Lanczos,C.,《通过最小化迭代求解线性方程组》,J.Res.Natl。伯尔。支架。,49, 33-53 (1952)
[11] O'Leary,D.P.,块共轭梯度算法及相关方法,线性代数应用。,13, 567-593 (1992)
[12] Saad,Y.,《稀疏线性系统的迭代方法》(1995),PWS:PWS Boston·Zbl 1002.65042号
[13] H.Sadok,K.Jbilou,全球Lanczos型方法及其应用,Appl。线性代数,提交。;H.Sadok,K.Jbilou,全球Lanczos型方法及其应用,Appl。线性代数,提交·Zbl 1308.65007号
[14] Simoncini,V.,块BiCG的稳定QMR版本,SIAM J.矩阵分析。申请。,18419-434(1997年)·Zbl 0872.65024号
[15] 西蒙西尼,V。;Gallopoulos,E.,《具有多个右手边的非对称系统的迭代方法》,SIAM J.Sci。计算。,16, 917-933 (1995) ·Zbl 0831.65037号
[16] 西蒙西尼,V。;Gallopoulos,E.,《具有多个右手边的非对称系统的混合块GMRES方法》,J.Compute。申请。数学。,66, 1-2, 457-469 (1996) ·Zbl 0859.65021号
[17] 西蒙西尼,V。;Gallopoulos,E.,求解具有多个右手边的系统的块GMRES的收敛性,(技术代表1316(1993年10月),超级计算研究与发展中心,伊利诺伊大学厄巴纳-香槟分校)·Zbl 0859.65021号
[18] Sonneveld,P.,CGS,非对称线性系统的快速Lanczos型解算器,SIAM J.Sci。统计计算。,10, 35-52 (1989) ·Zbl 0666.65029号
[19] van der Vorst,H.A.,Bi-CGSTAB:非对称线性系统解的Bi-CG的一个快速且平滑收敛的变体,SIAM J.Sci。统计计算。,13, 631-644 (1992) ·Zbl 0761.65023号
[20] van der Vorst,H.A.,利用对称正定矩阵的Krylov子空间信息求解(A)x=b的迭代方法,J.Compute。申请。数学。,18, 249 (1987) ·Zbl 0621.65022号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。