Abstract:The Multilevel Monte Carlo (MLMC) method has proven to be an effective variance-reduction statistical method for Uncertainty Quantification (UQ) in Partial Differential Equation (PDE) models, combining model computations at different levels to create an accurate estimate. Still, the computational complexity of the resulting method is extremely high, particularly for 3D models, which requires advanced algorithms for the efficient exploitation of High Performance Computing (HPC). In this article we present a new implementation of the MLMC in massively parallel computer architectures, exploiting parallelism within and between each level of the hierarchy. The numerical approximation of the PDE is performed using the finite element method but the algorithm is quite general and could be applied to other discretization methods as well, although the focus is on parallel sampling. The two key ingredients of an efficient parallel implementation are a good processor partition scheme together with a good scheduling algorithm to assign work to different processors. We introduce a multiple partition of the set of processors that permits the simultaneous execution of different levels and we develop a dynamic scheduling algorithm to exploit it. The problem of finding the optimal scheduling of distributed tasks in a parallel computer is an NP-complete problem. We propose and analyze a new greedy scheduling algorithm to assign samples and we show that it is a 2-approximation, which is the best that may be expected under general assumptions. On top of this result we design a distributed memory implementation using the Message Passing Interface (MPI) standard. Finally we present a set of numerical experiments illustrating its scalability properties.

A Fully Parallelized and Budgeted Multi-Level Monte Carlo Method and the Application to Acoustic Waves

A Fully Parallelized and Budgeted Multilevel Monte Carlo Method and the Application to Acoustic Waves

A Fully Parallelizable Space-Time Multilevel Monte Carlo Method for Stochastic Differential Equations with Additive Noise

$p$-multilevel Monte Carlo for acoustic scattering from large deviation rough random surfaces

High-performance modeling acoustic and elastic waves using the Parallel Dichotomy Algorithm

A Massively Parallel Implementation of Multilevel Monte Carlo for Finite Element Models

Scheduling massively parallel multigrid for multilevel Monte Carlo methods

A performance analysis of a mimetic finite difference scheme for acoustic wave propagation on GPU platforms

Multilevel Fast Multipole Boundary Element Method for 2D Acoustic Problems and Its Applications

A Parallel 2D/3D Aero-acoustics Numerical Algorithm and Implementation

Parallel Numerical Simulation of Weakly Range-Dependent Ocean Acoustic Waveguides by Adiabatic Modes Based on a Spectral Method

Parallel WNAD Method and Its Wave Field Simulations

Parallel implementation for three-dimensional acoustic field computation in a penetrable wedge by image source method

Multicore Parallel Programming of Acoustic Equation Numerical Modeling

Application of MPI in forward modeling of acoustic equation.

Dynamic hybrid parallel computing of the Ray Model for solving underwater acoustic fields in vast sea

A Low Overhead Heterogeneous Parallel Optimization Method Based on 3-D Elastic Wave Numerical Simulation

Massively parallel nodal discontinous Galerkin finite element method simulator for room acoustics

A Low Overhead Heterogeneous Parallel Optimization Method Based on Three-Dimensional Elastic Wave Numerical Simulation

Performance Optimization and Parallelization of a Parabolic Equation Solver in Computational Ocean Acoustics on Modern Many-core Computer

Parallelization of Kinetic Theory Simulations