×

具有近似最优位移的Helmholtz方程的半矩阵无网格预处理。 (英语) Zbl 1518.35247号

小结:由于它对波现象的重要性,在亥姆霍兹方程预条件的设计中付出了相当大的努力。设计预条件器的一种选择是对移位算子应用多重网格方法。在这种方法中,波数被一些虚数所移位。这一步的动机是,与标准亥姆霍兹方程相比,迭代求解器可以更有效地处理移位问题。然而,到目前为止,选择移位参数的最佳策略尚不明确。众所周知,一个好的移位参数敏感地取决于波数和离散化参数,如阶数和网格大小。因此,我们研究了如何选择一个接近最优的复数移位,以使灵活的广义最小残差(FGMRES)解算器以较少的迭代收敛。我们的目标是提供一个映射,该映射根据波数和网格大小返回预处理器的近似最佳偏移。为了计算这张图,我们考虑了一种数据驱动的方法:我们首先生成许多样本,然后在第二步中,我们对这些数据进行非线性回归。有了这个代表性的映射,通过简单的评估可以获得接近最佳的偏移。我们的预条件器基于应用于移位问题的双网格V循环,允许我们实现一种半无矩阵方法,其中只需要存储粗网格和边界矩阵。在精细网格上,只计算应用于向量的矩阵的作用,而不组装全局矩阵。这可以有效地使用计算资源,并允许在以前受可用内存限制的规模上解决问题。我们的预处理FGMRES求解器的性能通过几个二维和三维非均匀波数的基准问题进行了说明。

理学硕士:

