During the forward and backward modeling in Reverse Time Migration (RTM), stencil computations constitute one of the main computationally intensive components. Their classic implementation based on Spatial Blocking (SB) is subject to performance limitation on modern multicore architectures due to several reasons, including non-uniform memory access, memory bandwidth starvation, load imbalance, and limited data locality. The Multicore Wavefront Diamond-tiling Temporal Blocking technique (MWD-TB) introduced in (Malas, PhD thesis 2015, Malas et al., SIAM SciCo 2015, Malas et al., ACM Trans 2017) aims at reducing the memory bandwidth requirement of stencil computations by increasing cache reuse within successive time steps. The authors in (Akbudak et al. IJHPCA 2020) integrate the MWD-TB technique into the modeling phase and the authors in (Qu et al., KAUST Tech Report 2020) eventually embed it into the full RTM using in-memory I/O operations snapshotting for the imaging condition and illustrate with the Salt3D dataset. In this paper, we further enable Out-Of-Core (OOC) I/O snapshotting operations on the Lustre parallel file system using the buffering strategy from MLBS (Alturkestani et al., EuroPar 2020). We present preliminary results using the Marmoussi 3D dataset.


Article metrics loading...

Loading full text...

Full text loading...


  1. Versteeg, Roelof
    [1994] The Marmousi experience: Velocity model determination on a synthetic complex data set. The Leading Edge13 (9): 927–936.
    [Google Scholar]
  2. Malas, T.
    [2015] Tiling and Asynchronous Communication Optimizations for Stencil Computations, PhD thesis, King Abdullah University of Science and Technology.
    [Google Scholar]
  3. 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]
  4. 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]
  5. Alturkestani, T., Tonellot, T., Ltaief, H., Abdelkhalak, R., Etienne, V., and Keyes, D.
    [2019] MLBS: Transparent Data Caching in Hierarchical Storage for Out-of-Core HPC Applications, IEEE 26th International Conference on High Performance Computing, Data, and Analytics (HiPC), Hyderabad, India, pp. 312–322, doi: 10.1109/HiPC.2019.00046.
    https://doi.org/10.1109/HiPC.2019.00046 [Google Scholar]
  6. Alturkestani, T., Ltaief, H. and Keyes, D.
    [2020] Maximizing I/O bandwidth for reverse time migration on heterogeneous large-scale systems. Parallel Processing. Euro-Par 2020. Lecture Notes in Computer Science, 12247.
    [Google Scholar]
  7. AkbudakK, LtaiefH, EtienneV, AbdelkhalakR, TonellotT, KeyesD.
    [2020] Asynchronous computations for solving the acoustic wave propagation equation. The International Journal of High Performance Computing Applications. 34(4):377–393. doi:10.1177/1094342020923027.
    https://doi.org/10.1177/1094342020923027 [Google Scholar]
  8. Qu, L., Abdelkhalak, R., Ltaief, H., Said, I., and Keyes, D.
    [2020] High Performance Asynchronous Reverse Time Migration for Oil and Gas Exploration, KAUST Technical Report http://hdl.handle.net/10754/665739.
    [Google Scholar]
  9. Tian, X. B., Kang, I. B., Kim, G. Y. and Zhang, H. S.
    [2008] An Improvement in the Absorbing Boundary Technique for Numerical Simulation of Elastic Wave Propagation. Journal of Geophysics and Engineering5: 203–209.
    [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