Abstract:Parabolic partial differential equations (PDEs) and backward stochastic differential equations (BSDEs) are key ingredients in a number of models in physics and financial engineering. In particular, parabolic PDEs and BSDEs are fundamental tools in pricing and hedging models for financial derivatives. The PDEs and BSDEs appearing in such applications are often high-dimensional and nonlinear. Since explicit solutions of such PDEs and BSDEs are typically not available, it is a very active topic of research to solve such PDEs and BSDEs approximately. In the recent article (E et al., Multilevel Picard iterations for solving smooth semilinear parabolic heat equations, arXiv:1607.03295) we proposed a family of approximation methods based on Picard approximations and multilevel Monte Carlo methods and showed under suitable regularity assumptions on the exact solution of a semilinear heat equation that the computational complexity is bounded by O(dε-(4+δ))\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$O( d \, {\varepsilon }^{-(4+\delta )})$$\end{document} for any δ∈(0,∞)\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$\delta \in (0,\infty )$$\end{document} where d is the dimensionality of the problem and ε∈(0,∞)\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$${\varepsilon }\in (0,\infty )$$\end{document} is the prescribed accuracy. In this paper, we test the applicability of this algorithm on a variety of 100-dimensional nonlinear PDEs that arise in physics and finance by means of numerical simulations presenting approximation accuracy against runtime. The simulation results for many of these 100-dimensional example PDEs are very satisfactory in terms of both accuracy and speed. Moreover, we also provide a review of other approximation methods for nonlinear PDEs and BSDEs from the scientific literature.

A Fully Parallelizable Space-Time Multilevel Monte Carlo Method for Stochastic Differential Equations with Additive Noise

Monte-Carlo Simulation of a Stochastic Differential Equation

An Antithetic Multilevel Monte Carlo-Milstein Scheme for Stochastic Partial Differential Equations

Stochastic methods for solving high-dimensional partial differential equations

Scheduling massively parallel multigrid for multilevel Monte Carlo methods

Multilevel Space-Time Additive Schwarz Methods for Parabolic Equations

A Fully Parallelized and Budgeted Multilevel Monte Carlo Method and the Application to Acoustic Waves

Parallel-in-time Multilevel Krylov Methods: A Prototype

A Massively Parallel Implementation of Multilevel Monte Carlo for Finite Element Models

Multilevel Monte Carlo methods for simulating forward-backward stochastic differential equations using neural networks

An Adaptive Random Bit Multilevel Algorithm for SDEs

Multilevel Monte Carlo EM scheme for MV-SDEs with small noise

A Fully Parallelized and Budgeted Multi-Level Monte Carlo Method and the Application to Acoustic Waves

A Multilevel Monte Carlo Ensemble and Hybridizable Discontinuous Galerkin Method for a Stochastic Parabolic Problem

Scalable parallel multilevel algorithms for solving partial differential equations

A multilevel Monte Carlo algorithm for SDEs driven by countably dimensional Wiener process and Poisson random measure

Analysis of a New Space-Time Parallel Multigrid Algorithm for Parabolic Problems

A backward Monte-Carlo method for solving parabolic partial differential equations

Multilevel Monte Carlo methods for stochastic convection-diffusion eigenvalue problems

On Multilevel Picard Numerical Approximations for High-Dimensional Nonlinear Parabolic Partial Differential Equations and High-Dimensional Nonlinear Backward Stochastic Differential Equations

Multilevel Monte Carlo methods for the Dean-Kawasaki equation from Fluctuating Hydrodynamics