×

子域通缩与局部AMG相结合:使用AMGCL库的案例研究。 (英语) Zbl 1451.65053号

摘要:本文提出了子域收缩方法和局部代数多重网格的组合,作为一种可扩展的分布式内存预处理程序,能够求解大型线性方程组。该算法的实现作为开源AMGCL库的一部分提供给社区。该解决方案针对同构(仅CPU)和异构(CPU/GPU)系统,前者采用混合MPI/OpenMP方法,后者采用MPI、OpenMP和CUDA的组合。OpenMP的使用最大限度地减少了MPI进程的数量,从而减少了通缩方法的通信开销,并提高了预处理程序的弱可伸缩性和强可伸缩性。为了评估所实现算法的性能,考虑了标量(每个网格节点的单自由度)、类泊松系统以及源自Navier-Stokes方程离散化的非标量问题的示例。与基于成熟的Trilinos ML包的传统全局AMG预处理器进行了比较。

MSC公司:

65层50 稀疏矩阵的计算方法
65层10 线性系统的迭代数值方法
65日元10 特定类别建筑的数值算法
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] A.Alexandrescu,《现代(C++)设计:应用的通用编程和设计模式》(Addison-Wesley,Reading,MA,2001)。
[2] 巴雷特·R。;贝里,M。;Chan,T.F。;德梅尔,J。;J.多纳托。;Dongarra,J。;艾伊霍特,V。;波佐,R。;罗明,Ch。;van der Vorst,H.,《线性系统解的模板:迭代方法的构建块》(1994),费城:SIAM,费城
[3] Brandt,A。;McCoruick,S。;Huge,J.,《稀疏性及其应用》(1985),剑桥:剑桥大学出版社,剑桥·Zbl 0544.00021号
[4] O·布洛克。;Grote,M.J.,《几何和代数多重网格的稀疏近似逆平滑器》,应用。数字。数学。,41, 61-80 (2002) ·Zbl 0995.65129号 ·doi:10.1016/S0168-9274(01)00110-6
[5] 克利里·A·J。;R·D·法尔古特。;范·埃姆登,H。;Jones,J.E。;Manteuffel,Th.A。;圣·F·麦考密克。;米兰达,G.N。;Ruge,J.W.,代数多重网格的健壮性和可扩展性,SIAM J.Sci。计算。,21, 1886-1908 (2000) ·Zbl 0959.65049号 ·doi:10.1137/S1064827598339402
[6] Dadvand,P。;罗西,R。;吉尔,M。;Martorell,X。;Cotela,J。;Juanpere,E。;Idelsohn,S.R。;Oñate,E.,通用多物理框架到HPC环境的迁移,计算。流体,80,301-309(2013)·Zbl 1426.76644号 ·doi:10.1016/j.compfluid.2012.02.004
[7] Dadvand,P。;罗西,R。;Oñate,E.,为多学科应用开发有限元代码的面向对象环境,Arch。计算。《工程方法》,17253-297(2010)·Zbl 1360.76130号 ·doi:10.1007/s11831-010-9045-2
[8] D.Demidov,“MGCL:用代数多重网格方法求解大型稀疏线性系统的C++库”(2017)。
[9] Demidov,D.,AMGCL:一种高效、灵活和可扩展的代数多重网格实现,Lobachevskii J.Math。,40, 535-546 (2019) ·兹比尔1452.65426 ·doi:10.1134/S1995080219050056
[10] J.Donéa。;Huerta,A.,《流动问题的有限元方法》(2003),纽约:威利出版社,纽约
[11] Elman,H.C.,《不可压缩流模型的预处理策略》,《科学杂志》。计算。,25, 347-366 (2005) ·Zbl 1203.76098号 ·doi:10.1007/s10915-004-4648-0
[12] Frank,J。;Vuik,C.,《基于通缩的预条件的构建》,SIAM J.Sci。计算。,23, 442-462 (2001) ·Zbl 0997.65072号 ·doi:10.1137/S1064827500373231
[13] Erich Gamma,《设计模式:可重用面向对象软件的要素》(Pearson Education,印度,1995年)·Zbl 0814.68050号
[14] M.W.Gee、C.M.Siefert、J.J.Hu、R.S.Tuminaro和M.G.Sala,“ML 5.0平滑聚合用户指南”,技术代表SAND2006-2649(Sandia Natl.Labor.,2006)。
[15] 格梅纳,B。;Huber,M。;约翰·L。;吕德,美国。;Wohlmuth,B.,《斯托克斯解算器在极端规模下的定量性能研究》,J.Compute。科学。,17, 509-521 (2016) ·doi:10.1016/j.jocs.2016.06.006
[16] Hénon,P。;Ramet,P。;Roman,J.,PASTIX:稀疏对称正定系统的高性能并行直接求解器,并行计算。,28, 301-321 (2002) ·Zbl 0984.68208号 ·doi:10.1016/S0167-8191(01)00141-7
[17] 霍格,J。;Scott,J.,多核架构的新型并行稀疏直接求解器,算法,6702-725(2013)·Zbl 1461.65063号 ·doi:10.3390/a6040702
[18] 林,P。;贝当古,M。;多米诺,S。;费希尔,T。;霍姆曼,M。;胡,J。;菲普斯,E。;普罗科彭科,A。;拉贾马尼卡姆,S。;Siefert,Ch.,《利用第二代trilinos实现低马赫数流体的极值模拟》,《并行程序》。莱特。,24, 1442005 (2014) ·doi:10.1142/S0129626414420055
[19] 马修斯,J。;Walker,R.L.,《物理数学方法》(1970),纽约:W.A.Benjamin,纽约
[20] 纳本,R。;Vuik,C.,通货紧缩与平衡预条件的比较,SIAM J.Sci。计算。,27, 1742-1759 (2006) ·Zbl 1105.65049号 ·电话:10.1137/040608246
[21] Ruge,J.W。;Stüben,K.,多重网格方法(1987),费城:SIAM,费城·兹伯利0659.65094
[22] Yousef Saad,稀疏线性系统的迭代方法(SIAM,费城,2003)·Zbl 1031.65046号
[23] B.Schäling,《Boost(C++)Libraries》(XML出版社,加利福尼亚州拉古纳山,2014年)。
[24] Sleijpen,G.L.G。;Fokkema,D.R.,BiCGstab(l),涉及具有复谱的非对称矩阵的线性方程组,电子。事务处理。数字。分析。,1, 11-32 (1993) ·Zbl 0820.65016号
[25] B.史密斯。;比约斯塔德,P。;Gropp,W.,《区域分解:椭圆偏微分方程的并行多层方法》(2004),剑桥:剑桥大学出版社,剑桥
[26] Stuben,K.,“代数多重网格(AMG):应用简介”,GMD报告70(1999),德国:GMD,德国圣奥古斯丁
[27] Tang,J.M。;纳本,R。;Vuik,C。;Erlangga,Y.A.,从通缩、区域分解和多重网格方法导出的两级预条件的比较,SIAM J.Sci。计算。,39, 340-370 (2009) ·Zbl 1203.65073号 ·doi:10.1007/s10915-009-9272-6
[28] J.M.Tang和C.Vuik,泡沫流动问题通货紧缩技术的新变体(代尔夫特科技大学,代尔夫脱,2006)·Zbl 1235.76043号
[29] 美国特罗滕贝格。;Oosterlee,C。;Schüller,A.,Multigrid(2001),伦敦:学术出版社,伦敦·Zbl 0976.65106号
[30] Stefan Turek,《不可压缩流问题的高效求解器:算法和计算方法》(Springer,Berlin,Heidelberg,1999)·Zbl 0930.76002号
[31] Vuik,C。;西格尔,A。;Meijerink,J.A.,一种有效的预处理CG方法,用于解决系数差异极大的一类分层问题,J.Compute。物理。,152, 385-403 (1999) ·Zbl 0945.76048号 ·doi:10.1006/jcph.1999.6255
[32] Meier Yang,U.,BoomerAMG:并行代数多重网格解算器和预条件器,应用。数字。数学。,41, 155-177 (2002) ·Zbl 0995.65128号 ·doi:10.1016/S0168-9274(01)00115-5
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。