High-Order Finite-differences on multi-threaded architectures using OCCA

High-order finite-difference methods are commonly used in wave propagator for industrial subsurface imaging algorithms. Computational aspects of the reduced linear elastic vertical transversely isotropic propagator are considered. Thread parallel algorithms suitable for implementing this propagator on multi-core and many-core processing devices are introduced. Portability is addressed through the use of the OCCA runtime programming interface. Finally, performance results are shown for various architectures on a representative synthetic test case.

[1]  Kevin Skadron,et al.  Scalable parallel programming , 2008, 2008 IEEE Hot Chips 20 Symposium (HCS).

[2]  R. P. Fletcher,et al.  A New Pseudo-acoustic Wave Equation for VTI Media , 2008 .

[3]  Moshe Reshef,et al.  A nonreflecting boundary condition for discrete acoustic and elastic wave equations , 1985 .

[4]  R. M. Alford,et al.  ACCURACY OF FINITE‐DIFFERENCE MODELING OF THE ACOUSTIC WAVE EQUATION , 1974 .

[5]  Peter Wittek,et al.  GPU Technology Conference , 2013 .

[6]  John E. Stone,et al.  OpenCL: A Parallel Programming Standard for Heterogeneous Computing Systems , 2010, Computing in Science & Engineering.

[7]  B. Fornberg.  The pseudospectral method: Comparisons with finite differences for the elastic wave equation , 1987 .

[8]  Vadim Lisitsa,et al.  Numerical study of the interface errors of finite-difference simulations of seismic waves , 2014 .

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

[10]  Timothy C. Warburton,et al.  OCCA: A unified approach to multi-threading languages , 2014, ArXiv.

[11]  Olav Holberg,et al.  COMPUTATIONAL ASPECTS OF THE CHOICE OF OPERATOR AND SAMPLING INTERVAL FOR NUMERICAL DIFFERENTIATION IN LARGE-SCALE SIMULATION OF WAVE PHENOMENA* , 1987 .

[12]  Bengt Fornberg,et al.  Classroom Note: Calculation of Weights in Finite Difference Formulas , 1998, SIAM Rev..

[13]  Albert Farrés,et al.  Roofline-based Optimizations for Elastic Propagation on Xeon , 2015 .