35J05型 拉普拉斯算子、亥姆霍兹方程(约化波动方程)、泊松方程
65号55 多重网格方法;含偏微分方程边值问题的区域分解
65F08个 迭代方法的前置条件
65英尺50英寸 稀疏矩阵的计算方法
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] 阿莫科静力学测试。https://software.seg.org/datasets/2D/Statics_1994/ (1994). [在线]。2021年1月26日访问
[2] 阿尔瓦拉多,A。;Castillo,P.,LDG方法应用于多面体域中时间调和Maxwell方程的计算性能,科学杂志。计算。,67, 2, 453-474 (2016) ·兹比尔1354.78013 ·doi:10.1007/s10915-015-0087-3
[3] 阿梅斯托,P。;Buttari,A。;L'优秀,J-Y;Mary,T.,多核架构上块低阶多额因分解的性能和可扩展性,ACM Trans。数学。柔软。,45, 1-26 (2019) ·Zbl 1471.65025号 ·数字对象标识代码:10.1145/3242094
[4] 安德森,R。;安德烈·J。;巴克,A。;Bramwell,J。;卡米尔,J-S;杜布雷夫,JCV;Y.杜杜伊特。;费希尔,A。;科列夫,T。;Pazner,W。;斯托维尔,M。;托莫夫,V。;阿克曼,I。;Dahm,J。;麦地那,D。;Zampini,S.,MFEM:模块化有限元库,计算。数学。申请。,81, 42-74 (2021) ·Zbl 1524.65001号
[5] Balay,S.、Abhyankar,S.,Adams,M.F.、Brown,J.、Brune,P.、Buschelman,K.、Dalcin,L.、Dener,A.、Eijkhout,V.、Gropp,W.D.、Karpeyev,D.、Kaushik,D.、Knepley,M.G.、May,D.A.、McInnes,L.C.、Mills,R.T.、Munson,T.、Rupp,K.,Sanan,P.,Smith,B.F.、Zampini,S.和Zhang,H.:PETSc网页。https://www.mcs.anl.gov/petsc(2019年)
[6] 贝利斯,A。;哥尔斯坦,CI;Turkel,E.,亥姆霍兹方程的迭代方法,J.Compute。物理。,49, 3, 443-457 (1983) ·Zbl 0524.65068号 ·doi:10.1016/0021-991(83)90139-0
[7] 博尔顿,M。;Rittich,H.,《用多重网格法对周期模板进行傅里叶分析》,SIAM J.Sci。计算。,40、3、A1642-A1668(2018)·兹比尔1402.65175 ·doi:10.1137/16M1073959
[8] Brandt,A.,Ta'asan,S.:几乎奇异和轻微不确定问题的多重网格方法。见:《多重网格方法II》,第99-121页。施普林格(1986)·Zbl 0605.65022号
[9] Briggs,W.L.,Henson,V.E.,McCormick,S.F.:多重网格教程。SIAM(2000年)·Zbl 0958.65128号
[10] 球拍,P-H;Gander,MJ,移动的亥姆霍兹预条件器需要多大的偏移才能通过多重网格进行有效反演?,SIAM J.科学。计算。,39、2、A438-A478(2017)·Zbl 1365.65269号 ·数字对象标识码:10.1137/15M102085X
[11] 球,P-H;MJ甘德;Xiang,X.,2D常系数亥姆霍兹问题有限差分离散化的包含实位移波数的闭合形式色散校正,SIAM J.Sci。计算。,43、1、A278-A308(2021)·Zbl 1464.35064号 ·doi:10.1137/19M1304350
[12] Cools,S.公司。;Vanroose,W.,《亥姆霍兹问题的复数移位拉普拉斯预条件的局部傅里叶分析》,数值。线性代数应用。,20, 4, 575-597 (2013) ·Zbl 1313.65320号 ·doi:10.1002/nla.1881
[13] Drzisga,D。;吕德,美国。;Wohlmuth,B.,混合网格上向量值偏微分方程的Stencil缩放及其在广义牛顿流体中的应用,SIAM J.Sci。计算。,42、6、B1429-B1461(2020)·Zbl 1508.65157号 ·doi:10.1137/19M1267891
[14] 埃尔曼,HC;安永,OG;O'leary,DP,离散亥姆霍兹方程的Krylov子空间迭代增强的多重网格方法,SIAM J.Sci。计算。,23, 4, 1291-1315 (2001) ·兹比尔1004.65134 ·doi:10.1137/S1064827501357190
[15] 埃朗加,YA;奥斯特利,CW;Vuik,C.,一种新型的基于多重网格的异构亥姆霍兹问题预处理程序,SIAM J.Sci。计算。,27, 4, 1471-1492 (2006) ·Zbl 1095.65109号 ·doi:10.1137/040615195
[16] 亚拉巴马州二郎加;Vuik,C。;Oosterlee,CW,关于解亥姆霍兹方程的一类预条件,应用。数字。数学。,50, 3-4, 409-425 (2004) ·Zbl 1051.65101号 ·doi:10.1016/j.apnum.2004.01.09
[17] Ernst,O.G.,Gander,M.J.:为什么用经典迭代方法求解亥姆霍兹问题很困难。摘自:计算科学与工程课堂讲稿,第325-363页。柏林施普林格出版社(2011)·Zbl 1248.65128号
[18] Ernst,O.G.,Gander,M.J.:亥姆霍兹问题的多重网格方法:使用标准组件的一维收敛方案。在:波传播和应用中的正问题和反问题。De Gruyer,第135-186页(2013年)·Zbl 1287.65133号
[19] 埃斯特哈齐,S。;Melenk,JM,《L2和负范数中亥姆霍兹方程离散化的分析》,计算。数学。申请。,67, 4, 830-853 (2014) ·兹比尔135065123
[20] MJ甘德;格雷厄姆,IG;Spence,EA,将GMRES应用于带移位拉普拉斯预处理的亥姆霍兹方程:保证波数无关收敛的最大移位是多少?,数字。数学。,131, 3, 567-614 (2015) ·Zbl 1328.65238号 ·doi:10.1007/s00211-015-0700-2
[21] MJ甘德;Zhang,H.,Helmholtz方程的一类迭代求解器:因式分解、扫描预条件、源传输、单层势、极化道和优化Schwarz方法,SIAM Rev.,61,1,3-76(2019)·Zbl 1417.65216号 ·doi:10.1137/16M109781X
[22] 灰色,SH;Marfurt,KJ,《地形迁移:改善近地表图像》,加拿大。J.探索。地球物理学。,31, 1-2, 18-24 (1995)
[23] Greenfeld,D.,Galun,M.,Basri,R.,Yavneh,I.,Kimmel,R.:学习优化多重网格PDE解算器。摘自:《机器学习研究论文集》,第97卷,第2415-2423页,美国加利福尼亚州长滩(2019年)。PMLR公司
[24] 海因莱因,A。;Klawonn,A。;Lanser,M。;Weber,J.,《自适应区域分解方法中的机器学习——预测约束的几何位置》,SIAM J.Sci。计算。,41、6、A3887-A3912(2019)·兹比尔1429.65065 ·doi:10.1137/18M1205364
[25] Kiefer,J.,最大值的序列极小极大搜索,Proc。美国数学。《社会学杂志》,第4、3、502页(1953年)·Zbl 0050.35702号 ·doi:10.1090/S0002-9939-1953-0055639-3
[26] 克伦比克勒,M。;Ljungkvist,K.,图形处理器上无矩阵高阶有限元计算的多重网格,ACM Trans。并行计算。,6, 1, 1-32 (2019) ·doi:10.1145/3322813
[27] 库马尔,P。;罗德里戈,C。;法国加斯帕;Oosterlee,CW,《关于具有跳跃系数和随机系数的偏微分方程多重网格方法的局部傅里叶分析》,SIAM J.Sci。计算。,41、3、A1385-A1413(2019)·Zbl 1435.65053号 ·doi:10.1137/18M1173769
[28] 李,C。;乔,Z.,大型腔体电磁散射的快速预处理迭代算法,J.Sci。计算。,53, 2, 435-450 (2012) ·Zbl 1432.78012号 ·doi:10.1007/s10915-012-9580-0
[29] 刘,F。;Ying,L.,三维亥姆霍兹方程的递归扫描预条件,SIAM J.Sci。计算。,38、2、A814-A832(2016)·Zbl 1376.65035号 ·doi:10.1137/15M1010154
[30] LRZ公司。SuperMUC-NG的硬件。https://doku.lrz.de/display/PUBLIC/硬件++SuperMUC-NG,2020年2月25日访问
[31] 卢,P。;Xu,X.,基于亥姆霍兹方程区域分解方法的鲁棒多级预处理程序,J.Sci。计算。,81, 1, 291-311 (2019) ·Zbl 1442.65430号 ·doi:10.1007/s10915-019-01015-z
[32] Luz,I.、Galun,M.、Maron,H.、Basri,R.、Yavneh,I.:使用图形神经网络学习代数多重网格。摘自:机器学习国际会议,第6489-6499页。PMLR(2020年)
[33] Melenk,J.M.:关于广义有限元方法。博士论文,由数学系指导研究。马里兰大学帕克分校(1995年)
[34] Paszke,A.等人:PyTorch:一个命令式、高性能的深度学习库。收录于:Wallach,H.、Larochelle,H.,Beygelzimer,A.、Buc,F.、Fox,E.、Garnett,R.(编辑)《神经信息处理系统进展》,第32卷,第8024-8035页。Curran Associates,Inc.(2019年)
[35] Ramos,L.G.,Nabben,R.:亥姆霍兹问题的两层移位拉普拉斯预条件:值场分析和波数无关收敛。arXiv预印本arXiv:2006.08750(2020)
[36] Saad,Y.,《一种灵活的内外预处理GMRES算法》,SIAM J.Sci。计算。,14, 2, 461-469 (1993) ·Zbl 0780.65022号 ·doi:10.1137/0914028
[37] 萨阿德,Y。;Schultz,MH,GMRES:求解非对称线性系统的广义最小残差算法,SIAM J.Sci。统计计算。,856-869年7月3日(1986年)·Zbl 0599.65018号 ·doi:10.1137/0907058
[38] 斯托克,CC;艾哈迈德,M。;Bhowmik,SK,带优化粗网格校正的亥姆霍兹方程多重网格法,SIAM J.Sci。计算。,36,6,A2819-A2841(2014)·Zbl 1310.65157号 ·数字对象标识码:10.1137/13092349X
[39] 孙,Y。;Wu,S。;Xu,Y.,波动方程数值方法的并行实现,J.Sci。计算。,90, 1, 1-31 (2022) ·Zbl 1483.65056号 ·doi:10.1007/s10915-021-01701-x
[40] 陶斯,M。;Zepeda-Nüñez,L。;休伊特,RJ;Demanet,L.,L-Sweeps:高频亥姆霍兹方程的可扩展并行预处理程序,J.Compute。物理。,420 (2020) ·Zbl 07506624号 ·doi:10.1016/j.jcp.2020.109706
[41] Trottenberg,U.,Oosterlee,C.W.,Schuler,A.:多重电网。Elsevier(2000)
[42] Versteeg,R.,The Marmousi experience:在合成复杂数据集上确定速度模型,Lead。Edge,13,9,927-936(1994)·数字对象标识代码:10.1190/1.1437051
[43] Wienands,R.,Joppich,W.:多重网格方法的实用傅里叶分析。查普曼和霍尔/CRC(2004)·Zbl 1062.65133号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。