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.

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

On Multiscale Methods in Petrov-Galerkin formulation

Multilevel higher order Quasi-Monte Carlo Bayesian Estimation

Sparse Gaussian processes for solving nonlinear PDEs

Parametric PDEs: Sparse or Low-Rank Approximations?

Multilevel CNNs for Parametric PDEs

Hierarchical Super-Localized Orthogonal Decomposition Method

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

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

Nonlinear discontinuous Petrov-Galerkin methods

Nonlinear compressive reduced basis approximation for multi-parameter elliptic problem

Multistage DPG time-marching scheme for nonlinear problems

Chebyshev HOPGD with sparse grid sampling for parameterized linear systems

A One-level Additive Schwarz Preconditioner for a Discontinuous Petrov-Galerkin Method

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

Randomized neural network with Petrov–Galerkin methods for solving linear and nonlinear partial differential equations

Encapsulated PGD Algebraic Toolbox Operating with High-Dimensional Data

Local surrogate models with reduced dimensionality via overlapping domain decomposition and proper generalized decomposition

Natural Factor Based Solvers

Analysis of the SUPG Method for the Solution of Optimal Control Problems

Eliminating the pollution effect in Helmholtz problems by local subscale correction