×
迈克尔·兰格；劳伦斯·米切尔；马修·克内普利。；Gerard J.戈尔曼。

使用PETSc DMPlex在firedrake中进行高效网格管理。 （英语） Zbl 1352.65613号

SIAM J.科学。计算。 38，编号5，S143-S155（2016）.
摘要：使用可组合抽象可以将新的和已建立的算法应用于广泛的问题，同时自动继承众所周知的性能优化的优点。这项工作强调了PETSc DMPlex域拓扑抽象与firedrake自动有限元系统的组合，以创建一个结合了表达性、灵活性和高性能的偏微分方程求解环境。我们描述了firedrake如何利用DMPlex提供有限元装配所需的间接映射，同时支持各种网格输入格式和运行时域分解。特别是，我们描述了DMPlex及其附带的数据结构如何允许通用地创建用户定义的离散化，同时利用数据布局优化来提高缓存一致性并确保在汇编计算期间进行重叠通信。

引用于11文件

MSC公司：

65牛顿50 涉及偏微分方程的边值问题的网格生成、细化和自适应方法
65N30型 含偏微分方程边值问题的有限元、Rayleigh-Ritz和Galerkin方法
65号55 多重网格方法；含偏微分方程边值问题的区域分解
第65年 并行数值计算
65日元 数值算法的封装方法

关键词：

网格；拓扑；分区；重新编号；火蜥蜴；PETSc公司；算法；性能优化；有限元；区域分解

软件：

Firedrake公司；OP2（操作2）；UFL公司；PyOP2型；PUMI公司；MOAB公司；DMPlex公司；符号控制；菲亚特；Gms小时；DOLFIN公司；PETSc公司；炒作；FEniCS公司
PDF格式BibTeX公司 XML格式引用
\textit{M.Lange}等人，SIAM J.Sci。计算。38，第5号，S143--S155（2016；Zbl 1352.65613）
全文： 内政部 arXiv公司

参考文献：

