1887

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.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.201600118
2016-04-11
2024-03-29
Loading full text...

Full text loading...

http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.201600118
Loading
This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error