Abstract:Parabolic partial differential equations (PDEs) and backward stochastic differential equations (BSDEs) are key ingredients a number of models physics and financial engineering. In particular, parabolic PDEs and BSDEs are fundamental tools the state-of-the-art pricing and hedging of financial derivatives. The PDEs and BSDEs appearing 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 this paper we introduce a family of new approximation methods for high-dimensional PDEs and BSDEs. A key idea of our methods is to combine multilevel approximations with Picard fixed-point approximations. Thereby we obtain a class of full history recursive multilevel Picard approximations. Our error analysis proves for one of the proposed approximation methods that if $epsilonin(0,infty)$ is the prescribed approximation accuracy and if $dinmathbb{N}$ is the dimension of the considered PDE, then for every $delta in (0,infty)$ it holds that the computational effort of the approximation method (number of function evaluations of the coefficient functions of the considered PDE and number of used independent scalar standard normal random variables) is at most $O(d,epsilon^{-(4+delta)})$, that is, the computational effort grows only linearly the dimension and up to an arbitrarily small order quartically the reciprocal of the prescribed approximation accuracy. We illustrate the efficiency of one of the proposed approximation methods by means of numerical simulations presenting approximation accuracy against runtime for several nonlinear PDEs from physics (such as the Allen-Cahn equation) and financial engineering (such as derivative pricing incorporating default risks) the case of $d=100$ space dimensions.

Generalized Dimension Truncation Error Analysis for High-Dimensional Numerical Integration: Lognormal Setting and Beyond

Generalized dimension truncation error analysis for high-dimensional numerical integration: lognormal setting and beyond

Application of dimension truncation error analysis to high-dimensional function approximation in uncertainty quantification

Engineering Analysis Error Estimation when Removing Finite-Sized Features in Nonlinear Elliptic Problems.

Quasi-Monte Carlo for partial differential equations with generalized Gaussian input uncertainty

Finite Element Error Analysis of Elliptic PDEs with Random Coefficients and Its Application to Multilevel Monte Carlo Methods

A Dimension-adaptive Combination Technique for Uncertainty Quantification

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

Analyticity and sparsity in uncertainty quantification for PDEs with Gaussian random field inputs

Optimal randomized changing dimension algorithms for infinite-dimensional integration on function spaces with ANOVA-type decomposition

Truncation Error Analysis for an Accurate Nonlocal Manifold Poisson Model with Dirichlet Boundary

The TRUNC element in any dimension and application to a modified Poisson equation

Constructive Approximation of High-Dimensional Functions with Small Efficient Dimension with Applications in Uncertainty Quantification

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

Nonasymptotic Convergence Rate of Quasi-Monte Carlo: Applications to Linear Elliptic PDEs with Lognormal Coefficients and Importance Samplings

Variational Monte Carlo - Bridging Concepts of Machine Learning and High Dimensional Partial Differential Equations

Quasi-Monte Carlo and discontinuous Galerkin

Error Analysis of Kernel/GP Methods for Nonlinear and Parametric PDEs

On the convergence of adaptive Galerkin FEM for parametric PDEs with lognormal coefficients