皮埃尔·乔利韦;罗曼（Jose E.Roman）。;斯特凡诺·扎皮尼 KSPHPDDM和PCHPDDM：使用先进的Krylov方法和稳健的多层重叠Schwarz预条件扩展PETSc。 （英语） Zbl 1524.65135号 计算。数学。申请。 84, 277-295 (2021). 概要：计算科学和工程的当代应用通常需要求解可能具有不同大小、形状和结构的线性系统。本文的目标是解释PETSc和HPDDM这两个库是如何接口的，以便为最终用户提供强大的重叠Schwarz预处理程序和先进的Krylov方法，这些方法具有回收和处理多个右手边的能力。实现的灵活性通过最简、易于运行和可复制的示例展示和解释，以便于将这些算法集成到更高级的框架中。提供的示例涵盖了特征分析、弹性、燃烧和电磁学的应用。 引用于10文件 理学硕士： 65层10 线性系统的迭代数值方法 65号55 多重网格方法；含偏微分方程边值问题的区域分解 2005年5月 并行数值计算 关键词：Krylov方法;区域分解预处理器;分布式内存并行计算 软件：github;FROSch公司;拜洛什;货币基金组织;PCBDDC公司;xSDK软件包;洛佩克。米;TetGen公司;阿美索斯2;超级LU DIST;BDDC公司;超级LU;自由Fem++;炒作;元素;UMFPACK公司;SLEPc公司;PETSc公司;特里利诺斯;流行性腮腺炎 PDF格式BibTeX公司 XML格式引用 \textit{P.Jolivet}等人，《计算》。数学。申请。84、277--295（2021；Zbl 1524.65135） 全文： 内政部 参考文献： [1] Bartlett，R。；德梅什科，I。；甘布林，T。；哈蒙德，G。；Heroux，医学硕士。；约翰逊，J。；Klinvex，A。；Li，X.S。；McInnes，L.C.公司。；莫尔顿，J.D。；Osei-Kuffuor，D。；萨里奇，J。；B.F.史密斯。；Willenbring，J。；Yang，U.M.，XSDK基金会：走向极端规模的科学软件开发工具包，超级计算机。前面。因诺夫。，2017年4月1日，URLhttps://xsdk.info网址 [2] 巴拉伊，S。；Abhyankar，S。；M.F.亚当斯。；Brown，J。；布鲁纳，P。；Buschelman，K。；达尔星。；Dener，A。；埃伊霍特，V。；格罗普，W.D。；卡佩耶夫，D。；考希克，D。；Knepley，M.G。；May，D.A。；McInnes，L.C.公司。；Mills，R.T。；Munson，T。；鲁普，K。；萨南，P。；B.F.史密斯。；扎皮尼，S。；张，H。；Zhang，H.，PETSc网页（2020年），URL网址：http://www.mcs.anl.gov/petsc [3] 巴拉伊，S。；Abhyankar，S。；M.F.亚当斯。；Brown，J。；布鲁纳，P。；Buschelman，K。；大蒜素。；Dener，A。；埃伊霍特，V。；格罗普，W.D。；卡佩耶夫，D。；考希克，D。；Knepley，M.G。；May，D.A。；McInnes，L.C.公司。；Mills，R.T。；Munson，T。；鲁普，K。；萨南，P。；B.F.史密斯。；扎皮尼，S。；张，H。；Zhang，H.，PETSC用户手册技术。代表ANL-95/11-3.13版（2020年），阿贡国家实验室 [4] 法尔古特，R。；Yang，U.M.，Hypre:高性能预处理程序库，（计算科学-ICCS，2002（2002）），632-641，URLhttps://www.llnl.gov/casc/hypre ·Zbl 1056.65046号 [5] Si，H.，TetGen：优质四面体网格生成器和3D DElaunay三角仪技术。代表，13（2013），URLhttp://wias-berlin.de/software/tetgen [6] 伯斯特德，C。；Wilcox，L.C.公司。；Ghattas，O.，：八叉树森林上并行自适应网格细化的可缩放算法，SIAM J.Sci。计算。，33、3、1103-1133（2011），网址网址：http://www.p4est.org ·Zbl 1230.65106号 [7] Jolivet，P。；Hecht，F。；F.纳塔夫。；Prud’homme，C.，异质椭圆问题的可伸缩区域分解预条件，（高性能计算、网络、存储和分析国际会议论文集，SC13（2013），ACM） [8] Jolivet，P。；Tournier，P.-H.，《块迭代方法和循环以改进线性求解器的可扩展性》（2016年高性能计算、网络、存储和分析国际会议论文集，SC16（2016），IEEE） [9] T.Hoefler，R.Belli，《并行计算系统的科学基准测试：报告性能结果时告知大众的十二种方法》，载于《2015年高性能计算、网络、存储和分析国际会议论文集》，SC15，2015年。 [10] 埃尔南德斯，V。；罗曼，J.E。；Vidal，V.，SLEPc：用于解决特征值问题的可扩展且灵活的工具包，ACM Trans。数学。软件，31，3，351-362（2005），URL网址：https://slepc.upv.es ·Zbl 1136.65315号 [11] 巴雷特·R。；贝里，M。；Chan，T.F。；德梅尔，J。；J.多纳托。；Dongarra，J。；埃伊霍特，V。；波佐，R。；罗明，C。；Van der Vorst，H.，《线性系统解的模板：迭代方法的构建块》（1994），SIAM [12] Saad，Y.，《稀疏线性系统的迭代方法》（2003），SIAM·Zbl 1002.65042号 [13] 萨阿德，Y。；Schultz，M.H.，GMRES:求解非对称线性系统的广义最小残差算法，SIAM J.Sci。统计计算。，856-869年7月3日（1986年）·Zbl 0599.65018号 [14] Hestenes，M.R。；Stiefel，E.，求解线性系统的共轭梯度方法，J.Res.Nat.Bureau Standards，49，6，409-436（1952）·Zbl 0048.09901号 [15] Van der Vorst，H.A.，Bi-CGSTAB：非对称线性系统解的Bi-CG的一种快速且平滑收敛的变体，SIAM J.Sci。统计计算。，13, 2, 631-644 (1992) ·Zbl 0761.65023号 [16] 贝克，A.H。；Jessup，E.R。；Manteuffel，T.，《加速重新启动的GMRES收敛的技术》，SIAM J.矩阵分析。申请。，26, 4, 962-984 (2005) ·Zbl 1086.65025号 [17] 埃赫尔，J。；Burrage，K。；Pohl，B.，通货紧缩预处理的重新启动的GMRES，J.计算。申请。数学。，69,2303-318（1996年）·Zbl 0854.65025号 [18] 华盛顿州Wakam。；Pacull，F.，由可压缩流产生的线性系统的内存高效混合代数解算器，计算与流体，80，158-167（2013）·Zbl 1284.76259号 [19] Saad，Y.，《一种灵活的内外预处理GMRES算法》，SIAM J.Sci。计算。，14, 2, 461-469 (1993) ·Zbl 0780.65022号 [20] 艾森斯塔特，S.C。；Elman，H.C。；Schultz，M.H.，非对称线性方程组的变分迭代方法，SIAM J.Numer。分析。，20, 2, 345-357 (1983) ·Zbl 0524.65019号 [21] Parks，M.L。；de Sturler，E。；Mackey，G。；约翰逊，D.D。；Maiti，S.，回收线性系统序列的Krylov子空间，SIAM J.Sci。计算。，28, 5, 1651-1674 (2006) ·Zbl 1123.65022号 [22] Gutknecht，M.H.，《具有多个右手边的线性系统的Block krylov空间方法：简介》，（Siddiqui，A.；Duff，I.；Christensen，O.，《真实世界系统的现代数学模型、方法和算法》（2006），420-447 [23] Frommer，A。；Lund，K。；Szyld，D.B.，矩阵函数的块krylov子空间方法，电子。事务处理。数字。分析。，47, 100-126 (2017) ·Zbl 1372.65138号 [24] 托尼尔，P.-H。；Aliferis，I。；博纳佐利，M。；de Buhan，M。；Darbas，M。；多林，V。；Hecht，F。；Jolivet，P。；El Kanfoud，I。；米利亚乔，C。；F.纳塔夫。；皮肖特，C。；Semenov，S.，使用高性能计算对脑血管意外进行微波断层成像，并行计算。，85, 88-97 (2019) [25] 卡兰齐斯，V。；马洛西，A.C.I。；Bekas，C。；A.库里奥尼。；加洛普洛斯，E。；Saad，Y.，《数据分析中多个右侧的可扩展迭代密集线性系统解算器》，并行计算。，74136-153（2018） [26] Knyazev，A.V.，走向最优预条件本征解算器：局部最优块预条件共轭梯度法，SIAM J.Sci。计算。，23, 2, 517-541 (2001) ·Zbl 0992.65028号 [27] 卡兰德拉，H。；格拉顿，S。；Langou，J。；Pinel，X。；Vasseur，X.，块重启GMRES方法的灵活变体及其在地球物理中的应用，SIAM J.Sci。计算。，34、2、A714-A736（2012）·兹比尔1248.65036 [28] 樱井，T。；塔达诺，H。；Kuramashi，Y.，块krylov子空间算法在格子QCD中具有多个右侧的wilson-Dirac方程中的应用，计算。物理学。Comm.，181，1，113-117（2010）·Zbl 1205.81131号 [29] Heroux，医学硕士。；Bartlett，R.A。；豪尔，V.E。；霍克斯特拉，R.J。；胡建杰。；科尔达·T·G。；Lehoucq，R.B。；Long，K.R。；Pawlowski，R.P。；Phipps，E.T.，《trilinos项目概述》，ACM Trans。数学。柔和。（TOMS），31，3，397-423（2005），网址https://trilinos.github.io ·Zbl 1136.65354号 [30] Bavier，E。；霍姆曼，M。；Rajamanickam，S。；Thornquist，H.、Amesos2和Belos:大型稀疏线性系统的直接和迭代求解器，科学。程序。，20, 3, 241-255 (2012) [31] Notay，Y.，《柔性共轭梯度》，SIAM J.Sci。计算。，22, 4, 1444-1460 (2000) ·Zbl 0980.65030号 [32] 卡瓦略，L.M。；格拉顿，S。；拉戈，R。；Vasseur，X.，《带内正交化和收缩重启的柔性广义共轭残差法》，SIAM J.Matrix Ana。申请。，32, 4, 1212-1235 (2011) ·兹比尔1244.65047 [33] O'Leary，D.P.，块共轭梯度算法及相关方法，线性代数应用。，29, 293-322 (1980) ·Zbl 0426.65011号 [34] 纪浩。；Li，Y.，一种无故障块共轭梯度法，BIT-Numer。数学。，57, 2, 379-403 (2017) ·Zbl 1392.65057号 [35] 萨阿德，Y。；杨，M。；埃赫尔，J。；Guyomarc'h，F.，共轭梯度算法的缩减版本，SIAM J.Sci。计算。，21, 5, 1909-1926 (2000) ·Zbl 0955.65021号 [36] Stathopoulos，A。；Abdel-Rehim，A。；Orginos，K.，QCD中具有多个右手边的反演的通货紧缩，（物理学杂志：会议系列，第180卷（2009），IOP出版） [37] 索达尔特，K.M。；Szyld，D.B。；Xue，F.，移位线性系统序列的Krylov子空间循环，应用。数字。数学。，81, 105-118 (2014) ·Zbl 1291.65108号 [38] 阿梅斯托，P。；达夫，I。；L'Excellent，J.-Y。；Koster，J.，《使用分布式动态调度的完全异步多前沿解算器》，SIAM J.Matrix Ana。申请。，23、1、15-41（2001），网址http://mumps.enseeiht.fr ·Zbl 0992.65018号 [39] Stathopoulos，A。；Wu，K.，具有恒定同步要求的块正交化程序，SIAM J.Sci。计算。，23, 6, 2165-2182 (2002) ·Zbl 1018.65050号 [40] Gutknecht，M.H。；Schmelzer，T.，更新块三对角矩阵和块hessenberg矩阵的（Q R）分解，应用。数字。数学。，58, 6, 871-883 (2008) ·Zbl 1143.65032号 [41] NVIDIA，cuSPARSE网页（2020），URLhttps://docs.nvidia.com/cuda/cusparse [42] 布卡拉姆，W。；Turkiyyah，G。；Keyes，D.，《GPU上的层次矩阵运算：矩阵向量乘法和压缩》，ACM Trans。数学。软件，45（2019）·Zbl 1471.65027号 [43] Ambartsumyan，I。；布卡拉姆，W。；Bui-Thanh，T。；O.加塔斯。；凯斯，D。；斯塔德勒，G。；Turkiyyah，G。；Zampini，S.，在偏微分方程控制的反问题中产生的hessian的层次矩阵近似，SIAM J.Sci。计算。，42、5、A3397-A3426（2020）·Zbl 1453.65289号 [44] Davis，T.A.，算法832:UMFPACK——一种非对称模式的多面方法，ACM Trans。数学。软件，30，2，196-199（2004）·Zbl 1072.65037号 [45] 英特尔，MKL网页（2020），URLhttps://software.intel.com/content/www/us/en/develop/tools/math-kernel-library.html [46] Li，X.S.，《超级逻辑单元概述：算法、实现和用户界面》，ACM Trans。数学。软件，31，3，302-325（2005）·Zbl 1136.65312号 [47] Li，X.S。；Demmel，J.，Superlu_dist:非对称线性系统的可扩展分布式内存稀疏直接求解器，ACM Trans。数学。软件，29，2，110-140（2003）·Zbl 1068.90591号 [48] Poulson，J。；标记，B。；范德盖因，R.A。；哈蒙德，J.R。；Romero，N.A.，Elemental:分布式内存密集矩阵计算的新框架，ACM Trans。数学。软件，39，2（2013）·Zbl 1295.65137号 [49] 海涅克，A。；亨利·G。；Hutchinson，M。；Pabst，H.，LIBXSMM：通过运行时代码生成加速小矩阵乘法，（2016年高性能计算、网络、存储和分析国际会议论文集，SC16（2016），IEEE），URLhttps://github.com/hfp/libxsmm网站 [50] Hecht，F.，《FreeFem++的新发展》，J.Numer。数学。，20、3-4、251-266（2012），网址http://freefem.org ·Zbl 1266.68090号 [51] 穆林，J。；Jolivet，P。；Marquet，O.，大型水动力稳定性分析的增广拉格朗日预条件，计算。方法应用。机械。工程，351，718-743（2019），URLhttps://github.com/prj-/moulin2019al ·Zbl 1441.76069号 [52] Stewart，G.W.，大型特征问题的krylov-schur算法，SIAM J.矩阵分析。申请。，23, 3, 601-614 (2002) ·Zbl 1003.65045号 [53] Benzi，M。；Olshanskii，医学硕士。；Wang，Z.，不可压缩Navier-Stokes方程的修正增广拉格朗日预条件，国际。J.数字。《液体方法》，66，4，486-508（2011）·Zbl 1421.76152号 [54] 樱井，T。；Sugiura，H.，《使用数值积分求解广义特征值问题的投影方法》，J.Compute。申请。数学。，159, 1, 119-128 (2003) ·兹比尔1037.65040 [55] 亚当斯，M。；Bayraktar，H.H.公司。；Keaveny，T.M。；Papadopoulos，P.，《超过5亿自由度的固体力学中的超可伸缩隐式有限元分析》（2004年ACM/IEEE超级计算会议论文集，SC04（2004），IEEE计算机学会），34:1-34:15 [56] B.F.史密斯。；比约尔斯塔德，P。；Gropp，W.D.，《区域分解：椭圆偏微分方程的并行多层方法》（2004），剑桥大学出版社 [57] Dolean，V。；Jolivet，P。；Nataf，F.，《区域分解方法简介：算法、理论和并行实现》（2015），SIAM·Zbl 1364.65277号 [58] 托塞利，A。；Widlund，O.B.，（领域分解方法：算法和理论。领域分解方法，算法和理论，计算数学系列，第34卷（2005），Springer）·Zbl 1069.65138号 [59] Pechstein，C.，多尺度问题的有限元和边界元撕裂和互连解算器，计算科学与工程讲义（2012），Springer [60] Quarteroni，A。；Valli，A.，偏微分方程的区域分解方法，第10卷（1999年），Clarendon出版社·兹伯利0931.65118 [61] Mathew，T.，偏微分方程数值解的区域分解方法，第61卷（2008），Springer Science&Business Media·Zbl 1147.65101号 [62] 蔡，X.C。；Sarkis，M.，一般稀疏线性系统的限制加性Schwarz预条件，SIAM J.Sci。计算。，21, 2, 792-797 (1999) ·Zbl 0944.65031号 [63] Mandel，J.，平衡域分解，Commun。数字。方法。工程师，9，3，233-241（1993）·Zbl 0796.65126号 [64] Zampini，S.，PCBDDC：PETSc中的一类稳健双精度方法，SIAM J.Sci。计算。，38，5，S282-S306（2016），网址https://www.mcs.anl.gov/petsc/petsc-master/docs/manualpages/PC/PCBDDC.html ·Zbl 1352.65632号 [65] 佩奇斯坦，C。；Dohrmann，C.R.，自适应BDDC的统一框架，Electron。事务处理。数字。分析。，46, 273-336 (2017) ·Zbl 1368.65043号 [66] 扎皮尼，S。；瓦西列夫斯基，P。；Dobrev，V。；Kolev，T.，通过约束算法平衡任意阶卷曲一致空间的区域分解，（区域分解方法国际会议（2017），Springer），103-116·Zbl 1442.65435号 [67] 哦，D.-S。；Widlund，O.B。；扎皮尼，S。；Dohrmann，C.，《拉维亚特·托马斯向量场的BDDC算法与原始约束的豪华缩放和自适应选择》，数学。公司。，87, 310, 659-692 (2018) ·兹比尔1380.65065 [68] 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号 [69] Gee，M.W。；西弗特，C.M。；胡建杰。；杜米纳罗，R.S。；Sala，M.G.，ML 5.0平滑聚合用户指南技术。众议员SAND2006-2649（2006），网址https://trilinos.github.io/ml.html [70] 海因莱因，A。；Klawonn，A。；Rheinbach，O.，基于Trilinos的能量最小化粗糙空间的两级重叠Schwarz方法的并行实现，SIAM J.Sci。计算。，38、6、C713-C747（2016）·Zbl 1355.65168号 [71] 谢斯特克，J。；曼德尔，J。；Sousedík，B。；Burda，P.，多层BDDC的并行实现，（数值数学与高级应用2011（2013），Springer），681-689，URLhttp://users.math.cas.cz/sistek/software/bddcml.html ·Zbl 1311.65178号 [72] 巴迪亚，S。；马汀，A.F。；Principe，J.，通过约束平衡域分解的高度可扩展并行实现，SIAM J.Sci。计算。，36、2、C190-C218（2014）·Zbl 1296.65177号 [73] Gander，M.J.，优化的Schwarz方法，SIAM J.Numer。分析。，44, 2, 699-731 (2006) ·Zbl 1117.65165号 [74] 甘德，M.J。；Van Criekingen，S.，使用PETSc对优化限制添加剂Schwarz进行新的粗校正，（区域分解方法国际会议（2019年），Springer） [75] 斯皮兰，N。；Dolean，V。；Hauret，P。；F.纳塔夫。；佩奇斯坦，C。；Scheichl，R.，PDE系统的鲁棒两级区域分解预条件，Compt。R.数学。，349, 23, 1255-1259 (2011) ·Zbl 1252.65201号 [76] 斯皮兰，N。；Dolean，V。；Hauret，P。；F.纳塔夫。；佩奇斯坦，C。；Scheichl，R.，通过重叠中的广义特征问题抽象PDE系统的鲁棒粗糙空间，Numer。数学。，126, 4 (2013), 741-770 ·Zbl 1291.65109号 [77] Marchand，P。；克莱斯，X。；Jolivet，P。；纳塔夫，F。；Tournier，P.-H.，具有GenEO的超奇异算子的\（H\）型边界元逼近的两级预处理，Numer。数学。，146, 597-628 (2020) ·Zbl 1451.65215号 [78] Brown，J。；Knepley，M.G。；May，D.A。；Curfman McInnes，L。；Smith，B.F.，《用于多物理的可组合线性解算器》，（2012年第11届并行和分布式计算国际研讨会（2012），IEEE），55-62 [79] 巴特勒，R。；Dodwell，T。；Reinarz，A。；Sandhu，A。；Scheichl，R。；Seelinger，L.，高绩效沙丘用于解决大规模强各向异性椭圆问题的模块，并应用于航空复合材料，计算。物理学。社区。，249 (2020) ·Zbl 07678488号 [80] D.A.May、P.Sanan、K.Rupp、M.G.Knepley、B.F.Smith，《PETSc内的极端尺度多重网格组件》，载于《高级科学计算平台会议论文集》，2016年。 [81] Knepley，M.G。；Karpeev，D.A.，《筛网I的PDE网格算法：网格分布》，科学。程序。，17, 3, 215-230 (2009) [82] Tang，J。；纳本，R。；Vuik，C。；Erlangga，Y.，从通缩、区域分解和多重网格方法导出的两级预条件的比较，J.Sci。计算。，39, 3, 340-370 (2009) ·Zbl 1203.65073号 [83] Bastian，P。；Wittum，G。；Hackbusch，W.，加法和乘法多重网格——比较，计算，60，4，345-364（1998）·兹比尔0908.65107 [84] 哈弗萨斯，R。；Jolivet，P。；Nataf，F.，狮子算法的加性Schwarz方法类型理论和对称优化的限制加性Schvarz方法，J.Sci。计算。，39、4、A1345-A1365（2017）·Zbl 1371.65129号 [85] 柯南，L。；Dolean，V。；克劳斯，R。；Nataf，F.，基于Dirichlet-to-Neumann算子的异质Helmholtz问题的粗糙空间，J.Compute。申请。数学。，271, 83-99 (2014) ·Zbl 1321.65172号 [86] Al Daas，H。；格里戈里，L。；乔利维特，P。；Tournier，P.-H.，基于鲁棒粗糙空间层次结构的多级Schwarz预处理器，J.Sci。计算。（2019），提交出版。统一资源定位地址https://github.com/prj-/aldaas2019multi [87] 斯皮兰，N。；Rixen，D.J.，《稳健有限元撕裂和互连以及平衡区域分解算法的自动谱粗糙空间》，国际。J.数字。方法工程，95，11，953-990（2013）·Zbl 1352.65553号 [88] Karypis，G。；Kumar，V.，《划分不规则图的快速高质量多级方案》，SIAM J.Sci。计算。，20，1，359-392（1998），网址http://glaros.dtc.umn.edu/gkhome/metis/parmetis/overview ·Zbl 0915.68129号 [89] 安德森，R。；安德烈·J。；巴克，A。；Bramwell，J。；卡米尔，J.-S。；塞维尼，J。；Dobrev，V。；Y.杜杜伊特。；费希尔，A。；科列夫，T。；Pazner，W。；斯托维尔，M。；托莫夫，V。；阿克曼，I。；Dahm，J。；麦地那，D。；Zampini，S.，MFEM：模块化有限元方法库，计算。数学。申请。，81，42-74（2021），网址网址：http://mfem.org ·Zbl 1524.65001号 [90] 周，Y。；Saad，Y.，大型对称特征值问题的块Krylov-Schur方法，Numer。算法，47，4，341-359（2008）·Zbl 1153.65330号 此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。