GPU parallelization of a three dimensional marine CSEM code

Authors:
M. Sommer;S. HöLz;M. Moorkamp;A. Swidinsky;B. Heincke;C. Scholl;M. Jegen
Affiliations:
-;-;-;-;-;-;-
Venue:
Computers & Geosciences
Year:
2013

Citing 6
Cited 0

An iterative solution method for solving f(A)x=b, using Krylov subspace information obtained for the symmetric positive definite matrix A

Journal of Computational and Applied Mathematics
Matrix computations (3rd ed.)

Matrix computations (3rd ed.)
Unitary Triangularization of a Nonsymmetric Matrix

Journal of the ACM (JACM)
Massively parallel forward modeling of scalar and tensor gravimetry data

Computers & Geosciences
High-order finite-element seismic wave propagation modeling with MPI on a large GPU cluster

Journal of Computational Physics
Iterative Krylov solution methods for geophysical electromagnetic simulations on throughput-oriented processing units

International Journal of High Performance Computing Applications

Quantified Score

Hi-index	0.00

Visualization

Abstract

One of the main problem of 3D time domain controlled source electromagnetic (CSEM) inversion is the high runtimes of forward modeling codes. We reduced the runtime of the 3D time domain finite difference CSEM code TEMDDD by the GPU-parallelization of expensive algorithms. The code solves the electromagnetic diffusion equation by discretization of spatial operators and subsequent calculation of eigenpairs. These eigenpairs are found by approximation of the eigenspace in a Krylov subspace using the spectral Lanczos decomposition Method. This algorithm was in its original form not parallelizable due to implementation of the upper boundary condition at the air-water interface. We show for the marine case that replacing the original boundary condition at the air-water surface by a discretized air layer allows GPU parallelization of every time consuming algorithm of the code in the marine case. Speedups between 20 and 60 have been achieved compared to the original code for a larger 3-D model. In this model the bathymetry from a survey area offshore Egypt is used as an example demonstrating that the parallelized version of the code is applicable to real survey scenarios.