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.

Regular polynomial interpolation and approximation of global solutions of linear partial differential equations

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

Sparse-grid polynomial interpolation approximation and integration for parametric and stochastic elliptic PDEs with lognormal inputs

Sparse polynomial approximations for affine parametric saddle point problems

Linear Scaling Algorithms for Solving High-Dimensional Nonlinear Parabolic Differential Equations

High-order empirical interpolation methods for real time solution of parametrized nonlinear PDEs

Construction of polynomial particular solutions of linear constant-coefficient partial 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 Iterations for Solving Smooth Semilinear Parabolic Heat Equations

Adaptive rational Krylov methods for exponential Runge--Kutta integrators

Superconvergent Non-Polynomial Approximations

Generalised multilevel Picard approximations

Solution of Polynomial Systems Derived from Differential Equations

Exploiting locality in sparse polynomial approximation of parametric elliptic PDEs and application to parameterized domains

Optimal approximation of elliptic problems by linear and nonlinear mappings II

Analytic regularity and solution approximation for a semilinear elliptic partial differential equation in a polygon

Jump-preserving polynomial interpolation in non-manifold polyhedra

Regularized dynamical parametric approximation

Optimal error estimates for non-conforming approximations of linear parabolic problems with minimal regularity

Order Preserving Interpolation for Summation-by-Parts Operators at Non-Conforming Grid Interfaces