Full text loading...
-
Cluster Implementation of Low-Rank Multifrontal Direct Solver for 3D Helmholtz Problem
- Publisher: European Association of Geoscientists & Engineers
- Source: Conference Proceedings, 7th EAGE Saint Petersburg International Conference and Exhibition, Apr 2016, cp-480-00035
- ISBN: 978-94-6282-179-8
Abstract
The modern methods of processing the geophysical data, such as Reverse Time Migration (RTM) and Full Waveform Inversion (FWI) require solving series of forward problems where the main step is solution of Systems of Linear Algebraic Equations (SLAE) of big size. For big sizes, it is time and memory consuming problem. In this paper, we present a parallel direct algorithm to solve boundary value problems for 3D Helmholtz equation discretized with help of finite differences. The memory consumption has been resolved due to Nested Dissection approach, low-rank approximation technique and HSS format. OpenMP parallelization is based on standard BLAS and LAPACK functionality. For MPI parallelization, we propose a novel algorithm that uses dynamical distribution of the elimination tree nodes across cluster nodes. Numerical experiments show performance benefits of the proposed cluster algorithm compared to the not parallel version and demonstrate significant memory advantages over direct solvers without low-rank approximation.