This paper describes two methods to improve the performance of a FDTD solver for the first order formulation of the 3D acoustic wave equation. Based on spatial and temporal cache blocking techniques, these methods enable to maximize bandwidth of the memory subsystem, while reducing data traffic in-between the memory hierarchy. On the one hand, the spatial blocking (SB) approach increases data reuse among cores within each iteration of the time integration. On the other hand, the multicore wavefront diamond temporal blocking (MWD-TB) technique further leverages the SB performance by intrinsically reusing freshly cached data solutions across iterations of the time integration. While SB achieves sixfold performance speedup against the naive implementation (without cache blocking), MWD-TB outperforms SB by up to 50\% on a two-socket 16-core Intel Haswell system.


Article metrics loading...

Loading full text...

Full text loading...


  1. Alford, R.M., Kelly, K.R. and Boore, D.M.
    [1974] Accuracy of finite-difference modeling of the acoustic wave equation. GEOPHYSICS, 39(6), 834–842.
    [Google Scholar]
  2. Chang, W. and McMechan, G.
    [1990] 3D acoustic prestack reverse-time migration. Geophysical Prospecting, 38(7), 737–755.
    [Google Scholar]
  3. Etienne, V., Tonellot, T., Malas, T., Ltaief, H., Kortas, S., Thierry, P. and Keyes, D.
    [2017] High-Performance Seismic Modeling with Finite-Difference Using Spatial and Temporal Cache Blocking. 3rd EAGE Workshop High Performance Computing for Upstream.
    [Google Scholar]
  4. Komatitsch, D. and Martin, R.
    [2007] An Unsplit Convolutional Perfectly Matched Layer Improved at Grazing Incidence for the Seismic Wave Equation. GEOPHYSICS, 72(5), SM155–SM167.
    [Google Scholar]
  5. Malas, T., Hager, G., Ltaief, H. and Keyes, D.E.
    [2017] Multidimensional Intratile Parallelization for Memory-Starved Stencil Computations. ACM Trans. Parallel Comput., 4(3), 12:1–12:32.
    [Google Scholar]
  6. Malas, T., Hager, G., Ltaief, H., Stengel, H., Wellein, G. and Keyes, D.
    [2015] Multicore Optimized Wavefront Diamond Blocking for Optimizing Stencil Updates. SIAM Journal on Scientific Computing, 37(4), 439–464.
    [Google Scholar]
  7. Virieux, J. and Operto, S.
    [2009] An overview of full-waveform inversion in exploration geophysics. Geophysics, 76(6), WCC1–WCC26.
    [Google Scholar]
  8. Zhang, W. and Luo, J.
    [2013] Full-waveform Velocity Inversion Based on the Acoustic Wave Equation. American Journal of Computational Mathematics, 03, 13–20.
    [Google Scholar]

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error