[1] M.S.Aln \aes、A.Logg、K.B.Ølgaard、M.E.Rognes和G.N.Wells，《统一形式语言：偏微分方程弱公式的领域特定语言》，ACM Trans。数学。软件，40（2014），9。http://dl.acm.org/citation.cfm？id=2566630doi:10.1145/2566630·兹比尔1308.65175
[2] S.Balay、S.Abhyankar、M.F.Adams、J.Brown、P.Brune、K.Buschelman、L.Dalcin、V.Eijkhout、W.Gropp、D.Kaushik、M.Knepley、L.C.McInnes、K.Rupp、B.Smith、S.Zampini和H.Zhang，{it PETSc用户手册}，技术报告ANL-95/11-3.6版，阿尔贡国家实验室数学和计算机科学部，伊利诺伊州莱蒙特，2015年。可在线访问http://www.mcs.anl.gov/petsc/petsc-3.6/docs/manual.pdf。
[3] J.Brown，《星森林作为平行传播模型》，技术报告，2011年。http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20111224/7bfa1e70/attachent.pdf。
[4] J.Brown、M.G.Knepley、D.A.May、L.C.McInnes和B.Smith，{多物理的可组合线性解算器}，《第十一届并行与分布式计算国际研讨会论文集》，IEEE计算机学会，华盛顿特区，2012年，第55-62页。http://ieeexplore.ieee.org/xpl/articleDetails.jsp？arnumber=6341494doi:10.1109/ISPDC.2012.16。
[5] J.Brown、M.G.Knepley和B.F.Smith，{仿真软件的运行时扩展性和库化}，IEEE Compute。科学。工程，17（2015），第38-45页。
[6] P.R.Brune、M.G.Knepley和L.R.Scott，{复杂和分级网格上二维和三维非结构化几何多重网格}，SIAM J.Sci。计算。，35（2013），第A173-A191页。http://epubs.siam.org/doi/abs/10.1137/110827077doi:10.1137/10827077·Zbl 1264.65203号
[7] E.Cuthill和J.McKee，{减少稀疏对称矩阵的带宽}，《第24届全国会议论文集》（ACM’69），ACM，纽约，1969年，第157-172页。
[8] R.D.Falgout、J.E.Jones和U.M.Yang，{\it《并行计算机上偏微分方程数值解中并行高性能预处理程序库hypre的设计与实现》，a.M.Bruaset和a.Tveito编辑，《计算讲义》。科学。Eng.51，Springer，柏林，2006年，第267-294页·Zbl 1097.65059号
[9] A.George和J.W.Liu，{大型稀疏正定系统的计算机解}，Prentice-Hall Ser。计算。数学。，普伦蒂斯·霍尔，恩格尔伍德悬崖，新泽西州，1981年·Zbl 0516.65010号
[10] C.Geuzaine和J.-F.Remacle，{\it Gmsh:一个内置预处理和后处理设施的三维有限元网格生成器}，国际。J.数字。方法工程师，79（2009），第1309-1331页。http://onlinelibrary.wiley.com/doi/10.1002/nme.2579/摘要doi:10.1002/nme.2579·Zbl 1176.74181号
[11] F.Guönther、M.Mehl、M.Po¨gl和C.Zenger，{it基于空间填充曲线的分层数据结构上PDE的缓存算法}，SIAM J.Sci。计算。，28（2006），第1634-1650页。http://epubs.siam.org/doi/abs/10.1137/0406078doi:10.1137/040604078·Zbl 1162.68406号
[12] G.Haase，M.Liebmann，and G.Plank，{非结构化稀疏矩阵的Hilbert阶乘法格式}，Inter。J.并行紧急配送系统。，22（2007），第213-220页·Zbl 1122.65037号
[13] D.A.Ibanez、E.S.Seol、C.W.Smith和M.S.Shephard，提交了《并行非结构化网格基础设施》·Zbl 1369.65155号
[14] R.C.Kirby，{算法839:Fiat，计算有限元基函数的新范式}，ACM Trans。数学。《软件》，30（2004），第502-516页。http://dl.acm.org/citation.cfm？doid=1039813.1039820doi:10.1145/1039813.1039820·Zbl 1070.65571号
[15] M.G.Knepley和D.A.Karpeev，{筛网I的PDE网格算法：网格分布}，科学。程序。，17（2009），第215-230页。
[16] M.G.Knepley、M.Lange和G.J.Gorman，｛\it并行的非结构化重叠网格分布｝，提交。预打印版本在线提供http://arxiv.org/abs/1506.06194arXiv:1506.06194[cs.MS]。
[17] M.Lange、M.G.Knepley和G.J.Gorman，《使用PETSc DMPlex}实现灵活、可伸缩的网格和数据管理》，《第三届Exascale应用程序和软件国际会议论文集》（EASC 2015），英国爱丁堡大学，2015年，第71-76页。可在线访问http://www.easc2015.ed.ac.uk/sites/default/files/attachments/EASC15Proceedings.pdf。
[18] A.Logg，{\it计算网格的有效表示}，国际。J.计算。科学。工程，4（2009），第283-295页。http://www.inderscience.com/offer.php？id=29164doi:10.10504/IJCSE.2009.029164。
[19] A.Logg、G.N.Wells和J.Hake，{it DOLFIN:A C++/Python有限元库}，《用有限元方法自动求解微分方程》，A.Logg，K.-A.Mardal和G.Wells.编辑，《计算讲义》。科学。工程84，施普林格，柏林，2012年，第173-225页·Zbl 1247.65105号
[20] G.R.Mudalige、M.B.Giles、I.Reguly、C.Bertolli和P.H.J.Kelly，{it OP2:在多核和多核架构上解决基于网格的非结构化应用程序的主动库框架}，载于《创新并行计算学报》（InPar），IEEE计算机学会，华盛顿特区，2012年，第1-12页。
[21] C.Ollivier-Gooch，{\it A mesh-database-independent edge-and-face-swapping tool}，载于美国航空航天局第44届航空航天科学会议和展览论文集，内华达州里诺，2006年，美国航空航天协会2006-533。
[22] C.Ollivier-Goch，L.Diachin，M.S.Shephard，T.Tautges，J.Kraftcheck，V.Leung，X.Luo，和M.Miller，{一种用于网格查询和操作的可互操作的数据结构中性组件}，ACM Trans。数学。《软件》，37（2010），29。http://dl.acm.org/citation.cfm？id=1864430doi:10.1145/1824801.1864430。
[23] F.Rathgeber、D.A.Ham、L.Mitchell、M.Lange、F.Luporini、A.T.T.McRae、G.-T.Bercea、G.R.Markall和P.H.J.Kelly，提交了《Firedrake:通过构成抽象实现有限元方法自动化》。预打印版本在线提供http://arXiv.org/abs/1501.01809arxiv:1501.01809[cs.MS]·Zbl 1396.65144号
[24] F.Rathgeber、G.R.Markall、L.Mitchell、N.Loriant、D.A.Ham、C.Bertolli和P.H.J.Kelly，{it PyOP2:非结构化网格上性能可移植模拟的高级框架}，《高性能计算、网络、存储和分析学报》（2012 SC Companion），IEEE Computer Society，华盛顿特区，第1116-1123页。http://ieeexplore.ieee.org/xpl/articleDetails.jsp？arnumber=6495916doi:10.1109/SC.Companion.2012.134。
[25] T.J.Tautges、R.Meyers、K.Merkley、C.Stimpson和C.Ernst，{\it MOAB:A Mesh-Oriented Database}，技术报告SAND2004-1592，Sandia National Laboratories，Albuquerque，NM，2004年。
[26] S.-E.Yoon、P.Lindstrom、V.Pascucci和D.Manocha，{\it Cache-obliovious mesh layouts}，ACM Trans。图表。，24（2005），第886-893页。http://dl.acm.org/citation.cfm？id=1073278doi:10.1145/1073204.1073278。
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。