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.

Picard Approximation of a Singular Backward Stochastic Nonlinear Volterra Integral Equation

Singular backward stochastic Volterra integral equations in infinite dimensional spaces

Singular mean-field backward stochastic Volterra integral equations in infinite dimensional spaces

Numerical approximation of ergodic BSDEs using non linear Feynman-Kac formulas

On full history recursive multilevel Picard approximations and numerical approximations for high-dimensional nonlinear parabolic partial differential equations and high-dimensional nonlinear backward stochastic differential equations

Extended backward stochastic Volterra integral equations, Quasilinear parabolic equations, and Feynman–Kac formula

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

A Fractional spectral method for weakly singular Volterra integro-differential equations with delays of the third-kind

Backward stochastic Volterra integral equations with time delayed generators

Total Variation Error Bounds for the Approximation of the Invariant Distribution of Parabolic Semilinear SPDEs Using the Standard Euler Scheme

A numerical study on the non‐smooth solutions of the nonlinear weakly singular fractional Volterra integro‐differential equations

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

Multilevel Picard algorithm for general semilinear parabolic PDEs with gradient-dependent nonlinearities

Time discretization of Quadratic Forward-Backward SDEs with singular drifts

Backward Doubly Stochastic Differential Equations with Polynomial Growth Coefficients

Large Deviation Principle for Backward Stochastic Differential Equations with a stochastic Lipschitz condition on $z$

Multivalued stochastic partial differential-integral equations via backward doubly stochastic differential equations driven by a Lévy process

A Monte Carlo method for backward stochastic differential equations with Hermite martingales

A Müntz-collocation spectral method for weakly singular Volterra delay-integro-differential equations

Temporal approximation of stochastic evolution equations with irregular nonlinearities