Research reports

Multi-Level Monte Carlo Finite Element Method for elliptic PDEs with stochastic coefficients

by A. Barth and Ch. Schwab and N. Zollinger

(Report number 2010-18)

Abstract
It is a well-known property of Monte Carlo methods that quadrupling the sample size halves the error. In the case of simulations of a stochastic partial differential equations, this implies that the total work is the sample size times the discretization costs of the equation. This leads to a convergence rate which is impractical for many simulations, namely in finance, physics and geosciences. With the Multi--level Monte Carlo method introduced herein, the overall work can be reduced to that of the discretization of the equation, which results in the same convergence rate as for the standard Monte Carlo method. 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. Numerical examples complete the theoretical analysis.

Keywords: Multi-Level Monte Carlo, Stochastic partial differential equations, stochastic finite element methods, multilevel approximations

BibTeX
@Techreport{BSZ10_52,
  author = {A. Barth and Ch. Schwab and N. Zollinger},
  title = {Multi-Level Monte Carlo Finite Element Method for elliptic PDEs with stochastic coefficients},
  institution = {Seminar for Applied Mathematics, ETH Z{\"u}rich},
  number = {2010-18},
  address = {Switzerland},
  url = {https://www.sam.math.ethz.ch/sam_reports/reports_final/reports2010/2010-18.pdf },
  year = {2010}
}

Disclaimer
© Copyright for documents on this server remains with the authors. Copies of these documents made by electronic or mechanical means including information storage and retrieval systems, may only be employed for personal use. The administrators respectfully request that authors inform them when any paper is published to avoid copyright infringement. Note that unauthorised copying of copyright material is illegal and may lead to prosecution. Neither the administrators nor the Seminar for Applied Mathematics (SAM) accept any liability in this respect. The most recent version of a SAM report may differ in formatting and style from published journal version. Do reference the published version if possible (see SAM Publications).

JavaScript has been disabled in your browser