The solution of the Marmousi problem on GPU's

During the PhD research of Hans Knibbe one of the first results (July 2011) obtained on the LGM is the solution of the Marmousi problem.

The Helmholtz equation was used to compute the wave propagation in a medium of constant density but variable velocity. About 100 frequencies were computed in parallel using MPI on 15 GPUs (1 GPU per node used) of the LGM. The computation took 30 minutes We used first order radiation boundary conditions.

The Marmousi model can be found at the following site.

The source is a 15Hz Ricker wavelet. The computational grid is the same as the model grid (2301 x 751 grid points) as the grid is fine enough.

Click on the figure below to see the movement of the seismic front in the Marmousi problem.

Layout of the Marmousi problem

The second movie has been produced using the central part of the Marmousi model. The horizontal extents are from 3068 meters to 6136. Like the original model, the vertical extents are from 0 to 3004 meters. The model was resampled so that dx=dy=4/3 meters.

The source was a 15 Hz Ricker wavelet placed in the middle at a depth of 2*4/3=2.6 meters.

Part of the funding for the LGM has been obtained by NWO under the project number 612.071.305.

For background information of the method used in this work we refer to:

H. Knibbe and C.W. Oosterlee and C. Vuik (pdf, bibtex)
GPU implementation of a Helmholtz Krylov solver preconditioned by a shifted Laplace multigrid method
Journal of Computational and Applied Mathematics, 236, pp. 281-293, 2011

A presentation has been given at the ENUMATH Conference 2011.

Contact information: Kees Vuik

Back to the home page of Kees Vuik