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.

Method Monte-Carlo for solving of non-linear integral equations

An unbiased Monte Carlo method to solve linear Volterra equations of the second kind

A simple Monte Carlo method for solving of Navier-Stokes Equations

An iterative Monte Carlo method to solve nonlinear second-order differential equations

Finite Element Error Analysis of Elliptic PDEs with Random Coefficients and Its Application to Multilevel Monte Carlo Methods

A global approximation method for second-kind nonlinear integral equations

Numerical Solution of Multiple Nonlinear Volterra Integral Equations

A meshless local Galerkin method for solving a class of nonlinear time-dependent mixed integral equations on non-rectangular 2D domains

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

Application of Monte-Carlo method to reactor calculation

An extended study of two multistep higher order convergent methods for solving nonlinear equations

Nonasymptotic Convergence Rate of Quasi-Monte Carlo: Applications to Linear Elliptic PDEs with Lognormal Coefficients and Importance Samplings

Numerical solutions of a class of linear and nonlinear Volterra integral equations of the third kind using collocation method based on radial basis functions

Monte Carlo method and the random isentropic Euler system

Uniform convergence guarantees for the deep Ritz method for nonlinear problems

Asymptotic properties of Monte Carlo methods in elliptic PDE-constrained optimization under uncertainty

Reliable numerical solution of a class of nonlinear elliptic problems generated by the Poisson-Boltzmann equation

Euler-Maruyama approximation for stochastic fractional neutral integro-differential equations with weakly singular kernel

Multilevel Monte Carlo methods for stochastic convection-diffusion eigenvalue problems

Numerical Solution of Nonclassical Boundary Value Problems

Goal-oriented adaptive finite element multilevel Monte Carlo with convergence rates