×

心脏力学可分级多级预处理的比较研究。 (英语) Zbl 07742907号

摘要:在这项工作中,我们在结构化和非结构化有限元网格上对心脏力学的平衡域约束分解(BDDC)和代数多重网格(AMG)预处理程序进行了性能比较。心肌的力学行为可以用三维有限弹性方程来描述,这些方程在空间上用有限元离散,得到一个大规模非线性代数系统的解。该问题用Newton-Krylov方法求解,其中Jacobian线性系统的解由BDDC/AMG预条件加速。我们深入研究了BDDC预处理器的主要参数,以使比较公平。我们重点关注:BDDC算法的局部问题和粗问题的不同直接求解器的性能;BDDC原始自由度不同选择的影响;以及有限元度的影响。可伸缩性测试是在Linux集群上执行的,最多1024个处理器,最后我们对一个真实的机电模拟进行了性能研究。

MSC公司:

65传真 数值线性代数
65牛顿 偏微分方程边值问题的数值方法
92Cxx码 生理、细胞和医学主题
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] Ambrosi,D。;阿里奥利,G。;Nobile,F。;Quarteroni,A.M.,《心脏动力学中的机电耦合:主动应变方法》,SIAM J.Appl。数学。,71005-621(2011年)·Zbl 1419.74174号
[2] Adams,M.,《固体力学中三维有限元问题的三种非结构化多重网格方法的评估》,国际J。数值。方法工程,55,5,519-534(2002)·Zbl 1076.74547号
[3] 中华人民共和国埃姆斯泰。;达夫,I.S。;L'Excellent,J.-Y。;Koster,J.,Mumps:一种通用分布式内存稀疏求解器,(国际应用并行计算研讨会(2000),施普林格),121-130
[4] Africa,P.C.,《复杂有限元问题数值解的灵活、高性能库》,SoftwareX,20,第101252条,pp.(2022)
[5] 奥古斯丁,C.M。;奈克,A。;利伯曼,M。;Prassl,A.J。;尼德勒,S.A。;Haase,G。;普朗克,G.,《人体整个心脏机电学的解剖学精确高分辨率建模:非线性变形的强可伸缩代数多重网格求解器方法》,J.Compute。物理。,305, 622-646 (2016) ·Zbl 1349.92032号
[6] Ambrosetti,A。;Prodi,G.,《非线性分析入门》,第34卷(1995),剑桥大学出版社·Zbl 0818.47059号
[7] 巴莱,S。;Abhyankar,S。;M.F.亚当斯。;Brown,J。;布鲁纳,P。;Buschelman,K。;达尔星。;Dener,A。;埃伊霍特,V。;格罗普,W。;卡佩耶夫,D。;考希克,D。;Knepley,M。;D.梅。;Curfman McInnes,L。;米尔斯,R。;Munson,T。;鲁普,K。;萨南,P。;B.史密斯。;赞皮尼,S。;张,H。;Zhang,H.(2021),阿贡国家实验室,PETSc用户手册。技术报告ANL-95/11-3.13版
[8] 拜耳,J.D。;R.C.布莱克。;普朗克,G。;Trayanova,N.A.,一种新的基于规则的算法,用于为计算心脏模型分配心肌纤维方向,Ann.Biomed。工程师,40,10,2243-2254(2012)
[9] 班格思,W。;哈特曼,R。;Kanschat,G.,交易。ii-通用面向对象有限元库,ACM Trans。数学。软质。,33、4(2007)、24秒·Zbl 1365.65248号
[10] 贝克,A.H。;Klawonn,A。;科列夫,T。;Lanser,M。;莱茵巴赫,O。;Yang,U.M.,弹性为50万并行任务的经典代数多重网格的可伸缩性(2016),Springer
[11] 布雷齐纳,M。;Tong,C。;Becker,R.,《结构力学的并行代数多重网格》,SIAM J.Sci。计算。,27, 5, 1534-1554 (2006) ·Zbl 1097.74061号
[12] 卡斯滕森,C。;Dolzmann,G.,小载荷下多凸面材料非线性弹性有限元离散的先验误差估计,Numer。数学。,97, 1, 67-80 (2004) ·Zbl 1057.74038号
[13] Colli Franzone,P。;Pavarino,L.F。;Scacchi,S.(《心脏电生理数学》,第13卷(2014年),施普林格出版社)·Zbl 1318.92002号
[14] Colli Franzone,P。;Pavarino,L.F。;Scacchi,S.,HPC架构上可扩展心脏机电解算器的数值研究,前沿。生理学。,9, 268 (2018)
[15] 骑士,C。;Pellegrini,F.,PT Scotch:一种高效并行图排序工具,并行计算。,34, 6-8, 318-331 (2008)
[16] da Veiga,L.B。;Pavarino,L.F。;Scacchi,S。;Widlund,O.B。;Zampini,S.,具有高级缩放的等几何BDDC预处理器,SIAM J.Sci。计算。,36、3、A1118-A1139(2014)·Zbl 1320.65047号
[17] Davis,T.A.,UMFPACK 5.2.0版用户指南(2007年),佛罗里达大学,25
[18] Davis,T.A。;Natarajan,E.P.,算法907:KLU,电路仿真问题的直接稀疏求解器,ACM Trans。数学。软质。,37, 3, 1-17 (2010) ·Zbl 1364.65066号
[19] Dohrmann,C.R.,基于约束能量最小化的子结构预条件,SIAM J.Sci。计算。,25, 1, 246-258 (2003) ·Zbl 1038.65039号
[20] De Sterck,H。;R·D·法尔古特。;Nolting,J.W。;Yang,U.M.,并行代数多重网格的距离二插值,数值。线性代数应用。,15, 2-3, 115-139 (2008) ·Zbl 1212.65139号
[21] De Sterck,H。;杨女士。;Heys,J.J.,降低并行代数多重网格预条件的复杂性,SIAM J.矩阵分析。申请。,27, 4, 1019-1039 (2006) ·Zbl 1102.65034号
[22] El Maliki,A。;Fortin,M。;北卡罗来纳州塔迪厄。;Fortin,A.,《三维线性和非线性弹性问题的迭代求解器:位移和混合公式》,国际期刊Numer。方法工程,83,13,1780-1802(2010)·Zbl 1202.74005号
[23] Farhat,C。;Lesoinne,M。;LeTallec,P。;Pierson,K。;Rixen,D.,FETI-DP:双重原始统一FETI方法第一部分:两级FETI方法的更快替代方案,Int.J.Numer。方法工程,50,7,1523-1544(2001)·Zbl 1008.74076号
[24] 弗朗西斯基尼,A。;Paduletto Magri,V.A。;Mazzucco,G。;斯皮齐亚,N。;Janna,C.,结构力学的鲁棒自适应代数多重网格线性求解器,计算。方法应用。机械。工程,352389-416(2019)·Zbl 1441.74302号
[25] R·D·法尔古特。;Yang,U.M.,Hypre:一个高性能预处理器库,(国际计算科学会议(2002),施普林格),632-641·Zbl 1056.65046号
[26] 加布里埃尔,E。;Fagg,G.E。;博西尔卡,G。;Angskun,T。;Dongarra,J.J。;斯奎尔斯,J.M。;Vishal,S。;坎巴杜尔,P。;巴雷特,B。;Lumsdaine,A。;卡斯坦,R.H。;Daniel,D.J。;格雷厄姆,R.L。;Woodall,T.S.,《开放式MPI:下一代MPI实现的目标、概念和设计》,(第11届欧洲PVM/MPI用户小组会议论文集,第11届欧盟PVM/MP1用户小组会议,匈牙利布达佩斯(2004年9月)),97-104
[27] 古奇奥内,J.M。;McCulloch,A.D。;Waldman,L.K.,《根据圆柱形模型测定完整心室心肌的被动材料特性》,J.Biomech。工程,113,1,42-55(1991)
[28] 格里贝尔,M。;Oeltz,D。;Schweitzer,M.A.,线性弹性的代数多重网格法,SIAM J.Sci。计算。,25, 2, 385-407 (2003) ·Zbl 1163.65336号
[29] Holzapfel,G.A.,《非线性固体力学:工程科学的连续方法》,麦加尼卡,37,4,489-490(2002)
[30] Heroux,医学硕士。;Willenbring,J.M.,《trilinos项目的新概述》,科学。程序。,20, 2, 83-88 (2012)
[31] 亨森,V.E。;杨。Boomeramg,U.M.,并行代数多重网格解算器和预条件器,应用。数字。数学。,41, 1, 155-177 (2002) ·Zbl 0995.65128号
[32] 姜瑜。;陈,R。;Cai,X.C.,非结构化网格上模拟左心室的高度并行隐式区域分解方法,计算。机械。,66, 6, 1461-1475 (2020) ·Zbl 1468.74066号
[33] Jones,麻省理工学院。;Plassmann,P.E.,《平行图着色启发式》,SIAM J.Sci。计算。,14, 3, 654-669 (1993) ·兹比尔0772.68046
[34] Kwack,J。;Bauer,G.H。;Koric,S.,《Blue Waters-Cray XE6/XK7系统上并行线性方程求解器的性能测试》,(Cray用户小组会议(CUG2016)(2016))
[35] 卡拉贝拉斯,E。;格塞尔,M.A.F。;Haase,G。;普朗克,G。;Augustin,C.M.,一种精确、稳健、高效的有限元框架,适用于各向异性、几乎不可压缩和完全不可压缩的弹性,Comput。方法应用。机械。工程,394,第114887条pp.(2022)·Zbl 1507.74483号
[36] Karypis,G。;Schloegel,K。;Kumar,V.,Parmetis:并行图划分和稀疏矩阵排序库(1997)
[37] Klawonn,A。;莱因巴赫,O。;Widlund,O.B.,《平面中不规则子域上FETI-DP算法的分析》,SIAM J.Numer。分析。,46, 5, 2484-2504 (2008) ·Zbl 1176.65135号
[38] Levick,J.R.,《心血管生理学导论》(2013),巴特沃斯·海尼曼
[39] Land,S。;古列夫。;Arens,S。;奥古斯丁,C.M。;Baron,L。;布莱克·R。;布拉德利,C。;卡斯特罗,S。;克罗齐(A.Crozier)。;法维诺,M。;Fastl,T.E。;弗里茨,T。;高,H。;Gizzi,A。;格里菲斯,B.E。;Hurtado,D.E。;克劳斯,R。;罗,X。;纳什,M.P。;Pezzuto,S。;普朗克,G。;罗西,S。;Ruprecht,D。;Seemann,G.等人。;N.P.史密斯。;Sundnes,J。;赖斯,J.J。;特拉扬诺娃,N。;王,D。;王振杰。;Niederer,S.A.,《心脏力学软件验证:测试主动和被动材料行为的基准问题和解决方案》,Proc。R.Soc.A,数学。物理学。工程科学。,4712184(2015年12月),允许不受限制地使用,前提是原始作者和来源被记入贷方
[40] Li,X.S.,《超级lu概述:算法、实现和用户界面》,ACM Trans。数学。软质。,31, 3, 302-325 (2005) ·Zbl 1136.65312号
[41] Luby,M.,最大独立集问题的简单并行算法,(第十七届ACM计算理论年会论文集(1985)),1-10
[42] 李,J。;Widlund,O.B.,FETI-DP,BDDC和block Cholesky方法,国际期刊数字。方法工程,66,2,250-271(2006)·Zbl 1114.65142号
[43] 米切尔,L。;Bishop,M。;Hötzl,E。;奈克,A。;利伯曼,M。;Haase,G。;普朗克,G.,《peta-flops计算时代器官水平的心脏电生理建模》(AIP会议论文集,第1281卷(2010),美国物理研究所),407-410
[44] 曼德尔,J。;Sousedík,B。;Dohrmann,C.R.,多空间和多级BDDC,计算,83,2-3,55-85(2008)·Zbl 1163.65091号
[45] 曼德尔,J。;Sousedík,B。;Simístek,J.,三维自适应bddc,数学。计算。模拟。,82, 10, 1812-1831 (2012) ·Zbl 1255.65225号
[46] 佩奇斯坦,C。;Dohrmann,C.R.,自适应BDDC的统一框架,Electron。变速器。数字。分析。,46273-336,3(2017)·Zbl 1368.65043号
[47] Pfaller,M.R。;赫尔曼,J.M。;Weigl,M。;Nagler,A。;Chabiniok,R。;Bertoglio,C。;Wall,W.A.,《心包对心脏生物力学的重要性:从生理学到计算建模》,《生物医学》。模型。机械双醇。,18, 2, 503-529 (2019)
[48] 帕瓦里诺,L.F。;Scacchi,S。;Zampini,S.,非线性心脏力学的Newton-Krylov-BDDC解算器,计算。方法应用。机械。工程,295,562-580(2015)·Zbl 1423.74913号
[49] Quarteroni,A.M。;Lassila,T。;罗西,S。;Ruiz-Baier,R.,用于模拟心脏功能的心电耦合多尺度和多物理模型,计算。方法应用。机械。工程,314,345-407(2017)·Zbl 1439.74208号
[50] Quarteroni,A.M。;Valli,A.,偏微分方程的数值逼近,第23卷(2008),Springer Science&Business Media·Zbl 1151.65339号
[51] 雷加佐尼,F。;萨尔瓦多,M。;中国,非洲。;费德勒,M。;Dedè,L。;Quarteroni,A.M.,《闭环血液循环的心脏机电模型与集总参数模型耦合》,J.Compute。物理。,457,第111083条pp.(2022)·Zbl 1515.92032号
[52] Saad,Y.,《稀疏线性系统的迭代方法》(2003),SIAM·Zbl 1002.65042号
[53] Smith,B.F.,线性弹性问题有限元解的最佳区域分解预处理程序,SIAM J.Sci。统计计算。,13, 1, 364-378 (1992) ·Zbl 0751.73059号
[54] N.P.史密斯。;Nickerson,D.P。;克兰平,E.J。;Hunter,P.J.,心脏的多尺度计算模型,《数值学报》。,13, 371-431 (2004) ·兹比尔1112.92021
[55] Stallman,R.M.,《使用GNU编译器集合》(4(02)2003),自由软件基金会
[56] Stüben,K.,代数多重网格综述,(数值分析:20世纪的历史发展(2001)),331-359·Zbl 0979.65111号
[57] 托塞利,A。;Widlund,O.,《区域分解方法——算法和理论》,第34卷(2004年),Springer Science&Business Media
[58] Usyk,T.P。;LeGrice,I.J。;McCulloch,A.D.,三维心脏机电力学计算模型,计算机。视觉。科学。,4、4、249-257(2002年7月)·Zbl 1001.92005号
[59] 赖特,S。;Nocedal,J.,《数值优化》,第35卷(67-68)(1999),施普林格科学·Zbl 0930.65067号
[60] 徐,J。;Zikatanov,L.,代数多重网格方法,《数值学报》。,26, 591-721 (2017) ·Zbl 1378.65182号
[61] Yang,U.M.,并行代数多重网格方法-高性能预条件,(Bruaset,Are Magnus;Tveito,Aslak,并行计算机上偏微分方程的数值解(2006),Springer Berlin Heidelberg:Springer Bolin Heitelberg Berlin,Heidelbrg),209-236·Zbl 1097.65125号
[62] Zampini,S.,PCBDDC:PETSc中的一类稳健的双精度方法,SIAM J.Sci。计算。,38、5、S282-S306(2016)·Zbl 1352.65632号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。