Monte Carlo versus multilevel Monte Carlo in weak error simulations of SPDE approximations

Abstract

The simulation of the expectation of a stochastic quantity E[Y] by Monte Carlo methods is known to be computationally expensive especially if the stochastic quantity or its approximation Yn is expensive to simulate, e.g., the solution of a stochastic partial differential equation. If the convergence of Yn to Y in terms of the error |E[Y - Yn]| is to be simulated, this will typically be done by a Monte Carlo method, i.e., |E[Y] - EN[Yn]| is computed. In this article upper and lower bounds for the additional error caused by this are determined and compared to those of |EN[Y - Yn]|, which are found to be smaller. Furthermore, the corresponding results for multilevel Monte Carlo estimators, for which the additional sampling error converges with the same rate as |E[Y - Yn]|, are presented. Simulations of a stochastic heat equation driven by multiplicative Wiener noise and a geometric Brownian motion are performed which confirm the theoretical results and show the consequences of the presented theory for weak error simulations.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…