Memory efficient algorithm for solving the inverse gravimetry problem of finding several boundary surfaces in multilayered medium

Elena Nikolaevna Akimova, Vladimir Evgenevich Misilov, Murat A. Sultanov, Rauan Zh. Turebekov

2019

Closed access Metadata (CSL-JSON)
DOI 10.1063/1.5136168
Data confidence 70 medium
How this score is calculated

It is built from checkable properties of the record itself: whether the article has a DOI, how many independent sources confirmed its data, whether the required fields are filled in, and how recent the latest confirmation is. Disagreement between sources and a missing identifier lower it. A value of 80 or above means the record was published without further checking; 60 to 79 means it is published at reduced confidence, with a manual check advisable. A record verified by an editor by hand is never pushed below 90 by a recalculation. This scores our data about the article, not the article itself.

Abstract

For solving the inverse gravimetry problem of finding several boundary surfaces in a multilayered medium, the parallel algorithm was constructed and implemented for multicore CPU using OpenMP technology. The algorithm is based on the modified nonlinear conjugate gradient method with weighting factors previously proposed by authors. To reduce the memory requirements and computation time, the modification was constructed on the basis of utilizing the Toeplitz-block-Toeplitz structure of the Jacobian matrix of the integral operator. The model problem of reconstructing three surfaces using the quasi-real gravitational data was solved on a large grid. It was shown that the proposed implementation reduces the computation time by 80% in comparison with the earlier algorithm based on calculating the entire matrix. The parallel algorithm shows good scaling of 94% on 8-core processor.