In Monte Carlo methods quadrupling the sample size halves the error. In simulations of stochastic partial differential equations (SPDEs), the total work is the sample size times the solution cost of an instance of the partial differential equation. A Multi-level Monte Carlo method is introduced which allows, in certain cases, to reduce the overall work to that of the discretization of one instance of the deterministic PDE. The model problem is an elliptic equation with stochastic coefficients. Multi-level Monte Carlo errors and work estimates are given both for the mean of the solutions and for higher moments. The overall complexity of computing mean fields as well as k-point correlations of the random solution is proved to be of log-linear complexity in the number of unknowns of a single Multi-level solve of the deterministic elliptic problem. Numerical examples complete the theoretical analysis.
%0 Journal Article
%1 barth2011multilevel
%A Barth, Andrea
%A Schwab, Christoph
%A Zollinger, Nathaniel
%D 2011
%J Numerische Mathematik
%K ians ians-uq myown
%P 123-161
%R 10.1007/s00211-011-0377-0
%T Multi-level Monte Carlo Finite Element Method for Elliptic PDEs with Stochastic Coefficients.
%U http://dx.doi.org/10.1007/s00211-011-0377-0
%V 119
%X In Monte Carlo methods quadrupling the sample size halves the error. In simulations of stochastic partial differential equations (SPDEs), the total work is the sample size times the solution cost of an instance of the partial differential equation. A Multi-level Monte Carlo method is introduced which allows, in certain cases, to reduce the overall work to that of the discretization of one instance of the deterministic PDE. The model problem is an elliptic equation with stochastic coefficients. Multi-level Monte Carlo errors and work estimates are given both for the mean of the solutions and for higher moments. The overall complexity of computing mean fields as well as k-point correlations of the random solution is proved to be of log-linear complexity in the number of unknowns of a single Multi-level solve of the deterministic elliptic problem. Numerical examples complete the theoretical analysis.
@article{barth2011multilevel,
abstract = {In Monte Carlo methods quadrupling the sample size halves the error. In simulations of stochastic partial differential equations (SPDEs), the total work is the sample size times the solution cost of an instance of the partial differential equation. A Multi-level Monte Carlo method is introduced which allows, in certain cases, to reduce the overall work to that of the discretization of one instance of the deterministic PDE. The model problem is an elliptic equation with stochastic coefficients. Multi-level Monte Carlo errors and work estimates are given both for the mean of the solutions and for higher moments. The overall complexity of computing mean fields as well as k-point correlations of the random solution is proved to be of log-linear complexity in the number of unknowns of a single Multi-level solve of the deterministic elliptic problem. Numerical examples complete the theoretical analysis.},
added-at = {2023-12-01T16:28:04.000+0100},
author = {Barth, Andrea and Schwab, Christoph and Zollinger, Nathaniel},
biburl = {https://puma.ub.uni-stuttgart.de/bibtex/2fac703aaa65db63449604bd71697d6b0/abarth},
coden = {NUMMA7},
doi = {10.1007/s00211-011-0377-0},
fjournal = {Numerische Mathematik},
interhash = {31a6ff3dcb767c423775f6e1abbd123b},
intrahash = {fac703aaa65db63449604bd71697d6b0},
journal = {Numerische Mathematik},
keywords = {ians ians-uq myown},
mrclass = {65N30 (35J25 35R60 60H15 65C05)},
mrnumber = {2824857 (2012i:65252)},
mrreviewer = {Erich Novak},
owner = {barthaa},
pages = {123-161},
timestamp = {2023-12-04T12:53:37.000+0100},
title = {Multi-level Monte Carlo Finite Element Method for Elliptic PDEs with Stochastic Coefficients.},
url = {http://dx.doi.org/10.1007/s00211-011-0377-0},
volume = 119,
year = 2011
}