Solution of Large-scale Seismic Modeling Problems

The grid-characteristic method is commonly used to accurately model elastic waves propagation within the application of mechanical stress tests, collision of solid bodies and others. In this note, we discuss the application of the grid-characteristic method to numerical simulation of seismic waves which is an important problem in geophysics. A computational code for this application has been developed in our team. We studied its low-level optimization with streaming SIMD instructions (SSE and AVX). Beyond, we parallelized our code with MPI and reached parallel efficiency of up to 70% when using 16’000 computational cores.