We present a systematic approach to achieving computationally scalable parallelism in the context of geostatistical inversion. Our experience shows that efficient utilization of hardware requires recursive application of domain decomposition, ranging from multi-process model on a cluster of workstations to multithreading on individual CPU cores. Actual run-times and the degree of realistically achievable scalable parallelism depend on a multitude of factors. The list includes the project area size, the probabilistic model complexity, coarse scale of a stochastic iteration of the Multigrid Monte Carlo scheme, and hardware specifications, to name a few. However, overall, we were able to achieve close to optimal multithreaded scalability (where theoretically possible) on up to 32 CPU cores. We also observe that efficiency of multithreading becomes limited by the memory bus bandwidth at some point.


Article metrics loading...

Loading full text...

Full text loading...


  1. Contreras, A., Torres-Verdin, C., Kvien, K., Fasnacht, T. and Chesters, W.
    [2005] AVA Stochastic Inversion of Pre-Stack Seismic Data and Well Logs for 3D Reservoir Modeling.67th EAGE Conference & Exhibition.
    [Google Scholar]
  2. Gilks, W. R., Richardson, S. and Spiegelhalter, D.
    , [1996] Markov Chain Monte Carlo in practice.Chapman & Hall/CRC Interdisciplinary Statistics.
    [Google Scholar]
  3. Goodman, J. and Sokal, A.
    [1989] Multigrid Monte Carlo method. Conceptual foundations.Physical Review D.40(6), 2035–2071.
    [Google Scholar]
  4. Sutter, H.
    [2013] GotW #6a Solution: Const-Correctness, Part 1. https://herbsutter.com/2013/05/24/gotw-6a-const-correctness-part-1-3/
    [Google Scholar]
  5. Tarantola, A.
    [2005] Inverse Problem Theory and Methods for Model Parameter Estimation.Society for Industrial and Applied Mathematics
    [Google Scholar]

Data & Media 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