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.

Approximation of arbitrarily high-order PDEs by first-order hyperbolic relaxation

A relaxation approach to modeling properties of hyperbolic-parabolic type models

Approximation of Solution Operators for High-dimensional PDEs

Finite difference approximations for the first-order hyperbolic partial differential equation with point-wise delay

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

Approximating Continuous Spectra of Hyperbolic Systems with Summation-by-Parts Finite Difference Operators

First-order greedy invariant-domain preserving approximation for hyperbolic problems: scalar conservation laws, and p-system

Basic Aspects of Hyperbolic Relaxation Systems

Numerical schemes for hyperbolic conservation laws with stiff relaxation terms

A hyperbolic relaxation system of the incompressible Navier-Stokes equations with artificial compressibility

New High-Order Numerical Methods for Hyperbolic Systems of Nonlinear PDEs with Uncertainties

Frequency domain approach for the stability analysis of a fast hyperbolic PDE coupled with a slow ODE

Hyperbolic extensions of constrained PDEs

RELAXATION LIMIT AND INITIAL-LAYERS FOR A CLASS OF HYPERBOLIC-PARABOLIC SYSTEMS

Singular Perturbations of First-Order Hyperbolic Systems with Stiff Source Terms

Large time asymptotics for partially dissipative hyperbolic systems without Fourier analysis: application to the nonlinearly damped p-system

Weakly Nonlinear-Dissipative Approximations of Hyperbolic-Parabolic Systems with Entropy

A dual-pairing summation-by-parts finite difference framework for nonlinear conservation laws

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

Relaxation Methods for Hyperbolic PDE Mixed-Integer Optimal Control Problems

Fast Numerical Approximation of Linear, Second-Order Hyperbolic Problems Using Model Order Reduction and the Laplace Transform