×

矩阵函数的块Krylov子空间方法。二: 修改块FOM。 (英语) Zbl 1441.65051号

摘要:我们分析了我们工作第一部分的广义块Krylov子空间框架的扩展[ETNA,Electron.Trans.Numer.Anal.47100-126(2017;Zbl 1372.65138号)]. 这种扩展允许使用投影到块Krylov子空间的矩阵的低秩修改,并且作为特殊情况,包含块GMRES方法和新的块Radau-Arnoldi方法。在这种一般设置下,我们给出了将插值性质从非块情形扩展到块情形的矩阵多项式插值性质的结果,并且我们将投影矩阵的特征值与这些矩阵多项式的潜在根联系起来。给出了求解线性系统的改进块FOM方法的误差界。然后我们展示如何化妆品在移位线性块系统族的情况下,可以保留残差。该结果用于推导块Krylov近似的计算上实用的重新启动算法,该算法同时计算矩阵函数对几个向量集的作用。我们证明了一些误差界,并给出了数值结果,表明对FOM的两种修改,即矩阵函数的块调和方法和块Radau-Arnoldi方法,可以显著改善收敛性能。

MSC公司:

65层60 矩阵指数和相似矩阵函数的数值计算
65层50 稀疏矩阵的计算方法
65层55 低阶矩阵逼近的数值方法;矩阵压缩
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] O.Abidi、M.Heyouni和K.Jbilou,关于扩展块和全局Arnoldi方法的一些性质及其在模型约简中的应用,Numer。算法,75(2017),第285-304页,https://doi.org/10.1007/s11075-016-0207-7。 ·Zbl 1364.65077号
[2] A.H.Al-Mohy和N.J.Higham,利用指数积分器计算矩阵指数的作用,SIAM J.Sci。计算。,33(2011),第488-511页,https://doi.org/10.1137/100788860。 ·Zbl 1234.65028号
[3] T.Bakhos、P.K.Kitanidis、S.Ladenheim、A.K.Saibaba和D.B.Szyld,移位系统的多修复GMRES,SIAM J.Sci。计算。,19122(2017),第S222-S247页,https://doi.org/10.1137/090750688。 ·Zbl 1392.65043号
[4] B.Beckermann、D.Kressner和M.Schweitzer,矩阵函数的低库更新,SIAM J.矩阵分析。申请。,39(2018),第539-565页,https://epubs.siam.org/doi/10.1137/17M1140108。 ·Zbl 1390.15024号
[5] B.Beckermann和L.Reichel,通过Faber变换对矩阵函数进行误差估计和评估,SIAM J.Numer。分析。,47(2009),第3849-3883页,https://doi.org/10.1137/080741744。 ·Zbl 1204.65041号
[6] F.P.A.Beik和D.K.Salkuyeh,关于求解一般耦合矩阵方程的全局Krylov子空间方法,计算。数学。申请。,62(2011),第4605-4613页,https://doi.org/10.1016/j.camwa.2011.10.043。 ·Zbl 1236.15031号
[7] S.Birk,Hermitian正定矩阵的收缩移位块Krylov子空间方法,博士论文,Fakulta¨t fu¨r Mathematik und Naturwissenschaften,Bergische Universita¨t Wupertal,2015,https://d-nb.info/1073127559/34。
[8] J.C.Bloch、T.Breu、A.Frommer、S.Heybrock、K.Scha-fer和T.Wettig,Krylov子空间方法和符号函数:非Hermitian情况下的多重移位和通缩,Proc。科学。,91 (2009), https://pos.sissa.it/091/043。 ·Zbl 1219.81096号
[9] R.Bouyouli,K.Jbilou,R.Sadaka,和H.Sadok,多线性系统某些块Krylov子空间方法的收敛性,J.Compute。申请。数学。,196(2006),第498-511页,https://doi.org/10.1016/j.cam.2005.09.017。 ·Zbl 1100.65024号
[10] T.A.Davis和Y.Hu,佛罗里达大学稀疏矩阵收集,ACM Trans。数学。软件,38(2011),第1-25页,https://doi.org/10.1145/2049662.2049663。 ·Zbl 1365.65123号
[11] J.E.Dennis、J.F.Traub和R.P.Weber,《关于矩阵多项式、Lambda-Matrix和块特征值问题》,技术报告CMU-CS-71-110,卡内基梅隆大学计算机科学系,1971年,https://ecommons.cornell.edu/handle/1813/5954。
[12] J.E.Dennis、J.F.Traub和R.P.Weber,矩阵多项式的代数理论,SIAM J.Numer。分析。,13(1976年),第831-845页,https://epubs.siam.org/doi/10.1137/0713065。 ·兹伯利0361.15013
[13] J.E.Dennis、J.F.Traub和R.P.Weber,矩阵多项式溶剂的算法,SIAM J.Numer。分析。,15(1978年),第523-533页,https://epubs.siam.org/doi/101137/0715034。 ·Zbl 0386.65012号
[14] V.L.Druskin和L.A.Knizhnerman,计算对称矩阵函数的两种多项式方法,U.S.S.R.计算。数学。数学。物理。,29(1989),第112-121页,https://doi.org/10.1016/S0041-5553(89)80020-5. ·Zbl 0719.65035号
[15] S.C.Eisenstat、H.C.Elman和M.H.Schultz,非对称线性方程组的变分迭代方法,SIAM J.Numer。分析。,20(1983年),第345-357页,https://epubs.siam.org/doi/10.1137/0720023。 ·Zbl 0524.65019号
[16] L.Elbouyahyaoui、A.Messaoudi和H.Sadok,块GMRES和块Arnoldi方法的代数性质,电子。事务处理。数字。分析。,33(2008),第207-220页,https://etna.ricam.eoaw.ac.at/vol.33.2008-2009/pp207-220.dir/pp207-220.pdf。 ·Zbl 1188.65031号
[17] A.Frommer和U.Gla¨ssner,移位线性系统的重新启动GMRES,SIAM J.Sci。计算。,19(1998),第15-26页,https://doi.org/10.1137/S1064827596304563。 ·Zbl 0912.65023号
[18] A.Frommer、S.Gu¨sken、T.Lippert、B.No-ckel和K.Schilling,《一击多质量:夸克传播者的经济计算》,国际。现代杂志。物理。C、 6(1995),第627-638页,https://doi.org/10.1142/S0129183195000538。
[19] A.Frommer、S.Guöttel和M.Schweitzer,矩阵Stieltjes函数重启Krylov子空间方法的收敛性,SIAM J.矩阵分析。申请。,35(2014),第1602-1624页,https://doi.org/10.1137/10973463。 ·Zbl 1316.65040号
[20] A.Frommer、S.Guöttel和M.Schweitzer,基于求积的矩阵函数的高效稳定Arnoldi重启,SIAM J.matrix Anal。申请。,35(2014),第661-683页,https://epubs.siam.org/doi/abs/10.1137/13093491X。 ·Zbl 1309.6500号
[21] A.Frommer、K.Lund、M.Schweitzer和D.B.Szyld,矩阵函数的Radau–Lanczos方法,SIAM J.matrix Ana。申请。,38(2017),第710-732页,http://epubs.siam.org/doi/10.1137/16M1072565。 ·Zbl 1371.65042号
[22] A.Frommer、K.Lund和D.B.Szyld,矩阵函数的块Krylov子空间方法,电子。事务处理。数字。分析。,47(2017),第100-126页,http://etna.mcs.kent.edu/vol.47.2017/pp100-126.dir/pp100-126.5df。 ·Zbl 1372.65138号
[23] C.Gattringer和C.B.Lang,《晶格上的量子色动力学:介绍性演讲》,Springer,柏林,2010年,https://doi.org/10.1007/978-3642-01850-3。
[24] I.Gohberg、P.Lancaster和L.Rodman,矩阵多项式,应用中的经典。数学。58,SIAM,费城,2009年·Zbl 1170.15300号
[25] G.H.Golub和G.Meurant,《矩阵、矩和求积及其应用》,普林斯顿大学出版社,新泽西州普林斯顿,2010年·兹比尔1217.65056
[26] M.H.Gutknecht,《具有多个右手边的线性系统的Block Krylov空间方法:现实世界系统的现代数学模型、方法和算法简介》,A.H.Siddiqi、I.S.Duff和O.Christensen编辑,阿纳马亚,新德里,2007年,第420-447页。
[27] N.Hale、N.J.Higham和L.N.Trefethen,《计算》(A^{alpha},log(A)),以及通过轮廓积分实现的相关矩阵函数,SIAM J.Numer。分析。,46(2008),第2505-2523页,https://doi.org/10.1137/070700607。 ·Zbl 1176.65053号
[28] P.Henrici,《应用和计算复杂性分析》,第2卷,John Wiley&Sons,纽约,1977年·兹比尔0363.30001
[29] M.Heyouni和A.Essai,具有多个右手边的线性系统的矩阵Krylov子空间方法,Numer。算法,40(2005),第137-156页,https://doi.org/10.1007/s11075-005-1526-2。 ·Zbl 1087.65028号
[30] N.J.Higham,《矩阵的函数:理论与计算》,SIAM,费城,2008年·Zbl 1167.15001号
[31] M.Hochbruck和A.Ostermann,指数积分器,Acta Numer。,19(2010),第209-286页,https://doi.org/10.1017/S0962492910000048。 ·兹比尔1242.65109
[32] K.Jbilou、A.Messaoudi和H.Sadok,矩阵方程的全局FOM和GMRES算法,应用。数字。数学。,31(1999),第49-63页,https://doi.org/10.1016/S0168-9274(98)00094-4. ·Zbl 0935.65024号
[33] M.D.Kent、Chebyshev、Krylov、Lanczos:矩阵关系和计算,加州斯坦福大学计算机科学系博士论文,1989年,https://dl.acm.org/doi/book/10.5555/76475。
[34] P.Lancaster,《Lambda-计量与振动系统》,英国牛津佩加蒙出版社,1966年·Zbl 0146.32003号
[35] J.Liesen和Z.Strakoš,《Krylov子空间方法:原理与分析》,牛津大学出版社,英国牛津,2013年·Zbl 1263.65034号
[36] H.Liu和B.Zhong,具有多个右手边的非对称系统的更简单块GMRES,Electron。事务处理。数字。分析。,30(2008),第1-9页,http://etna.mcs.kent.edu/vol.30.2008/pp1-9.dir/pp1-9.pdf。 ·Zbl 1188.65034号
[37] K.Lund,《一个新的块Krylov子空间框架及其对作用于多个向量的矩阵函数的应用》,Temple大学数学系博士论文,以及Fakulta­t Mathematik und Naturwissenschaften der Bergischen Universita­t Wuppertal,2018年,https://digital.library.temple.edu/digital/collection/p245801coll10/id/493125/rec/1, http://elpub.bib.uni-woppertal.de/edocs/dokumente/fbc/mathematik/dis2018/lund?lang=en。
[38] D.P.O'Leary,块共轭梯度算法及相关方法,线性代数应用。,29(1980),第293-322页,https://doi.org/10.1016/0024-3795(80)90247-5. ·Zbl 0426.65011号
[39] C.C.Paige、B.N.Parlett和H.A.van der Vorst,Krylov子空间的近似解和特征值界,Numer。线性代数。申请。,2(1995),第115-133页,https://doi.org/10.1002/nla.1680020205。 ·Zbl 0831.65036号
[40] S.Rashedi、G.Ebadi、S.Birk和A.Frommer,关于具有多个右手边的线性系统的短递归Krylov型方法,J.Compute。申请。数学。,300(2016),第18-29页,https://doi.org/10.1016/j.cam.2015.11.040。 ·Zbl 1382.65092号
[41] Y.Saad,矩阵指数算子的Krylov子空间近似分析,SIAM J.Numer。分析。,29(1992),第209-228页,https://epubs.siam.org/doi/abs/10.1137/0729014。 ·Zbl 0749.65030号
[42] Y.Saad,稀疏线性系统的迭代方法,第2版,SIAM,费城,2003年·Zbl 1031.65046号
[43] Y.Saad和M.H.Schultz,GMRES:求解非对称线性系统的广义最小残差算法,SIAM J.Sci。统计计算。,7(1986),第856-869页,http://epubs.siam.org/doi/10.1137/0907058。 ·Zbl 0599.65018号
[44] A.K.Saibaba、T.Bakhos和P.K.Kitanidis,用于振荡液压断层扫描的移动系统的灵活Krylov解算器,SIAM J.Sci。计算。,35(2013),第A3001-A3023页,https://doi.org/10.1137/120902690。 ·Zbl 1288.65041号
[45] V.Simoncini,Ritz和Pseudo-Ritz值使用矩阵多项式,线性代数应用。,241-243(1996),第787-801页,https://doi.org/10.1016/0024-3795(95)00682-6. ·Zbl 0856.65036号
[46] V.Simoncini和E.Gallopoulos,块GMRES和矩阵多项式的收敛性,线性代数应用。,247(1996),第97-119页,https://doi.org/10.1016/0024-3795(95)00093-3. ·Zbl 0861.65023号
[47] K.M.Soodhalter,不相关右手边移位系统的块Krylov子空间循环,SIAM J.Sci。计算。,38(2016),第A302-A324页,http://epubs.siam.org/doi/pdf/10.1137/10998214。
[48] G.Starke,非对称椭圆问题预处理迭代方法的场值分析,Numer。数学。,78(1997),第103-117页,https://doi.org/10.1007/s002110050306。 ·Zbl 0888.65037号
[49] D.-l.Sun,T.-Z.Huang,Y.-F.Jing,B.Carpentieri,求解具有多个移位和多个右手边的线性系统的带放气重启的块GMRES方法,Numer。线性代数。申请。,25(2018),e2148,https://doi.org/10.1002/nla.2148。 ·Zbl 1513.65085号
[50] J.van den Eshof、A.Frommer、T.Lippert、K.Schilling和H.A.van der Vorst,QCDd重叠算子的数值方法:I.符号函数和误差界,计算。物理。社区。,146(2002),第203-224页,https://doi.org/10.1016/S0010-4655(02)00455-1. ·Zbl 1008.81508号
[51] B.Vital,E⁄tude quelques me⁄thodes de re⁄solution de problèmes line⁄aires de grand taille sur multiprocesseur,博士论文,雷恩大学,1990年,http://www.sudoc.fr/044024320。
[52] G.Wu,Y.-c.Wang,X.-Q.Jin,多阻尼因子PageRank问题的预处理和移位GMRES算法,SIAM J.Sci。计算。,34(2012),第A2558-A2575页,http://epubs.siam.org/doi/abs/10.1137/10834585。 ·Zbl 1263.65037号
[53] J.Zhang、H.Dai和J.Zhao,《具有多个右手边的线性系统的新全局方法家族》,J.Compute。申请。数学。,236(2011),第1562-1575页,https://doi.org/10.1016/j.cam.2011.09.20。 ·Zbl 1252.65072号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。