High performance cluster based low-rank solution of helmholtz problem for applying in full waveform inversion

S. A. Solovyev, V. I. Kostin

Research output: Chapter in Book/Report/Conference proceedingConference contribution


We describe a MPI-based multifrontal direct solver for Helmholtz problem in 3D heterogeneous media. To reduce memory consumption and computational time at the factorization step, intermediate data are compressed with help of Low-Rank approximation and Hierarchically Semi-Separable format. The cluster version is based on performance efficient approach: the factorization of various parts of matrix are distributed through cluster nodes in advance. It allows us highly parallelize the major jobs of factorization process. i.e. low-rank approximation and computing the Schur complement. Such improvements make it possible to factorize in acceptable time (∼1 hour) a system of more than 10 ∧ 8 equations corresponding to a realistic geophysical velocity model on ∼1200km ∧ 2 square. The inversion of factorized system takes less than 1 second per one right hand side. It allows us to use this solver in solving series of forward problems in Full Waveform Inversion (FWI) process where the main step is solution of Systems of Linear Algebraic Equations (SLAE) of big size. Numerical experiments show well scalability on moderate number of cluster nodes and demonstrate high performance and memory compressibility both on homogeneous and high contrast heterogeneous velocity models for various frequencies.
Original languageEnglish (US)
Title of host publicationSaint Petersburg 2018
PublisherEAGE Publications BV
ISBN (Print)9789462822474
StatePublished - Apr 9 2018
Externally publishedYes

Bibliographical note

KAUST Repository Item: Exported on 2022-06-28
Acknowledgements: The research described was partially supported by RFBR grants 16-05-00800, 18-55-15008, 18-01-00225, 18-01-00295 and the Russian Academy of Sciences 3rogram ”Arctic”. We also appreciate KAUST authorities for providing access to Shaheen II supercomputer.
This publication acknowledges KAUST support, but has no KAUST affiliated authors.


Dive into the research topics of 'High performance cluster based low-rank solution of helmholtz problem for applying in full waveform inversion'. Together they form a unique fingerprint.

Cite this