Multilevel Picard iterations for solving smooth semilinear parabolic heat equations

Abstract

We introduce a new family of numerical algorithms for approximating solutions of general high-dimensional semilinear parabolic partial differential equations at single space-time points. The algorithm is obtained through a delicate combination of the Feynman-Kac and the Bismut-Elworthy-Li formulas, and an approximate decomposition of the Picard fixed-point iteration with multilevel accuracy. The algorithm has been tested on a variety of semilinear partial differential equations that arise in physics and finance, with very satisfactory results. Analytical tools needed for the analysis of such algorithms, including a semilinear Feynman-Kac formula, a new class of semi-norms and their recursive inequalities, are also introduced. They allow us to prove for semilinear heat equations with gradient-independent nonlinearity that the computational complexity of the proposed algorithm is bounded by $O(d{\epsilon }^{-(4+\delta )})$ for any $\delta \in (0,\infty )$ under suitable assumptions, where $d\in \mathbb{N}$ is the dimensionality of the problem and $\epsilon \in (0,\infty )$ is the prescribed accuracy.

References in corpus (10)

• A regression-based Monte Carlo method to solve backward stochastic differential equations
• Formulae for the derivatives of heat semigroups
• On multilevel Picard numerical approximations for high-dimensional nonlinear parabolic partial differential equations and high-dimensional nonlinear backward stochastic differential equations
• Statistical Romberg extrapolation: A new variance reduction method and applications to option pricing
• A proof that rectified deep neural networks overcome the curse of dimensionality in the numerical approximation of semilinear heat equations
• On a perturbation theory and on strong convergence rates for stochastic ordinary and partial differential equations with non-globally monotone coefficients
• Overcoming the curse of dimensionality in the numerical approximation of semilinear parabolic partial differential equations
• Overcoming the curse of dimensionality in the approximative pricing of financial derivatives with default risks
• Overcoming the curse of dimensionality in the numerical approximation of Allen-Cahn partial differential equations via truncated full-history recursive multilevel Picard approximations
• Numerical simulation of BSDEs using empirical regression methods: theory and practice

Cited by in corpus (31)

• Solving high-dimensional partial differential equations using deep learning
• Deep learning-based numerical methods for high-dimensional parabolic partial differential equations and backward stochastic differential equations
• Tackling the Curse of Dimensionality with Physics-Informed Neural Networks
• Algorithms for Solving High Dimensional PDEs: From Nonlinear Monte Carlo to Machine Learning
• A proof that artificial neural networks overcome the curse of dimensionality in the numerical approximation of Black-Scholes partial differential equations
• Convergence of the Deep BSDE Method for Coupled FBSDEs
• On multilevel Picard numerical approximations for high-dimensional nonlinear parabolic partial differential equations and high-dimensional nonlinear backward stochastic differential equations
• Deep splitting method for parabolic PDEs
• An overview on deep learning-based approximation methods for partial differential equations
• Overcoming the curse of dimensionality in the approximative pricing of financial derivatives with default risks
• Overcoming the curse of dimensionality in the numerical approximation of Allen-Cahn partial differential equations via truncated full-history recursive multilevel Picard approximations
• Hutchinson Trace Estimation for High-Dimensional and High-Order Physics-Informed Neural Networks
• Numerical simulations for full history recursive multilevel Picard approximations for systems of high-dimensional partial differential equations
• Numerical methods for backward stochastic differential equations: A survey
• On nonlinear Feynman-Kac formulas for viscosity solutions of semilinear parabolic partial differential equations
• Overcoming the curse of dimensionality in the numerical approximation of high-dimensional semilinear elliptic partial differential equations
• Solving non-linear Kolmogorov equations in large dimensions by using deep learning: a numerical comparison of discretization schemes
• Multilevel Picard approximations for McKean-Vlasov stochastic differential equations
• Generalised multilevel Picard approximations
• Learning the random variables in Monte Carlo simulations with stochastic gradient descent: Machine learning for parametric PDEs and financial derivative pricing
• Deep learning algorithms for solving high dimensional nonlinear backward stochastic differential equations
• Multilevel Picard approximations for high-dimensional semilinear second-order PDEs with Lipschitz nonlinearities
• Monte Carlo for high-dimensional degenerated Semi Linear and Full Non Linear PDEs
• Multi-level Picard approximations of high-dimensional semilinear parabolic differential equations with gradient-dependent nonlinearities
• On existence and uniqueness properties for solutions of stochastic fixed point equations
• Full history recursive multilevel Picard approximations for ordinary differential equations with expectations
• Numerical computation of probabilities for nonlinear SDEs in high dimension using Kolmogorov equation
• A learning scheme by sparse grids and Picard approximations for semilinear parabolic PDEs
• Strong $L^p$-error analysis of nonlinear Monte Carlo approximations for high-dimensional semilinear partial differential equations
• Strong convergence rate of Euler-Maruyama approximations in temporal-spatial Hölder-norms
• A Numerical Scheme For High-dimensional Backward Stochastic Differential Equation Based On Modified Multi-level Picard Iteration