An algorithm for solving systems of linear algebraic equations based on the Gaussian elimination method is proposed. The algorithm is aimed to solve boundary value problems for the Helmholtz equation in 3D heterogeneous media. In order to solve linear systems raised from geophysical applications, we developed a parallel version targeted on heterogeneous high-performance computing clusters (MPP and SMP architecture). Using the low-rank approximation technique and the HSS format allows us to solve problems larger than by the use of traditional direct solvers with saving the L-factor in full rank (FR). Using the proposed approach reduces computation time; it is the key-point of 3D geophysical problems. Numerical experiments demonstrate a number of advantages of the proposed low-rank approach in comparison with direct solvers (FR-approaches).

Keywords: Helmholtz equation, algorithms for sparse systems of linear algebraic equations, Gaussian elimination method, low-rank approximation, HSS matrix representation, distributed parallel systems, heterogeneous high-performance computing systems.

• PDF (Full text in Russian, References in English)   (519KB)



• PDF. zip (Full text in Russian, References in English)   (452KB)