The current trend in seismic imaging aims at using an improved physical model, considering that the Earth is not rigid but an elastic body. This new model takes simulations closer to the real physics of the problem, at the cost of raising the needed computational resources. On the hardware front, recently developed high-performing devices, called accelerators or co-processors, have shown that can outperform their general purpose counterparts by orders of magnitude in terms of performance per watt. These new alternatives may then provide the necessary resources for making possible to represent complex wave physics in a reasonable time. There might be, however, a penalty associated to the usage of such devices, as some portion of the simulation code might need some re-writing or new optimization strategies explored and applied. In this work we will show some optimization strategies evaluated and applied to an elastic propagator based on a Fully Staggered Grid, running on the Intel® Xeon Phi(TM) coprocessor. It is important to remark, that the propagator is able to reproduce elastic wave propagation, even for an arbitrary anisotropy.


Article metrics loading...

Loading full text...

Full text loading...


  1. Aki, K., Richards, P.G.
    , 2002. Quantitative Seismology. University Science Books, Sausalito, California.
    [Google Scholar]
  2. Araya-Polo, M., Cabezas, J., Hanzich, M., Pericas, M., Rubio, F., Gelado, I., Shafiq, M., Morancho, E., Navarro, N., Ayguade, E., CelaJosé, M., Valero, M.
    , 2011. Assessing accelerator-based HPC reverse time migration. IEEE Trans. Parallel Distrib. Syst.22 (1), 147–162.
    [Google Scholar]
  3. Davydycheva, S., Druskin, V., Habashy, T.
    , 2003. An efficient finite-difference scheme for electromagnetic logging in 3D anisotropic inhomogeneous media. Geophysics68(5), 1525–1536.
    [Google Scholar]
  4. Klemm, M., Duran, A., Tian, X., Saito, H., Caballero, D., Martorell, X.
    , 2012. Extending OpenMP* with Vector Constructs for Modern Multicore SIMD Architectures. IWOMP‘12 Mercurium C/C++ Source-to-source Compiler. Accessed: 05/2015. http://pm.bsc.es/mcxx.
    [Google Scholar]
  5. Caballero, D., RoyuelaS., FerrerR., DuranA., Martorell, X.
    , 2015. Optimizing Overlapped Memory Accesses in User-directed Vectorization. International Conference on Supercomputing.
    [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