Abstract
The parallel performance of a numerical solution method for the scalar 2D Helmholtz equation written for inhomogeneous media is studied. The numerical solution is obtained by an iterative method applied to the preconditioned linear system which has been derived from a finite difference discretization. The preconditioner is approximately inverted using multigrid iterations. Parallel execution is implemented using the MPI library. Only a few iterations are required to solve numerically the so-called full Marmousi problem [Bourgeois, A., et al. in The Marmousi Experience, Proceedings of the 1990 EAEG Workshop on Practical Aspects of Seismic Data Inversion: Eur. Assoc. Expl. Geophys., pp. 5-16 (1991)] for the high frequency range.
Original language | English |
---|---|
Pages (from-to) | 139-146 |
Number of pages | 8 |
Journal | Computing and Visualization in Science |
Volume | 11 |
Issue number | 3 |
DOIs | |
Publication status | Published - May 2008 |
Externally published | Yes |