斯特凡·德斯科姆;马克斯·杜阿尔特;蒂埃里·杜蒙特;吉列，托马斯;紫罗兰百叶窗;马克·马索特 多核结构上时空多尺度反应扩散系统的基于任务的自适应多分辨率。 （英语） 兹伯利1416.65184 SMAI J.计算。数学。 3, 29-51 (2017). 摘要：最近引入了一种具有时空自适应和误差控制功能的新型求解器，用于求解刚性反应扩散系统的数值解。该策略基于算子分裂、有限体积自适应多分辨率和高阶时间积分器，每个算子具有特定的稳定性特性，在标准体系结构（如功能强大的工作站）上进行大规模多维计算时具有较高的计算效率。然而，基于指针树的原始实现的数据结构为提高效率提供了有限的机会，同时在并行编程和负载平衡方面提出了严重的挑战。本文提出了包括Radau5和ROCK4在内的整套数值方法的一种新的实现方法，它依赖于完全不同的数据结构，并使用特定的库TBB来实现共享内存、基于任务的并行工作密封。在多核和多核架构上，我们的实现通过一系列在二维和三维上难度不断增加的测试用例进行了性能评估，证明了其具有很高的可扩展性。 引用于2文件 MSC公司： 65升04 刚性方程的数值方法 2005年5月 并行数值计算 65平方米 偏微分方程初值和初边值问题的线法 65M50型 涉及偏微分方程初值和初边值问题数值解的网格生成、细化和自适应方法 35K57型 反应扩散方程 软件：RAMSES公司;TBB（待定）;英特尔TBB;AMROC公司;皮亚诺;p4测试;浣熊;雷达5;梅子;西尔克;帕拉梅什;罗德斯 PDF格式BibTeX公司 XML格式引用 \textit{S.Descombes}等人，SMAI J.计算。数学。3、29-51（2017；Zbl 1416.65184） 全文： 内政部 arXiv公司 参考文献： [1] A.Abdulle，“具有递推关系的四阶切比雪夫方法”，SIAM J.Sci。计算23（2002）第6号，第2041-2054页（电子版）·Zbl 1009.65048号 ·doi:10.1137/S1064827500379549 [2] S.B.Baden、N.P.Chrisochoides、D.B.Gannon和M.L.Norman（编辑），结构化自适应网格细化（SAMR）网格方法，数学中的IMA卷及其应用117，Springer-Verlag，纽约，2000年，明尼苏达州明尼阿波利斯市明尼苏打大学举行的IMA研讨会论文，1997年3月12日至13日·Zbl 0930.00061号 ·doi:10.1007/978-1-4612-1252-2 [3] D.S.Balsara和C.D.Norton，“使用基于并行语言的方法进行高度并行结构自适应网格细化”，《并行计算》27（2001）第1-2期，第37-70页·Zbl 0971.68017号 ·doi:10.1016/S0167-8191（00）00088-0 [4] J.Bell、M.Berger、J.Saltzman和M.Welcome，“双曲守恒律的三维自适应网格细化”，SIAM J.Sci。计算15（1994）第1号，第127-138页·Zbl 0793.65072号 ·doi:10.137/0915008 [5] R.D.Blumofe、C.F.Joerg、B.C.Kuszmaul、C.E.Leiserson、K.H.Randall和Y.Zhou，“Cilk：高效的多线程运行时系统”，J.Parallel District Com.37（1996）第1期，第55-69页·doi:10.1006/jpdc.1996.0107 [6] R.D.Blumofe和C.E.Leiserson，“通过工作窃取调度多线程计算”，J.ACM46（1999）第5期，第720-748页·Zbl 1065.68504号 ·doi:10.1145/324133.324234 [7] K.Brix，S.Melian，S.Müller&M.Bachmann，自适应多分辨率方法：数据结构、实现和并行化的实际问题，多分辨率和自适应网格细化方法暑期学校，ESAIM Proc。34，EDP科学。，Les Ulis，2011年，第151-183页·Zbl 1302.65219号 [8] K.Brix，S.S.Melian，S.Müller&G.Schieffer，使用空间填充曲线对基于多尺度的网格自适应进行并行化，收录于ESAIM:Proc。，EDP科学，2009年，第108-129页·兹比尔1181.65118 [9] C.Burstede、L.Wilcox和O.Ghattas，“p4est:八叉树森林上并行自适应网格细化的可缩放算法”，SIAM J.Sci。计算33（2011）第3号，第1103-1133页·Zbl 1230.65106号 ·数字对象标识代码：10.1137/100791634 [10] A.Cohen、S.M.Kaber、S.Muller和M.Postel，“守恒定律的完全自适应多分辨率有限体积格式”，数学。计算72（2003），第183-225页·Zbl 1010.65035号 ·doi:10.1090/S0025-5718-01-01391-6 [11] W.Crutchfield和M.L.Welcome，“自适应网格细化算法的面向对象实现”，科学杂志。程序。（1993），第145-156页 [12] W.Dahmen、T.Gotzen、S.S.Melian-Flamand和S.Müller，使用自适应多尺度技术进行冷却气体喷射的数值模拟，法国几何与实践数学研究所，2011年 [13] R.Deiterding，面向对象C++中块结构自适应网格优化的通用框架，技术报告，网址：http://amroc.sourceforge.net, 2003 [14] R.Deiterding，块结构自适应网格细化-理论、实现和应用，多分辨率和自适应网格细化方法暑期学校，ESAIM Proc。34，EDP科学。，Les Ulis，2011年，第97-150页·Zbl 1302.65220号 ·doi:10.1051/proc/201134002 [15] S.Descombes，M.Duarte，T.Dumont，F.Laurent，V.Louvet&M.Massot，“非线性反应扩散方程非共形域中算子分裂的分析。预混合火焰动力学的应用”，SIAM J.Numer。分析52（2014）第3期，第1311-1334页·Zbl 1320.65080号 ·doi:10.137/130926006 [16] S.Descombes和T.Dumont，“中风数值模拟：计算问题和方法”，Prog。生物物理学。《分子生物学》97（2008）第1期，第40-53页·doi:10.1016/j.pbiomolbio.2007.10.003 [17] S.Descombes，T.Dumont，V.Louvet&M.Massot，“关于具有高空间梯度的反应扩散方程分裂近似的局部和全局误差”，Int.J.Comput。数学84（2007）第6期，第749-765页·Zbl 1122.65061号 ·doi:10.1080/00207160701458716 [18] S.Descombes、T.Dumont和M.Massot，刚性非线性反应扩散系统的算子分裂：阶数减少和螺旋波、模式和波的应用（圣彼得堡，2002），AkademPrint，圣彼得堡 [19] S.Descombes和M.Massot，“具有熵结构的非线性反应扩散系统的算子分裂：奇异摄动和降阶”，数值。数学97（2004）第4期，第667-698页·Zbl 1060.65105号 ·doi:10.1007/s00211-003-0496-3 [20] J.Dreher&R.Grauer，“浣熊：双曲守恒律的平行网格自适应框架”，《平行计算》31（2005）第8-9期，第913-932页·doi:10.1016/j.parco.2005.04.011 [21] M.-A.Dronne，J.-P.Boissel&E.Grenier，“中风期间灰质中离子运动的数学模型”，J.Theor。《生物学》240（2006）第4期，第599-615页·Zbl 1447.92087号 ·doi:10.1016/j.jtbi.2005.10.023 [22] F.Drui，A.Fikl，P.Kestener，S.Kokh，A.Larat，V.Le Chenadec&M.Massot，“两相流AMR模拟p4est库的实验”，ESAIM:Proc。Surv.53（2016），第232-247页·Zbl 1408.76380号 ·doi:10.1051/proc/201653014 [23] M.Duarte，《用于多尺度反应前沿模拟的时空自适应数值方法》，博士论文，巴黎中央大学，2011年 [24] M.Duarte、Z.Bonaventura、M.Massot、A.Bourdon、S.Descombes和T.Dumont，“多尺度拖缆放电模拟中具有时空自适应性和误差控制的新数值策略”，J.Compute。Phys.231（2012）第3期，第1002-1019页·Zbl 1406.76063号 ·doi:10.1016/j.jcp.2011.07.002 [25] M.Duarte、S.Descombes、C.Tenaud、S.Candel和M.Massot，“燃烧锋面模拟的时空自适应数值方法”，Combust。Flame160（2013）第6期，第1083-1101页·doi:10.1016/j.com.bustframe.2013.01.013 [26] M.Duarte，M.Massot，S.Descombes，C.Tenaud，T.Dumont，V.Louvet&F.Laurent，使用时间算子分裂和空间自适应多分辨率的多尺度反应波新分辨率策略：在人类缺血性中风中的应用，多分辨率和自适应网格精化方法暑期学校，ESAIM Proc。34，EDP科学。，Les Ulis，2011年，第277-290页·Zbl 1302.65211号 ·doi:10.1051/proc/201134006 [27] M.Duarte，M.Massot，S.Descombes，C.Tenaud，T.Dumont，V.Louvet&F.Laurent，“使用时间算子分裂、空间自适应多分辨率和专用高阶隐式/显式时间积分器的多尺度反应波新分辨率策略”，SIAM J.Sci。计算结果34（2012）第1号，第A76-A104页·兹比尔1243.65107 ·doi:10.1137/100816869 [28] T.Dumont、M.Duarte、S.Descombes、M.-A.Dronne、M.Massot和V.Louvet，“真实3D几何中人体缺血性中风的模拟”，Commun。非线性科学。数字。Simul.18（2013）第6期，第1539-1557页·兹比尔1322.92019 ·doi:10.1016/j.cnsns.2012.10.002 [29] W.Eckhardt&T.Weinzierl，《动态自适应PDE求解器的多核架构阻塞策略》，D.Hutchison等人主编，《并行处理和应用数学》，施普林格-柏林-海德堡出版社，2010年，第567-575页 [30] M.Essadki、S.de Chaisemartin、M.Massot、F.Laurent和A.Larat，“模拟多分散蒸发喷雾的自适应网格细化和高阶几何矩方法”，石油天然气科学。Technol.71（2016），第1-25页·doi:10.2516/ogst/2016012 [31] C.J.Forster，具有相变的多相流的并行小波自适应直接数值模拟，博士论文，佐治亚理工学院，2016 [32] T.Gautier，J.V.F.Lima，N.Maillard&B.Raffin，《XKaapi：异构架构上数据流任务编程的运行时系统》，并行分布式处理（IPDPS），2013年IEEE第27届国际研讨会，2013年，第1299-1308页·doi:10.1109/IPDPS.2013.66 [33] P.Gray和S.K.Scott，《化学振荡和不稳定性：非线性化学动力学》，国际化学专著系列21，牛津大学出版社，1994年 [34] E.Hairer，“Fortran和Matlab代码”， [35] E.Hairer，S.P.Nørsett和G.Wanner，求解常微分方程。一、 计算数学中的Springer级数8，Springer-Verlag，柏林，1993，非刚性问题·兹比尔0789.65048 [36] E.Hairer和G.Wanner，《求解常微分方程》。二、 计算数学中的Springer级数14，Springer-Verlag，柏林，1996，Stiff和微分代数问题·Zbl 0859.65067号 [37] A.Harten，“双曲守恒律数值解的多分辨率算法”，《纯粹应用数学通讯》48（1995），第1305-1342页·Zbl 0860.65078号 ·doi:10.1002/cpa.3160481201 [38] B.Hejazialhosseini，D.Rossinelli，M.Bergdorf&P.Koumoutsakos，“用于模拟冲击波相互作用的多核结构上具有局部时间步进的小波自适应网格上的高阶有限体积方法”，J.Compute。Phys.229（2010）第22号，第8364-8383页·Zbl 1381.76218号 ·doi:10.1016/j.jcp.2010.07.021 [39] 英特尔公司，“英特尔线程构建块”， [40] W.Jahnke，W.E.Skaggs&A.T.Winfree，“Belousov-Zhabotinskii反应和双变量Oregonator模型中的化学涡旋动力学”，J.Phys。《化学》93（1989）第2期，第740-749页·doi:10.1021/j100339a047 [41] H.Ji，F.-S.Lien和E.Yee，“一种新的自适应网格细化数据结构及其在爆炸中的应用”，J.Compute。Phys.229（2010）第23号，第8981-8993页·Zbl 1207.80023号 ·doi:10.1016/j.jcp.2010.08.023 [42] S.Jones&A.Lichtl，GPU到火星，SpaceX火星火箭发动机的全尺寸模拟，GPU技术会议记录，GPU按需技术会议，2015年 [43] R.Keppens、Z.Meliani、A.J.van Marle、P.Delmont、A.Vlasis和B.van der Holst，“相对论流体力学和磁流体力学的并行网格自适应方法”，J.Compute。Phys.231（2012）第3期，第718-744页·Zbl 1426.76385号 ·doi:10.1016/j.jp.2011.01.020 [44] P.MacNeice、K.M.Olson、C.Mobarry、R.de Fainchtein和C.Packer，“PARAMESH:并行自适应网格优化社区工具包”，计算。物理学。Commun.126（2000）第3号，第330-354页·Zbl 0953.65088号 ·doi:10.1016/S0010-4655（99）00501-9 [45] 面向计算机的大地测量数据库G.M.Morton；和文件排序的新技术，技术报告，IBM，1966年 [46] S.Müller，《保护法的自适应多尺度方案》27，斯普林格-Verlag，海德堡，2003年·Zbl 1016.76004号 [47] A.Nejadmalayeri、A.Vezolainen、E.Brown-Dymkoski和O.V.Vasilyev，“PDE的并行自适应小波配置方法”，J.Compute。《物理学》298（2015），第237-253页·Zbl 1349.65652号 ·doi:10.1016/j.jcp.2015.05.028 [48] Netlib，“压缩行存储（CRS）”， [49] R.M.Noyes、R.Field和E.Koros，“化学系统中的振荡。I.显示时间振荡的系统中的详细机制”，J.Ame。化学。Soc.94（1972）第4号，第1394-1395页·doi:10.1021/ja00759a080 [50] L.Oliker&R.Biswas，“PLUM:自适应非结构化网格的并行负载平衡”，J.Parallel District Com.52（1998）第2期，第150-177页·Zbl 0920.68016号 ·doi:10.1006/jpdc.1998.1469 [51] OpenMP Architecture Review Board，“OpenMP Application Program Interface.Version 3.1”，2011年7月，在线阅读： [52] J.Reinders，Intel Threading Building Blocks，O'Reilly&Associates，Inc.，加利福尼亚州塞巴斯托波尔，美国 [53] C.A.Rendleman、V.E.Beckner、M.Lijewski、W.Crutchfield和J.B.Bell，“结构化分层自适应网格细化算法的并行化”，计算。视觉。科学3（2000）第3号，第147-157页·Zbl 0971.65089号 ·doi:10.1007/PL00013544 [54] D.Rossinelli、B.Hejazialhosseni、M.Bergdorf和P.Koumoutsakos，“复杂系统仿真多核架构上的小波自适应解算器”，Concurr。计算：实际。Exper.23（2011）第2期，第172-186页·Zbl 1381.76218号 ·doi:10.1002/cpe.1639 [55] O.Roussel、K.Schneider、A.Tsigulin和H.Bockhorn，“抛物线偏微分方程的保守全自适应多分辨率算法”，J.Compute。《物理学》188（2003）第2期，第493-523页·Zbl 1022.65093号 ·doi:10.1016/S0021-9991（03）00189-X [56] H.Sagan，《空间填充曲线》，Universitext，Springer-Verlag，纽约，1994年·Zbl 0806.01019号 ·doi:10.1007/978-1-4612-0871-6 [57] K.Schneider和O.V.Vasilyev，“计算流体动力学中的小波方法”，Annu。《流体力学修订版》42（2010）第1期，第473-503页·Zbl 1345.76085号 ·doi:10.1146/annurev-fluid-121108-145637 [58] M.Schreiber、T.Weinzierl和H.-J.Bungartz，《动态自适应网格模拟的聚类优化和并行化》，D.Hutchison et al.，ed.，Euro Par 2013 Parallel Processing，Springer Berlin Heidelberg，2013，第484-496页 [59] G.Strang，“关于差分格式的构造和比较”，SIAM J.Numer。分析5（1968年），第506-517页·Zbl 0184.38503号 ·doi:10.1137/0705041 [60] R.Teyssier，“具有自适应网格细化的宇宙流体动力学：一种称为RAMSES的新的高分辨率代码”，Astron。天体物理学385（2002）第1期，第337-364页·doi:10.1051/0004-6361:20011817 [61] A.I.Volpert、V.A.Volpert和V.A.Volpert，抛物线系统的行波解，美国数学学会，普罗维登斯，RI，1994年，从詹姆斯·海达的俄罗斯手稿翻译而来·Zbl 0805.35143号 [62] T.Weinzierl和M.Mehl，“类八进制自适应笛卡尔多尺度网格的Peano-A遍历和存储方案”，SIAM J.Sci。计算单位33（2011）第5号，第2732-2760页·Zbl 1245.65169号 ·doi:10.1137/100799071 [63] S.Williams、A.Waterman和D.Patterson，“Roofline:多核架构的有见地的视觉性能模型”，Commun。ACM52（2009）第4号，第65-76页·数字对象标识代码：10.1145/1498765.1498785 [64] A.M.Wissink、R.D.Hornung、S.R.Kohn、S.S.Smith和N.Elliott，《使用SAMRAI框架进行大规模并行结构AMR计算》，超级计算，ACM/IEEE 2001年会议，2001年，第22-22页·doi:10.1109/SC.2001.10029 [65] 是的。B.Zel′dovich，G.I.Barenblatt，V.B.Librovich&G.M.Makhviladze，燃烧和爆炸的数学理论，顾问局[全体会议]，纽约，1985年，Donald H.McNeill从俄语翻译 此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。