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
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
The Marmousi model can be found at the following
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
Part of the funding for the LGM has been obtained by NWO under the project number
For background information of the method used in this work we refer to:
H. Knibbe and C.W. Oosterlee and C. Vuik
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
has been given at the ENUMATH Conference 2011.
Back to the
of Kees Vuik