×

移位和反转Krylov子空间的精确重启计算非对称矩阵的矩阵指数作用。 (英语。俄文原件) 兹比尔1469.65087

计算。数学。数学。物理学。 第561684-698号(2021); Zh的翻译。维奇尔。Mat.Mat.Fiz公司。61,第5期,706-722(2021)。
摘要:提出了一种精确的剩余时间(AccuRT)重启动方法,用于利用移位和反转(SAI)Krylov子空间方法计算非对称矩阵的矩阵指数作用。所提出的重启方法是最近提出的RT(剩余时间)重启的扩展,旨在避免常规RT重启中可能出现的精度损失。SAI Krylov方法的一个昂贵部分是用移位矩阵求解线性系统。由于AccuRT算法调整移位值,我们讨论了如何通过移位矩阵的单个LU分解(或预条件设置)来实现所建议的重启。数值实验表明,该方法提高了精度和效率。

MSC公司:

65层60 矩阵指数和相似矩阵函数的数值计算
15甲16 矩阵的指数函数和相似函数
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] Hochbruck,M。;Ostermann,A.,指数积分器,《数值学报》。,19, 209-286 (2010) ·Zbl 1242.65109号 ·doi:10.1017/S0962492910000048
[2] de la Cruz Cabrera,俄亥俄州。;马塔尔,M。;Reichel,L.,通过矩阵指数分析有向网络,J.Compute。申请。数学。,355, 182-192 (2019) ·Zbl 1417.90041号 ·doi:10.1016/j.cam.2019.01.015
[3] P.Kürschner,“大型系统在有限时间间隔内的平衡截断模型降阶”,高级计算。数学。44 (2018). doi:10.1007/s10444-018-9608-6·Zbl 1453.65093号
[4] Frommer,A。;Simoncini,V.,《模型降阶:理论、研究方向和应用》(2008),柏林:施普林格出版社,柏林·Zbl 1203.65079号
[5] Bergamaschi,L。;Vianello,M.,大型稀疏对称矩阵指数算子的高效计算,Numer。线性代数应用。,7, 27-45 (2000) ·Zbl 0983.65058号 ·doi:10.1002/(SICI)1099-1506(200001/02)7:1<27::AID-NLA185>3.0.CO;2-4
[6] Al-Mohy,A.H。;Higham,N.J.,计算矩阵指数的作用,以及指数积分器的应用,SIAM J.Sci。计算。,33, 488-511 (2011) ·Zbl 1234.65028号 ·doi:10.1137/100788860
[7] 莫雷特,I。;Novati,P.,矩阵指数的RD有理逼近,BIT,44595-615(2004)·兹比尔1075.65062 ·doi:10.1023/B:BITN.0000046805.27551.3b
[8] van den Eshof,J。;Hochbruck,M.,将Lanczos近似预处理为矩阵指数,SIAM J.Sci。计算。,27, 1438-1457 (2006) ·Zbl 1105.65051号 ·数字对象标识代码:10.1137/040605461
[9] 德鲁斯金,V。;Knizhnerman,L。;Zaslavsky,M.,使用带优化移位的有理Krylov子空间解决大规模进化问题,SIAM J.Sci。计算。,31, 3760-3780 (2009) ·Zbl 1204.65042号 ·数字对象标识代码:10.1137/080742403
[10] S.Güttel,博士论文(弗赖堡理工大学Bergakademie,2010年)。www.guettel.com。
[11] S.Güttel,“矩阵函数的有理Krylov近似:数值方法和最优极点选择”,GAMM Mitteilungen 36(1),8-31(2013)。www.guettel.com·Zbl 1292.65043号
[12] Celledoni,E。;Moret,I.,ODE系统的Krylov投影方法,应用。数字。数学。,24, 365-378 (1997) ·Zbl 0885.65073号 ·doi:10.1016/S0168-9274(97)00033-0
[13] Tal-Ezer,H.,关于w=f(A)\({v}\)的Krylov近似的重启和误差估计,SIAM J.Sci。计算。,29, 2426-2441 (2007) ·Zbl 1154.65320号 ·doi:10.1137/040617868
[14] Afanasjew,M。;艾尔曼,M。;Ernst,O.G。;Güttel,S.,矩阵函数计算的重新启动Krylov子空间方法的实现,线性代数应用。,429, 2293-2314 (2008) ·Zbl 1153.65042号 ·doi:10.1016/j.laa.2008.06.029
[15] 艾尔曼,M。;Ernst,O.G。;Güttel,S.,矩阵函数的放气重启,SIAM J.矩阵分析。申请。,32, 621-641 (2011) ·Zbl 1264.65070号 ·doi:10.1137/090774665
[16] Frommer,A。;Güttel,S。;Schweitzer,M.,基于求积的矩阵函数的高效和稳定Arnoldi重启,SIAM J.矩阵分析。申请。,35, 661-683 (2014) ·Zbl 1309.6500号 ·数字对象标识码:10.1137/13093491X
[17] T.Jawecki、W.Auzinger和O.Koch,“一类矩阵指数和(φ)-函数的Krylov近似的可计算严格上界”(2018)。https://arxiv.org/pdf/1809.03369。 ·Zbl 1432.65053号
[18] Hochbruck,M。;Lubich,C.,量子经典分子动力学的指数积分器,BIT,39,620-645(1999)·兹比尔0985.81002 ·doi:10.1023/A:1022335122807
[19] Hochbruck,M。;帕祖尔,T。;Schulz,A.,线性波动方程间断Galerkin近似的有效时间积分,Z.Angew。数学。机械。,95, 237-259 (2015) ·Zbl 1322.65095号 ·doi:10.1002/zamm.201300306
[20] 伯纳,R.-U。;Ernst,O.G。;Güttel,S.,使用有理Krylov方法的三维瞬态电磁建模,地球物理学。《国际期刊》,2022025-2043(2015)·doi:10.1093/gji/ggv224
[21] 博切夫,医学硕士。;汉斯,A.M。;Uppu,R.,用于建模单色光源多频光学响应的指数Krylov时间积分,J.Compute。申请。数学。,340, 474-485 (2018) ·Zbl 1391.65124号 ·doi:10.1016/j.cam.2017.12.014
[22] M.A.Botchev和L.A.Knizhnerman,“ART:Krylov子空间矩阵指数评估的自适应剩余时间重启”,《计算杂志》。申请。数学。364(112311)(2020).doi:10.1016/j.cam.2019.06.027·Zbl 1427.65060号
[23] J.Niehoff,博士论文(Mathematisch-Naturwissenschaftlichen Fakultät der Heinrich-Heine-Unie-Universityät Düusseldorf,2006)。
[24] Golub,G.H。;Van Loan,C.F.,《矩阵计算》(1996),巴尔的摩:约翰霍普金斯大学出版社·Zbl 0865.65009号
[25] 莫勒,C.B。;Van Loan,C.F.,《计算矩阵指数的十九种可疑方法》,25年后,SIAM Rev.,45,3-49(2003)·Zbl 1030.65029号 ·doi:10.1137/S00361445024180
[26] Higham,N.J.,《矩阵的函数:理论与计算》(2008),宾夕法尼亚州费城:SIAM,宾夕法尼亚州,费城·Zbl 1167.15001号 ·数字对象标识代码:10.1137/1.9780898717778
[27] 博切夫,医学硕士。;格林·V。;Hochbruck,M.,矩阵指数的残差、重启和Richardson迭代,SIAM J.Sci。计算。A、 351376-1397(2013)·Zbl 1278.65052号 ·数字对象标识代码:10.1137/10820191
[28] Parlett,B.N.,《对称特征值问题》(1998),费城,宾夕法尼亚州:SIAM,费城·兹伯利0885.65039 ·doi:10.1137/1.9781611971163
[29] van der Vorst,H.A.,《大型线性系统的迭代Krylov方法》(2003),剑桥:剑桥大学出版社,剑桥·Zbl 1023.65027号 ·doi:10.1017/CBO9780511615115
[30] Saad,Y.,《稀疏线性系统的迭代方法》(2003),费城,宾夕法尼亚州:SIAM,费城·Zbl 1002.65042号 ·doi:10.137/1.9780898718003
[31] Druskin,V.L。;Knizhnerman,L.A.,计算对称矩阵函数的两种多项式方法,苏联计算。数学。数学。物理。,291112-121(1989年)·兹比尔0719.65035 ·doi:10.1016/S0041-5553(89)80020-5
[32] Knizhnerman,L.A.,使用Arnoldi方法计算非对称矩阵的函数,苏联计算。数学。数学。物理。,31, 1-9 (1991) ·Zbl 0774.65021号
[33] Hochbruck,M。;Lubich,C.,《关于矩阵指数算子的Krylov子空间逼近》,SIAM J.Numer。分析。,34, 1911-1925 (1997) ·Zbl 0888.65032号 ·doi:10.1137/S0036142995280572
[34] Druskin,V.L。;格林鲍姆,A。;Knizhnerman,L.A.,《在矩阵函数计算中使用非正交Lanczos向量》,SIAM J.Sci。计算。,19, 38-54 (1998) ·Zbl 0912.65021号 ·doi:10.1137/S1064827596303661
[35] Hundsdorfer,W。;Verwer,J.G.,时间相关对流扩散反应方程的数值解(2003),柏林:施普林格出版社,柏林·Zbl 1030.65100号 ·数字对象标识代码:10.1007/978-3-662-09017-6
[36] 巴雷特·R。;贝里,M。;Chan,T.F.,《线性系统解的模板:迭代方法的构建块》(1994),宾夕法尼亚州费城:SIAM,宾夕法尼亚州,费城·doi:10.1137/1.9781611971538
[37] 萨阿德,Y。;Schultz,M.H.,GMRES:求解非对称线性系统的广义最小残差算法,SIAM J.Sci。统计计算。,7, 856-869 (1986) ·Zbl 0599.65018号 ·doi:10.1137/0907058
[38] Y.Saad,《稀疏线性系统的迭代方法》,2000年出版。www-users.cs.umn.edu/saad/books.html。
[39] Sidje,R.B.,Expokit:计算矩阵指数的软件包,ACM Trans。数学。软件,24130-156(1998)·Zbl 0917.65063号 ·数字对象标识代码:10.1145/285861.285868
[40] Krukier,L.A.,隐式差分格式和求解某类拟线性方程组的迭代方法,Sov。数学。,23, 43-55 (1979) ·Zbl 0476.65067号
[41] 塔夫罗夫,A。;Hagness,S.C.,《计算电动力学:有限差分时域方法》(2005年),马萨诸塞州波士顿:阿特奇大厦·Zbl 0963.78001号
[42] 科尔,J.S。;Figge,M.T。;De Raedt,H.,求解含时Maxwell方程的无条件稳定算法,物理学。版本E,64,066705(2001)·doi:10.1103/PhysRevE.64.066705
[43] Verwer,J.G。;Botchev,M.A.,麦克斯韦方程的无条件稳定积分,线性代数应用。,431300-317(2009年)·兹比尔1170.78004 ·doi:10.1016/j.laa.2008.12.036
[44] Botchev,M.A.,光子晶体建模中Maxwell方程的Krylov子空间指数时域解,J.Compute。申请。数学。,293, 24-30 (2016) ·Zbl 1322.65056号 ·doi:10.1016/j.cam.2015.04.022
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。