Parallel algorithm for 3D modeling of monochromatic acoustic field using integral equations

Mikhail S. Malovichko, Nikolay I. Khokhlov, Nikolay B. Yavich, Michael S. Zhdanov

Research output: Contribution to journalArticlepeer-review

1 Citation (Scopus)

Abstract

We present a parallel algorithm for solution of the three-dimensional Helmholtz equation in the frequency domain using volume-integral equations. The algorithm is applied to seismic forward modeling. The method of integral equations reduces the size of the problem by dividing the geologic model into the anomalous and background parts, but leads to a dense system matrix. Tolerable memory consumption and numerical complexity were achieved by applying an iterative solver, accompanied by an effective matrix-vector multiplication operation, based on the fast Fourier transform. We used OpenMP to speed up the matrix-vector multiplication, while MPI was used to speed up the equation system solver, and also for parallelizing across multiple sources. Practical examples and efficiency tests are presented.

Original languageEnglish
Pages (from-to)74-78
Number of pages5
JournalSupercomputing Frontiers and Innovations
Volume3
Issue number4
DOIs
Publication statusPublished - 2016
Externally publishedYes

Keywords

  • Acoustics
  • Integral equations
  • MPI
  • OpenMP
  • Seismics

Fingerprint

Dive into the research topics of 'Parallel algorithm for 3D modeling of monochromatic acoustic field using integral equations'. Together they form a unique fingerprint.

Cite this