Title
Recycling Samples in the Multigrid Multilevel (Quasi-)Monte Carlo Method
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 Matern covariance function show that sample recycling is more efficient when the input random field is nonsmooth.
Year
DOI
Venue
2019
10.1137/18M1194031
SIAM JOURNAL ON SCIENTIFIC COMPUTING
Keywords
Field
DocType
multilevel Monte Carlo,multigrid,PDEs with random coefficients
Applied mathematics,Covariance function,Mathematical optimization,Monte Carlo method,Random field,Uncertainty quantification,Quasi-Monte Carlo method,Solver,Variance reduction,Multigrid method,Mathematics
Journal
Volume
Issue
ISSN
41
5
1064-8275
Citations 
PageRank 
References 
0
0.34
0
Authors
3
Name
Order
Citations
PageRank
Pieterjan Robbe110.71
Dirk Nuyens216817.97
Stefan Vandewalle350162.63