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 test of backward stochastic differential equations solver for solving semilinear parabolic differential equations in 1D and 2D

A backward differential deep learning-based algorithm for solving high-dimensional nonlinear backward stochastic differential equations

A forward differential deep learning-based algorithm for solving high-dimensional nonlinear backward stochastic differential equations

Deep Learning-Based Numerical Methods for High-Dimensional Parabolic Partial Differential Equations and Backward Stochastic Differential Equations

Numerical methods for backward stochastic differential equations: A survey

Gradient Convergence of Deep Learning-Based Numerical Methods for BSDEs

Convergence of the Deep BSDE Method for Coupled FBSDEs

Deep multi-step mixed algorithm for high dimensional non-linear PDEs and associated BSDEs

Quasi-Regression Monte-Carlo scheme for semi-linear PDEs and BSDEs with large scale parallelization on GPUs

The One Step Malliavin scheme: new discretization of BSDEs implemented with deep learning regressions

Towards Robust and Stable Deep Learning Algorithms for Forward Backward Stochastic Differential Equations

SEMI-LINEAR SYSTEMS OF BACKWARD STOCHASTIC PARTIAL DIFFERENTIAL EQUATIONS IN ℝn

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

FBSDE based Neural Network Algorithms for High-Dimensional Quasilinear Parabolic PDEs

On the Homotopy Analysis Method for Backward/forward-Backward Stochastic Differential Equations

Convergence of the Backward Deep BSDE Method with Applications to Optimal Stopping Problems

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

Local discontinuous Galerkin method for nonlinear BSPDEs of Neumann boundary conditions with deep backward dynamic programming time-marching

Deep Operator BSDE: a Numerical Scheme to Approximate the Solution Operators

Gradient boosting-based numerical methods for high-dimensional backward stochastic differential equations

A Convolution Method for Numerical Solution of Backward Stochastic Differential Equations Based on the Fractional FFT