Minimod: A Finite Difference solver for Seismic Modeling

This article introduces a benchmark application for seismic modeling using finite difference method, which is namedMiniMod, a mini application for seismic modeling. The purpose is to provide a benchmark suite that is, on one hand easy to build and adapt to the state of the art in programming models and changing high performance hardware landscape. On the other hand, the intention is to have a proxy application to actual production geophysical exploration workloads for Oil & Gas exploration, and other geosciences applications based on the wave equation. From top to bottom, we describe the design concepts, algorithms, code structure of the application, and present the benchmark results on different current computer architectures.

[1]  M. A. Dablain,et al.  The application of high-order differencing to the scalar wave equation , 1986 .

[2]  Jean Virieux,et al.  An overview of full-waveform inversion in exploration geophysics , 2009 .

[3]  J. Virieux.  P-SV wave propagation in heterogeneous media: Velocity‐stress finite‐difference method , 1986 .

[4]  Samuel Williams,et al.  Roofline: an insightful visual performance model for multicore architectures , 2009, CACM.

[5]  Mauricio Araya-Polo,et al.  Towards a Multi-Level Cache Performance Model for 3D Stencil Computation , 2011, ICCS.

[6]  Felix J. Herrmann,et al.  Interactive comment on “ Devito ( v 3 . 1 . 0 ) : an embedded domain-specific language for finite differences and geophysical exploration , 2018 .

[7]  D. Komatitsch,et al.  An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation , 2007 .

[8]  David B. Whalley,et al.  Tools for application-oriented performance tuning , 2001, ICS '01.

[9]  K. Yee.  Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[10]  Jean Virieux,et al.  SH-wave propagation in heterogeneous media; velocity-stress finite-difference method , 1984 .

[11]  Robert W. Graves,et al.  Simulating seismic wave propagation in 3D elastic media using staggered-grid finite differences , 1996, Bulletin of the Seismological Society of America.

[12]  Wei Liu,et al.  On the instability in second-order systems for acoustic VTI and TTI media , 2012 .

[13]  Henri Calandra,et al.  Experiences with OpenMP, PGI, HMPP and OpenACC Directives on ISO/TTI Kernels , 2012, 2012 SC Companion: High Performance Computing, Networking Storage and Analysis.

[14]  Mauricio Hanzich,et al.  Assessing Accelerator-Based HPC Reverse Time Migration , 2011, IEEE Transactions on Parallel and Distributed Systems.

[15]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[16]  Henri Calandra,et al.  Performance portability in reverse time migration and seismic modelling via OpenACC , 2017, Int. J. High Perform. Comput. Appl..

[17]  Pradeep Dubey,et al.  3.5-D Blocking Optimization for Stencil Computations on Modern CPUs and GPUs , 2010, 2010 ACM/IEEE International Conference for High Performance Computing, Networking, Storage and Analysis.

[18]  Chau-Wen Tseng,et al.  Tiling Optimizations for 3D Scientific Computations , 2000, ACM/IEEE SC 2000 Conference (SC'00).

[19]  V. I. Lebedev,et al.  Difference analogues of orthogonal decompositions, basic differential operators and some boundary problems of mathematical physics. I , 1964 .

[20]  Mauricio Araya-Polo,et al.  Algorithm 942 , 2014 .