MPI- and CUDA- implementations of modal finite difference method for P-SV wave propagation modeling

Among different discretization approaches, Finite Difference Method (FDM) is widely used for acoustic and elastic full-wave form modeling. An inevitable deficit of the technique, however, is its sever requirement to computational resources. A promising solution isparallelization, where the problem is broken into several segments, and the calculations are distributed over different processors. For thepresent FD routines, however, such parallelization techniqueinevitably needs domain-decomposition and inter-core data exchange, due to the coupling of the governing equations. In this study, a newFD-based procedure for seismic wave modeling, named as ‘ModalFinite Difference Method (MFDM) is introduced, which deals with the simulation in the decoupled modal space; thus, neither domaindecomposition nor inter-core data exchange is anymore required, which greatly simplifies parallelization for both MPI- and CUDA implementations over CPUs and GPUs. With MFDM, it is alsopossible to simply cut off less-significant modes and run the routine for just the important ones, which will effectively reduce computation and storage costs. The efficiency of the proposed MFDM is shown by some numerical examples.

Keywords:

Finite difference method , Graphics Processing Unit (GPU) , , Message Passing Interface (MPI) , Modal , Wave propagation