Efficiency of High Order Spectral Element Methods on Petascale Architectures

High order methods for the solution of PDEs expose a trade-off between computational cost and accuracy on a per degree of freedom basis. In many cases, the cost increases due to higher arithmetic intensity while affecting data movement minimally. As architectures tend towards wider vector instructions and expect higher arithmetic intensities, the best order for a particular simulation may change.

[1]  O A Hurricane,et al.  Panel 3 Report: Implosion Hydrodynamics , 2012 .

[2]  A. Patera.  A spectral element method for fluid dynamics: Laminar flow in a channel expansion , 1984 .

[3]  Michel Schanen,et al.  On the Strong Scaling of the Spectral Element Solver Nek5000 on Petascale Systems , 2016, EASC.

[4]  H.M. Tufo,et al.  Terascale Spectral Element Algorithms and Implementations , 1999, ACM/IEEE SC 1999 Conference (SC'99).

[5]  John M. Levesque,et al.  An MPI/OpenACC implementation of a high-order electromagnetics solver with GPUDirect communication , 2016, Int. J. High Perform. Comput. Appl..

[6]  Chun Chen,et al.  Speeding up Nek5000 with autotuning and specialization , 2010, ICS '10.

[7]  Rémi Abgrall,et al.  High‐order CFD methods: current status and perspective , 2013 .

[8]  Paul Linden,et al.  On the structure of salt fingers , 1973 .

[9]  Robert A. van de Geijn,et al.  Anatomy of high-performance matrix multiplication , 2008, TOMS.

[10]  G. Henry,et al.  LIBXSMM: A High Performance Library for Small Matrix Multiplications , 2015 .

[11]  Leonhard Rannabauer,et al.  High-Order ADER-DG Minimizes Energy- and Time-to-Solution of SeisSol , 2015, ISC.

[12]  M. Zingale,et al.  Direct Numerical Simulations of Type Ia Supernovae Flames. II. The Rayleigh-Taylor Instability , 2004, astro-ph/0401247.

[13]  Michael Zingale,et al.  A comparative study of the turbulent Rayleigh-Taylor instability using high-resolution three-dimensional numerical simulations: The Alpha-Group collaboration , 2004 .

[14]  Erwin Laure,et al.  OpenACC acceleration of the Nek5000 spectral element code , 2015, Int. J. High Perform. Comput. Appl..

[15]  Maxwell Hutchinson.  Direct numerical simulation of single mode three-dimensional Rayleigh-Taylor experiments , 2015 .

[16]  Paul F. Fischer,et al.  Hybrid Multigrid/Schwarz Algorithms for the Spectral Element Method , 2005, J. Sci. Comput..

[17]  Erwin Laure,et al.  Evaluation of Parallel Communication Models in Nekbone, a Nek5000 Mini-Application , 2015, 2015 IEEE International Conference on Cluster Computing.