×

通过有理Krylov子空间方法加速轮廓积分技术。 (英语) Zbl 1372.65139号

摘要:我们建议对指数积分器中出现的矩阵函数和向量的乘积进行有理Krylov子空间近似。我们考虑在左复半平面扇区中具有值场的矩阵。我们的方法的极点的选择是通过基于扇区周围双曲线的轮廓积分的固定有理逼近提出的。与固定近似相比,我们的有理Krylov子空间方法显示出一个加速且更稳定的收敛阶(mathcal{O}(e^{-Cn}))。

MSC公司:

65层60 矩阵指数和相似矩阵函数的数值计算
65升04 刚性方程的数值方法
65升70 常微分方程数值方法的误差界
65平方米 含偏微分方程初值和初边值问题离散方程的数值解
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Skaflestad,B。;W.M.Wright,与指数相关的矩阵函数的缩放和修正平方法,应用。数字。数学。,59, 3-4, 783-799 (2009) ·Zbl 1160.65034号
[2] Krogstad,S.,刚性偏微分方程的广义积分因子方法,J.Compute。物理。,203, 1, 72-88 (2005) ·Zbl 1063.65097号
[3] 斯特雷梅尔,K。;Weiner,R.、Linear-implizite Runge-Kutta Methoden und ihre Anwendungen(1992)、Teubner:Teubner Leipzig·Zbl 0759.65047号
[4] Hochbruck,M。;Ostermann,A.,指数积分器,《数值学报》。,19, 209-286 (2010) ·兹比尔1242.65109
[5] 贝克曼,B。;Reichel,L.,通过Faber变换对矩阵函数进行误差估计和评估,SIAM J.Numer。分析。,47, 5, 3849-3883 (2009) ·Zbl 1204.65041号
[6] López-Fernández,M。;巴伦西亚,C。;Schädle,A.,用于反转扇形拉普拉斯变换的谱序方法,SIAM J.Numer。分析。,44, 3, 1332-1350 (2006) ·Zbl 1124.65120号
[7] 贝克曼,B。;Güttel,S.,矩阵函数逼近的有理Arnoldi方法的超线性收敛,Numer。数学。,121, 2, 205-236 (2012) ·Zbl 1271.65059号
[8] 德鲁斯金,V。;雷米斯,R。;Zaslavsky,M.,《模拟无界域中波传播的扩展Krylov子空间模型降阶技术》,J.Compute。物理。,272608-618(2014)·Zbl 1349.65545号
[9] 德鲁斯金,V。;Knizhnerman,L。;Zaslavsky,M.,使用带优化移位的有理Krylov子空间解决大规模进化问题,SIAM J.Sci。计算。,31, 5, 3760-3780 (2009) ·Zbl 1204.65042号
[10] 德鲁斯金,V。;Zaslavsky,M.,关于时不变自共轭动力系统的Krylov子空间逼近的收敛性,线性代数应用。,436, 10, 3883-3903 (2012) ·兹比尔1247.65082
[12] Frommer,A。;Simoncini,V.,厄米特矩阵和对称矩阵有理矩阵函数的停止准则,SIAM J.Sci。计算。,30, 3, 1387-1412 (2008) ·Zbl 1203.65079号
[13] Frommer,A。;Güttel,S。;Schweitzer,M.,基于求积的矩阵函数的高效和稳定Arnoldi重启,SIAM J.矩阵分析。申请。,35, 2, 661-683 (2014) ·Zbl 1309.6500号
[14] 哥克勒,T。;Grimm,V.,用于指数积分器中算子函数逼近的轴向Krylov aubspace方法的收敛性分析,SIAM J.Numer。分析。,51, 4, 2189-2213 (2013) ·兹比尔1278.65076
[15] Grimm,V.,算子函数的分解Krylov子空间近似,BIT,52,3,639-659(2012)·Zbl 1258.65052号
[16] 格林·V。;Hochbruck,M.,三角算子的有理逼近,BIT,48,2,215-229(2008)·Zbl 1148.65075号
[17] 格林·V。;Gugat,M.,用预解级数逼近半群和相关算子函数,SIAM J.Numer。分析。,48, 5, 1826-1845 (2010) ·Zbl 1222.65052号
[18] Güttel,S.,矩阵函数的有理Krylov逼近:数值方法和最优极点选择,GAMM-Mitt。,36, 1, 8-31 (2013) ·Zbl 1292.65043号
[19] 洛佩兹,L。;Simoncini,V.,有理函数逼近矩阵指数的投影方法分析,SIAM J.Numer。分析。,44, 2, 613-635 (2006) ·兹比尔1158.65031
[20] Moret,I.,关于分数进化问题的Krylov方法的注释,Numer。功能。分析。最佳。,34, 5, 539-556 (2013) ·Zbl 1272.65061号
[21] 莫雷特,I。;Novati,P.,《矩阵指数的RD-有理逼近》,BIT,44,595-615(2004)·Zbl 1075.65062号
[22] Novati,P.,使用指数积分器的限制分母有理Arnoldi方法,SIAM J.矩阵分析。申请。,32, 4, 1537-1558 (2011) ·兹比尔1247.65093
[23] 波波利齐奥,M。;Simoncini,V.,矩阵指数算子近似的加速技术,SIAM J.矩阵分析。申请。,30, 2, 657-683 (2008) ·Zbl 1168.65021号
[24] van den Eshof,J。;Hochbruck,M.,将Lanczos近似预处理为矩阵指数,SIAM J.Sci。计算。,27, 4, 1438-1457 (2006) ·Zbl 1105.65051号
[25] Hochbruck,M。;帕日尔,T。;舒尔茨,A。;Thawian,E。;Wieners,C.,线性波动方程间断Galerkin近似的有效时间积分,ZAMM Z.Angew。数学。机械。,95,3237-259(2015)·Zbl 1322.65095号
[26] Güttel,S。;Knizhnerman,L.,Cauchy-Stieltjes矩阵函数的黑盒有理Arnoldi变量,BIT,53,3,595-616(2013)·Zbl 1276.65026号
[27] Schmelzer,T。;Trefethen,L.N.,通过Carathéodory-Fejér近似和轮廓积分评估指数积分器的矩阵函数,Electron。事务处理。数字。分析。,29, 1-18 (2007-2008) ·Zbl 1186.65092号
[28] Ruhe,A.,特征值计算的有理Krylov序列方法,线性代数应用。,58, 391-405 (1984) ·Zbl 0554.65025号
[29] Ruhe,A.,非对称特征值问题的有理Krylov算法。三、 实矩阵的复数移位,BIT,34,1,165-176(1994)·Zbl 0810.65032号
[30] 洛佩斯·费尔南德斯,M。;Palencia,C.,关于某些全纯映射的拉普拉斯变换的数值反演,Appl。数字。数学。,51, 2-3, 289-303 (2004) ·Zbl 1059.65120号
[31] 黑尔,N。;新泽西州海姆。;Trefethen,L.N.,《计算》(A^\alpha,\log(A)),以及通过轮廓积分实现的相关矩阵函数,SIAM J.Numer。分析。,46, 5, 2505-2523 (2008) ·Zbl 1176.65053号
[32] Trefethen,L.N。;魏德曼,J.A.C。;Schmelzer,T.,Talbot象限和有理近似,BIT,46,3653-670(2006)·Zbl 1103.65030号
[33] Sadeghi,A。;Ismail,A.I.M.,使用梯形规则近似矩阵的第(p)个根,国际数学杂志。数学。科学。,13(2012),第634698条·Zbl 1242.65083号
[34] Stenger,F.,通过Whittaker基数函数的近似,J.近似理论,17,3,222-240(1976)·Zbl 0332.41013号
[35] Stenger,F.,基于Whittaker基数或sinc函数的数值方法,SIAM Rev.,23,2,165-224(1981)·Zbl 0461.65007号
[36] 魏德曼,J.A.C。;Trefethen,L.N.,计算Bromwich积分的抛物线和双曲线,数学。公司。,76, 259, 1341-1356 (2007) ·兹比尔1113.65119
[37] 恩格尔,K.-J。;Nagel,R.,线性发展方程的单参数半群(2000),Springer-Verlag:Springer-Verlag纽约-柏林,海德堡·Zbl 0952.47036号
[38] 哥克勒,T。;Grimm,V.,用带简单极点的有理Krylov子空间方法对指数积分器中(varphi)函数的一致逼近,SIAM J.矩阵分析。申请。,35, 4, 1467-1489 (2014) ·Zbl 1416.65192号
此参考列表基于出版商或数字数学图书馆提供的信息。它的项目与zbMATH标识符启发式匹配,并且可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。