pith. sign in

arxiv: 1806.05619 · v1 · pith:OFLMQQ7Gnew · submitted 2018-06-14 · 🧮 math.NA

Recycling Samples in the Multigrid Multilevel (Quasi-)Monte Carlo Method

classification 🧮 math.NA
keywords methodcarlomontemultilevelrandomwhenapproachcoarse
0
0 comments X
read the original abstract

The Multilevel Monte Carlo method is an efficient variance reduction technique. It uses a sequence of coarse approximations to reduce the computational cost in uncertainty quantification applications. The method is nowadays often considered to be the method of choice for solving PDEs with random coefficients when many uncertainties are involved. When using Full Multigrid to solve the deterministic problem, coarse solutions obtained by the solver can be recycled as samples in the Multilevel Monte Carlo method, as was pointed out by Kumar, Oosterlee and Dwight [Int. J. Uncertain. Quantif., 7 (2017), pp. 57--81]. In this article, an alternative approach is considered, using Quasi-Monte Carlo points, to speed up convergence. Additionally, our method comes with an improved variance estimate which is also valid in case of the Monte Carlo based approach. The new method is illustrated on the example of an elliptic PDE with lognormal diffusion coefficient. Numerical results for a variety of random fields with different smoothness parameters in the Mat\'ern covariance function show that sample recycling is more efficient when the input random field is nonsmooth.

This paper has not been read by Pith yet.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.