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

• Published in: Proceedings, Saint Petersburg 2018
• Main Authors: Solovyev, S. A., Kostin, V. I.
• Other Authors: IPGG SB RAS, Russian Federation
• Format: Conference Object
• Language: unknown
• Published: EAGE Publications BV 2018
• Subjects: Arctic
• Online Access: http://hdl.handle.net/10754/679403; https://doi.org/10.3997/2214-4609.201800210

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. 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.