×

三角算子的有理逼近。 (英语) Zbl 1148.65075号

摘要:我们考虑用三角积分器数值求解波动方程时出现的三角算子函数的近似。众所周知,当步长大于算子的范数时,无指数衰减矩阵函数的Krylov子空间方法表现出超线性收敛性。因此,对于无界算子,Krylov近似可能无法收敛。
本文提出并分析了一种有理Krylov子空间方法,它不仅对微分算子的有限元或有限差分逼近收敛,而且对抽象无界算子也收敛。与标准Krylov方法相比,收敛性与算子的范数无关,因此与空间离散化无关。我们讨论了有限元离散化的有效实现,并用数值实验说明了我们的分析。

理学硕士:

65M60毫米 涉及偏微分方程初值和初边值问题的有限元、Rayleigh-Ritz和Galerkin方法
35升05 波动方程
65层10 线性系统的迭代数值方法
65个M12 含偏微分方程初值和初边值问题数值方法的稳定性和收敛性
PDF格式BibTeX公司 XML格式引用
全文: 内政部 链接

参考文献:

[1] L.Bergamaschi和M.Vianello,大型稀疏对称矩阵指数算子的高效计算,数值。线性代数应用。,7(2000),第27-45页·Zbl 0983.65058号 ·doi:10.1002/(SICI)1099-1506(200001/02)7:1<27::AID-NLA185>3.0.CO;2-4
[2] D.Cohen、E.Hairer和C.Lubich,通过调制傅里叶展开对非线性扰动波动方程进行长期分析,Arch。定额。机械。分析。,187(2008),第341–368页·Zbl 1145.35087号 ·doi:10.1007/s00205-007-0095-z
[3] D.Cohen、E.Hairer和C.Lubich,《非线性波动方程数值离散化中的能量守恒、动量守恒和作用》,发表于Numer。数学。(2008). ·兹比尔1163.65066
[4] D.Cohen、E.Hairer和C.Lubich,《非线性波动方程数值离散化的长时间能量守恒》,发表于《Equadiff07学报》(2008)·Zbl 1163.65066号
[5] V.A.Dougalis和S.M.Serbin,关于余弦的一类有理逼近的评论,BIT,20(1980),第204-211页·Zbl 0425.65006号 ·doi:10.1007/BF01933192
[6] V.L.Druskin和L.A.Knizhnerman,扩展Krylov子空间:矩阵平方根和相关函数的近似,SIAM J.矩阵分析。申请。,19(1998年),第755-771页·Zbl 0912.65022号 ·doi:10.1137/S0895479895292400
[7] M.Eiermann和O.Ernst,用于评估矩阵函数的重新启动的Krylov子空间方法,SIAM J.Numer。分析。,44(2006),第2481-2504页·Zbl 1129.65019号 ·doi:10.1137/050633846
[8] A.Frommer和V.Simoncini,厄米特矩阵和对称矩阵有理矩阵函数的停止准则,SIAM J.Sci。计算。,30(2008),第1387-1412页·Zbl 1203.65079号 ·数字对象标识代码:10.1137/070684598
[9] B.García-Archilla、J.Sanz-Serna和R.Skeel,振荡微分方程的长时间步长方法,SIAM J.Sci。计算。,30(3)(1998),第930-963页·Zbl 0927.65143号 ·doi:10.1137/S1064827596313851
[10] G.H.Golub和C.F.van Loan,《矩阵计算》,第二版。,约翰·霍普金斯大学出版社,马里兰州巴尔的摩,1989年·Zbl 0733.65016号
[11] V.Grimm,关于振荡二阶微分方程的Gautschi型方法的注释,Numer。数学。,102(2005),第61-66页·Zbl 1082.65080号 ·doi:10.1007/s00211-005-0639-9
[12] V.Grimm,应用于振荡二阶微分方程的Gautschi型指数积分器的误差界,数值。数学。,100(2005),第71-89页·Zbl 1074.65097号 ·doi:10.1007/s00211-005-0583-8
[13] V.Grimm,《关于波方程的高斯型指数积分器的使用》,《数值数学与高级应用》,ENUMATH 2005,A.Bermüdez de Castro,D.Gómez,P.Quintela,and P.Salgado,eds.,Springer,Berlin,Heidelberg,2006年,第557–563页·Zbl 1119.65353号
[14] V.Grimm和M.Hochbruck,振荡二阶微分方程指数积分器的误差分析,J.Phys。A、 数学。Gen.,39(2006),第5495–5507页·Zbl 1093.65078号 ·doi:10.1088/0305-4470/39/19/S10
[15] E.Hairer和C.Lubich,振荡微分方程数值方法的长时间能量守恒,SIAM J.Numer。分析。,38(2000),第414-441页·Zbl 0988.65118号 ·doi:10.1137/S00361429993594
[16] E.Hairer和C.Lubich,长期弱非线性波动方程的谱半离散化,发现。计算。数学。,8(2008),第319-334页·Zbl 1156.35007号 ·doi:10.1007/s10208-007-9014-9
[17] E.Hairer和C.Lubich,《数值哈密顿系统中的长时间振荡》,发表在艾萨克·牛顿研究所的《HOP程序学报》上,2008年·Zbl 1156.35007号
[18] M.Hochbruck和C.Lubich,关于矩阵指数算子的Krylov子空间逼近,SIAM J.Numer。分析。,34(1997),第1911–1925页·Zbl 0888.65032号 ·doi:10.1137/S0036142995280572
[19] M.Hochbruck和C.Lubich,振荡二阶微分方程的Gautschi型方法,数值。数学。,83(1999),第403-426页·Zbl 0937.65077号 ·doi:10.1007/s002110050456
[20] L.Lopez和V.Simoncini,有理函数逼近矩阵指数的投影方法分析,SIAM J.Numer。分析。,44(2)(2006),第613-635页·兹比尔1158.65031 ·数字对象标识代码:10.1137/05062590
[21] I.Moret和P.Novati,矩阵指数的RD-有理逼近,BIT,44(2004),第595-615页·Zbl 1075.65062号 ·doi:10.1023/B:BITN.0000046805.27551.3b
[22] I.Moret和P.Novati,拟核多项式零点上矩阵的插值函数,数值。线性代数应用。,11(2005),第337-353页·Zbl 1164.65387号 ·doi:10.1002/nla.393
[23] O.Nevanlinna,线性方程迭代的收敛性,Birkhäuser,巴塞尔,1993年·Zbl 0846.47008号
[24] P.P.Petrushev和V.A.Popov,实函数的有理逼近,剑桥大学出版社,1987年·Zbl 0644.41010号
[25] M.Popolizio和V.Simoncini,矩阵指数算子逼近的加速技术,SIAM J.矩阵分析。申请。,30(2008),第657–683页·Zbl 1168.65021号 ·doi:10.1137/060672856
[26] A.Ruhe,特征值计算的有理Krylov序列方法,线性代数应用。,58(1984年),第391-405页·Zbl 0554.65025号 ·doi:10.1016/0024-3795(84)90221-0
[27] Y.Saad,矩阵指数算子的Krylov子空间逼近分析,SIAM J.Numer。分析。,29(1992),第209-228页·Zbl 0749.65030号 ·doi:10.1137/0729014
[28] J.M.Sanz-Serna,高振荡微分方程的Mollified脉冲方法,SIAM J.Numer。分析。,46(2008),第1040–1059页·Zbl 1169.65081号 ·数字对象标识代码:10.1137/070681636
[29] G.R.Sell和Y.You,进化方程动力学,Springer,纽约,柏林,海德堡,2002年·Zbl 1254.37002号
[30] A.F.Timan,实变量函数逼近理论,佩加蒙出版社,牛津,1963年·Zbl 0117.2901号
[31] J.van den Eshof和M.Hochbruck,矩阵指数的Lanczos近似预处理,SIAM J.Sci。计算。,24(4)(2006),第1438-1457页·Zbl 1105.65051号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。