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.

Parametric PDEs: Sparse or Low-Rank Approximations?

Sparse polynomial approximations for affine parametric saddle point problems

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

A space-time adaptive low-rank method for high-dimensional parabolic 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

Expression Rates of Neural Operators for Linear Elliptic PDEs in Polytopes

Adaptive Low-Rank Methods: Problems on Sobolev Spaces

Nonlinear compressive reduced basis approximation for multi-parameter elliptic problem

Measurability and continuity of parametric low-rank approximation in Hilbert spaces: linear operators and random variables

A Dimension-adaptive Combination Technique for Uncertainty Quantification

Quasi-Monte Carlo for Bayesian design of experiment problems governed by parametric PDEs

Dynamical low-rank tensor approximations to high-dimensional parabolic problems: existence and convergence of spatial discretizations

Nonparametric Estimation via Partial Derivatives

Tensor-Sparsity of Solutions to High-Dimensional Elliptic Partial Differential Equations

Multi-level Compressed Sensing Petrov-Galerkin discretization of high-dimensional parametric PDEs

Randomized quasi-optimal local approximation spaces in time

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

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

Optimal approximation of elliptic problems by linear and nonlinear mappings II

Goal-oriented error estimation and adaptivity for elliptic PDEs with parametric or uncertain inputs