
Full text loading...
Improving both the accuracy and computational performance of simulation tools is a major challenge for seismic imaging, and generally requires specialized algorithms and computational implementations to make full use of modern hardware architectures. We present a computational strategy based on a high-order discontinuous Galerkin time-domain method. Our implementation can be run on several architectures thanks to a unified multi-threading programming framework, and exhibits a good load balancing and minimum data movements. A first benchmark validates this implementation and confirms the interest of accelerators in computational geophysics.