×

不定矩阵ILU预条件的实验研究。 (英语) Zbl 0891.65028号

不完全LU分解预条件已成功地应用于许多一般非对称和不定矩阵的情况。然而,为了使它们对一般矩阵有用,仍然存在一些数值困难,包括非常小的枢轴、不稳定的三角形解算器、下降引起的不准确性以及零枢轴。此外,这四个问题通常同时出现。因此,有必要更好地了解LU预条件并帮助提高其可靠性。
本文展示了这些问题是如何表现出来的,如何检测这些问题,以及有时如何通过旋转、缩放、重新排序、扰动对角元素和保持对称结构来规避这些问题。数值实验表明,所提出的技术是有效的。

MSC公司:

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

参考文献:

[1] Axelsson,O.,《迭代求解方法》(1994),剑桥大学出版社:剑桥大学出版社·Zbl 0795.65014号
[2] O.阿克塞尔森。;Kolotilina,L.Yu。,对角补偿还原和相关预处理方法,编号。线性代数应用。,1, 155-177 (1994) ·Zbl 0837.65023号
[3] O.阿克塞尔森。;Lindskog,G.,关于一类预处理方法的特征值分布,Numer。数学。,48679-498(1986年)·Zbl 0564.65016号
[4] O.阿克塞尔森。;Lindskog,G.,关于预处理共轭梯度法的收敛速度,Numer。数学。,48, 499-523 (1986) ·Zbl 0564.65017号
[5] 博塔,E.F.F。;van der Ploeg,A。;Wubs,F.W.,嵌套网格ILU分解(NGILU),J.Comp。申请。数学。,66, 515-526 (1996) ·Zbl 0858.65028号
[6] 布鲁塞特,A.M。;特维托,A。;Winther,R.,关于松弛不完全LU分解的稳定性,数学。公司。,54, 701-719 (1990) ·Zbl 0701.65021号
[7] Buleev,N.I.,(代表BNL-TR-551(1973),布鲁克海文国家实验室:布鲁克海温国家实验室,纽约厄普顿),英语翻译·Zbl 0105.10401号
[8] 邦奇,J.R。;Kaufman,L.,计算惯性和求解对称线性系统的一些稳定方法,数学。公司。,31, 162-179 (1977) ·Zbl 0355.65023号
[9] Chan,T.F。;Van der Vorst,H.A.,近似和不完全因子分解,(技术报告871(1994),数学系:乌得勒支大学数学系)·Zbl 0865.65015号
[10] 查普曼,A。;萨阿德,Y。;Wigton,L.,CFD问题的高阶ILU预处理程序,(技术报告UMSI 96/14(1996),明尼苏达大学明尼阿波利斯分校明尼苏达到超级计算机研究所:明尼苏那大学明尼苏打超级计算机研究院)·Zbl 0959.76077号
[11] 周,E。;Saad,Y.,ILUS:稀疏天际线格式的不完全LU预处理程序,国际。J.数字。液体方法,25739-748(1997)·Zbl 0896.76037号
[12] D'Azevdo,E.F。;福赛斯,P.A。;Tang,W.-P.,《朝着一种具有成本效益的高填方ILU预处理剂迈进》,BIT,32,442-463(1992)·Zbl 0761.65017号
[13] 达夫,I.S。;新墨西哥州古尔德。;Reid,J.K。;Scott,J.A.,稀疏不定矩阵的因式分解,IMA J.Numer。分析。,11, 181-204 (1991) ·Zbl 0739.65018号
[14] 达夫,I.S。;Meurant,G.A.,排序对预处理共轭梯度的影响,BIT,29635-657(1989)·兹伯利0687.65037
[15] Dutto,L.C.,用于求解可压缩Navier-Stokes方程的预处理GMRES算法的排序影响,国际。J.数字。方法工程,36,457-497(1993)·Zbl 0767.76026号
[16] Elman,H.C.,不完全LU因子分解的稳定性分析,数学。公司。,47, 191-217 (1986) ·Zbl 0621.65024号
[17] Elman,H.C.,非自伴线性系统的松弛和稳定不完全因子分解,BIT,29890-915(1989)·Zbl 0694.65011号
[18] Engelman,M.,FIDAP:示例手册,6.0版(1991),流体动力学国际:流体动力学国际,伊利诺伊州埃文斯顿
[19] Engelman,医学硕士。;Hasbani,I.,有限元环境下无矩阵解算法(技术报告88-1(1988),流体动力学国际:伊利诺伊州埃文斯顿流体动力学国际)
[20] 乔治,A。;Liu,J.W.,大型稀疏正定系统的计算机解(1981),Prentice-Hall:Prentice-Hall Englewood Cliffs,NJ·Zbl 0516.65010号
[21] Gustafsson,I.,一类一阶因式分解方法,BIT,18,142-156(1978)·Zbl 0386.65006号
[22] A.詹宁斯。;Malik,G.M.,部分消除,J.Inst.数学应用。,20, 307-316 (1977) ·Zbl 0367.65019号
[23] Kershaw,D.S.,线性方程组迭代解的不完全Cholesky共轭梯度法,J.Compute。物理。,26, 43-65 (1978) ·Zbl 0367.65018号
[24] Manteuffel,T.A.,正定线性系统的不完全因式分解技术,数学。公司。,34, 473-497 (1980) ·Zbl 0422.65018号
[25] Meijerink,J.A。;Van der Vorst,H.A.,系数矩阵为对称M矩阵的线性系统的迭代求解方法,数学。公司。,31, 148-162 (1977) ·Zbl 0349.65020号
[26] Meijerink,J.A。;Van der Vorst,H.A.,《在解决实际问题中出现的线性方程组时使用不完全分解的指南》,J.Compute。物理。,44, 134-155 (1981) ·Zbl 0472.65028号
[27] Munksgaard,N.,通过预处理共轭梯度求解线性方程的稀疏对称集,ACM Trans。数学。软件,6206-219(1980)·Zbl 0438.65035号
[28] Oliphant,T.A.,求解二维含时扩散问题的隐式数值方法,Quart。申请。数学。,19, 221-229 (1961) ·Zbl 0105.10702号
[29] Oliphant,T.A.,解线性系统的外推过程,夸特。申请。数学。,20, 257-267 (1962) ·Zbl 0109.34601号
[30] Ortega,J.M.,《线性系统的并行和矢量解导论》(1988),Plenum出版社:Plenum Press New York·兹伯利0669.65017
[31] Ø斯特比,O。;Zlatev,Z.,(稀疏矩阵的直接方法,计算机科学课堂讲稿,第157卷(1983年),施普林格:施普林格-柏林)·兹伯利0516.65011
[32] Parter,S.V.,《线性图在高斯消去中的应用》,SIAM Rev.,3119-130(1961)·Zbl 0102.11302号
[33] Robert,Y.,实正定矩阵的正则不完全因式分解,线性代数应用。,48105-117(1982年)·Zbl 0502.65018号
[34] Saad,Y.,《不定和非对称线性系统的预处理技术》,J.Compute。申请。数学。,24, 89-105 (1988) ·Zbl 0662.65028号
[35] Saad,Y.,ILUT:双阈值不完全ILU因式分解,Numer。线性代数应用。,1, 387-402 (1994) ·Zbl 0838.65026号
[36] Saad,Y.,CFD应用的预处理Krylov子空间方法,(Habashi,W.G.,Proc.Internat.大型CFD问题解决技术研讨会,Proc.Internat.大型CFD解决技术研讨会
[37] Van der Vorst,H.A.,PDE问题产生的非对称矩阵的某些稀疏线性系统的迭代求解方法,J.Compute。物理。,44, 1-19 (1981) ·Zbl 0472.65027号
[38] Van der Vorst,H.A.,稳定的不完全LU-分解作为切比雪夫迭代的预处理,(Evans,D.J.,预处理方法:分析与应用(1983),Gordon and Breach:Gordon和Breach New York)·Zbl 0767.65024号
[39] Varga,R.S.,因式分解和归一化迭代方法,(Langer,R.E.,微分方程中的边界问题(1960),威斯康星大学出版社:威斯康辛大学出版社,加州圣巴巴拉),121-142·Zbl 0100.12501号
[40] 瓦尔加,R.S。;Saff,E.B。;Mehrman,V.,矩阵的不完全分解和与H-矩阵的连接,SIAM J.Numer。分析。,17, 787-793 (1980) ·Zbl 0477.65020号
[41] Watts-,J.W.,油藏模拟压力方程迭代解的共轭梯度截断直接法,Soc.Pet。《工程师杂志》,21,345-353(1981)
[42] Wittum,G.,《关于ILU平滑的稳健性》,SIAM J.Sci。统计计算。,10, 699-717 (1989) ·Zbl 0677.65096号
[43] A.余。Yeremin,私人通信,1995年。;A.余。Yeremin,私人通信,1995年。
[44] Zlatev,Z。;巴克,V.A。;Thomsen,P.G.,SLEST:用于求解稀疏线性方程组的Fortran IV子例程,用户指南,(NI-78-01(1978)技术报告,Numerisk Institute:丹麦林比Numerisk-Institute Lyngby)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。