×

使用部分矩阵估计的无矩阵预处理。 (英语) Zbl 1105.65046号

摘要:我们考虑无矩阵求解器环境,其中有关底层矩阵的信息只能通过矩阵向量计算获得,而这些计算无法访问完全组装的矩阵。我们引入了部分矩阵估计的概念,用于构造在这种无矩阵环境下Krylov迭代方法中使用的好的代数预条件,并提出了三个新的部分矩阵估计图着色问题。使用这些公式之一的数值实验证明了该方法的可行性。

MSC公司:

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

参考文献:

[1] M.A.Ajiz和A.Jennings,《稳健的不完全Choleski共轭梯度算法》,国际期刊Numer。《方法工程》,20(1984),第949–966页·Zbl 0541.65019号 ·doi:10.1002/nme.1620200511
[2] O.Axelsson,《迭代求解方法》,剑桥大学出版社,剑桥,1994年·Zbl 0795.65014号
[3] O.Axelsson和L.Kolotilina,对角补偿还原和相关预处理方法,数值。林氏代数应用。,1(1994),第155–177页·Zbl 0837.65023号 ·doi:10.1002/nla.1680010207
[4] M.W.Benson,大尺度线性系统的迭代解,硕士论文,莱克黑德大学,雷霆湾,1973年。
[5] M.Benzi和M.Tuma,非对称线性系统的稀疏近似逆预条件,SIAM J.Sci。计算。,19(1998),第968–994页·Zbl 0930.65027号 ·doi:10.1137/S1064827595294691
[6] E.Chow,并行稀疏近似逆预条件的先验稀疏模式,SIAM J.Sci。计算。,21(2000),第1804-1822页·兹比尔0957.65023 ·doi:10.1137/S106482759833913X
[7] E.Chow,ParaSails用户指南,技术代表UCRL-MA-137863,劳伦斯·利弗莫尔国家实验室,加利福尼亚州利弗莫尔市,2000年。
[8] E.Chow,具有先验稀疏模式的稀疏近似逆的并行实现和实际使用,国际高性能计算杂志。申请。,15(2001),第56-74页·doi:10.1177/10943420101500106
[9] T.F.Coleman和J.Cai,稀疏Hessian矩阵的循环着色问题和估计,SIAM J.离散算法。方法,7(1986),第221-235页·Zbl 0613.65066号 ·doi:10.1137/0607026
[10] T.F.Coleman、B.S.Garbow和J.J.Moré,稀疏雅可比矩阵估计软件,ACM Trans。数学。软质。,10(1984年),第329-345页·Zbl 0548.65006号 ·数字对象标识代码:10.1145/1271.1610
[11] T.F.Coleman和J.J.Moré,稀疏雅可比矩阵的估计和图着色问题,SIAM J.Numer。分析。,20(1983年),第187-209页·Zbl 0527.65033号 ·doi:10.1137/0720013
[12] T.F.Coleman和J.J.Moré,稀疏Hessian矩阵的估计和图着色问题,数学。程序。,28(1984年),第243-270页·Zbl 0572.65029号 ·doi:10.1007/BF02612334
[13] T.F.Coleman和A.Verma,《结构和有效雅可比计算》,载于《计算微分:技术、应用和工具》,Berz等人,第149-159页,SIAM,费城,1996年·Zbl 0869.65012号
[14] T.F.Coleman和A.Verma,使用自动微分有效计算稀疏雅可比矩阵,SIAM J.Sci。计算。,19(1998),第1210–1233页·Zbl 0913.65041号 ·doi:10.1137/S1064827595295349
[15] A.R.Curtis、M.J.D.Powell和J.K.Reid,关于稀疏雅可比矩阵的估计,数学研究所。申请。,13(1974年),第117-119页·Zbl 0273.65036号
[16] T.Davis,稀疏矩阵收集,《北美文摘》,42(1994),http://www.cise.ufl.edu/avis/sparese/。
[17] T.F.Dupont,R.P.Kendall和H.H.Rachford,解自伴椭圆差分方程的近似因子分解过程,SIAM J.Numer。分析。,5(1968年),第559-573页·Zbl 0174.47603号 ·doi:10.1137/0705045
[18] A.H.Gebremedhin、F.Manne和A.Pothen,你的Jacobian是什么颜色的?计算导数的图形着色,SIAM Rev.,47(2005),第629–705页·Zbl 1076.05034号 ·doi:10.1137/S0036144504444711
[19] K.Georg,无矩阵数值延拓与分岔,数值。功能。分析。优化,22(2001),第303–320页·Zbl 0998.65055号 ·doi:10.1081/NFA-100105106
[20] A.Griewank,《不累加雅可比数的牛顿步长直接计算》,载于《大尺度数值优化》,T.F.Coleman和Y.Li,eds.,第115–137页,SIAM,费城,1990年·Zbl 0726.90080号
[21] A.Griewank,梯度复杂性的一些边界,Jacobians和Hessians,《数值优化中的复杂性》,P.M.Pardalos主编,世界科学出版社,新泽西州River Edge,1993年·Zbl 0968.90524号
[22] I.Gustafsson,一类一阶因式分解方法,BIT,18(1978),第142-156页·Zbl 0386.65006号 ·doi:10.1007/BF01931691
[23] S.Hossain,《关于稀疏雅可比矩阵和牛顿步长的计算》,技术代表第146号,卑尔根大学信息学系,1998年。
[24] S.Hossain和T.Steihaug,图着色和分段列稀疏雅可比矩阵的估计,技术代表第72号,卑尔根大学信息学系,1997年。
[25] S.Hossain和T.Steihaug,雅可比矩阵计算中的稀疏性问题,载于《2002年符号和代数计算国际研讨会论文集》,T.Mora主编,第123–130页,ACM出版社,纽约,2002年·Zbl 1072.68671号
[26] S.Hossain和T.Steihaug,《减少计算稀疏雅可比矩阵的AD通过次数》,载于《自动微分:从模拟到优化》,G.Corliss、C.Faure、a.Griewank、L.HascoéT和U.Naumann编辑,《计算机与信息科学》,纽约斯普林格出版社,2002年。
[27] T.Huckle,矩阵逆和预处理的近似稀疏模式,载于1997年第15届IMACS科学计算、建模和应用数学世界大会论文集,第二卷,第569-574页,IMACS,新不伦瑞克,1997年。
[28] D.Hysom和A.Pothen,不完全因子预处理的可扩展并行算法,SIAM J.Sci。计算。,22(2001),第2194-2215页·Zbl 0986.65048号 ·doi:10.1137/S1064827500376193
[29] D.Hysom和A.Pothen,ILU(k)预条件的高效并行计算,《超级计算学报》,1999年,第1-19页,IEEE出版社,新泽西州皮斯卡塔韦,1999年。
[30] G.Karypis和V.Kumar,基于并行阈值的ILU因子分解,技术代表TR-96-061,明尼苏达大学,1996年。
[31] C.T.Kelley,用牛顿法求解非线性方程,SIAM,费城,2003年·Zbl 1031.65069号
[32] D.E.Keyes,偏微分方程的Terascale隐式方法,Contemp。数学。,306(2001),第29-84页·Zbl 1018.65121号
[33] D.A.Knoll和D.E.Keyes,《无雅可比牛顿-Krylov方法:方法和应用综述》,J.Comput。物理。,193(2004),第357-397页·Zbl 1036.65045号 ·doi:10.1016/j.jcp.2003.08.010
[34] L.Y.Kolotilina、A.A.Nikishin和A.Yu。Yeremin,分解稀疏近似逆预处理。IV: 提高效率的简单方法,Numer。线性代数应用。,6(1999年),第515-531页·Zbl 0983.65062号 ·doi:10.1002/(SICI)1099-1506(199910/11)6:7<515::AID-NLA176>3.0.CO;2-0
[35] G.Meurant,《大型线性系统的计算机解决方案,数学及其应用研究》,第28期,爱思唯尔,阿姆斯特丹-洛桑-纽约-牛津-香农-新加坡-东京,1999年·兹比尔0954.65025
[36] J.L.Morales和J.Nocedal,有限内存准Newton更新的自动预处理,SIAM J.Optim。,10(2000),第1079–1096页·Zbl 1020.65019号 ·doi:10.1137/S1052623497327854
[37] J.L.Morales和J.Nocedal,算法809:PREQN:Fortran 77共轭梯度法预处理,ACM Trans。数学。软质。,27(2001),第83-91页·Zbl 1070.65538号 ·数字对象标识代码:10.1145/382043.382343
[38] P.E.Plassman,多处理器上的稀疏雅可比估计和因式分解,载于《大尺度优化》,T.F.Coleman和Y.Li编辑,第152-179页,SIAM,费城,1990年·Zbl 0726.90078号
[39] Y.Saad,一种灵活的内外预处理GMRES算法,SIAM J.Sci。统计计算。,14(1993),第461-469页·Zbl 0780.65022号 ·doi:10.1137/0914028
[40] Y.Saad,ILUT:双阈值不完全LU因式分解,Numer。线性代数应用。,1(1994年),第387-402页·Zbl 0838.65026号 ·doi:10.1002/nla.1680010405
[41] Y.Saad,《稀疏线性系统的迭代方法》,第2版,SIAM,费城,2003年·Zbl 1031.65046号
[42] P.L.Toint,关于线性方程的稀疏对称矩阵更新,数学。计算。,31(1977年),第954-961页·Zbl 0379.65034号 ·doi:10.1090/S0025-5718-1977-0455338-4
[43] D.A.Korzekwa和A.David,《Truchas固体力学能力的现状》,技术代表LA-UR-04-4092,洛斯阿拉莫斯国家实验室,2004年。
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。