Massively parallel stochastic solution of the geophysical gravity problem

TL;DR

This paper presents a stochastic algorithm for solving the geophysical gravity problem that is highly parallelizable, mesh-free, and effective for mineral exploration, providing accurate results as an alternative to traditional methods.

Contribution

The paper introduces a stochastic, mesh-free method for solving the Poisson equation in gravity problems, suitable for parallel computing environments.

Findings

Accurately computes gravitational potential and acceleration.

Demonstrates suitability for multi-core and GPU clusters.

Offers an efficient alternative to deterministic methods.

Abstract

In this paper, we report the advantages of using a stochastic algorithm in the context of mineral exploration based on gravity measurements. This approach has the advantage over deterministic methods in that it allows one to find the solution of the Poisson equation in specified, isolated points without the need of meshing the computational domain and solving the Poisson equation over the entire domain. Moreover, the stochastic approach is embarrassingly parallelizable and therefore suitable for an implementation on multi-core compute clusters with or without GPUs. Benchmark tests are carried out that show that the stochastic approach can yield accurate results for both the gravitational potential and the gravitational acceleration and could hence provide an alternative to existing deterministic methods used in mineral exploration.