Solving Elliptic Equations with Brownian Motion: Bias Reduction and Temporal Difference Learning

Abstract

The Feynman-Kac formula provides a way to understand solutions to elliptic partial differential equations in terms of expectations of continuous time Markov processes. This connection allows for the creation of numerical schemes for solutions based on samples of these Markov processes which have advantages over traditional numerical methods in some cases. However, naïve numerical implementations suffer from statistical bias and sampling error. We present methods to discretize the stochastic process appearing in the Feynman-Kac formula that reduce the bias of the numerical scheme. We also propose using temporal difference learning to assemble information from random samples in a way that is more efficient than the traditional Monte Carlo method.

arXiv comment18 pages, 6 figures

References in corpus (14)

• Artificial Neural Networks for Solving Ordinary and Partial Differential Equations
• DGM: A deep learning algorithm for solving partial differential equations
• Solving high-dimensional partial differential equations using deep learning
• Hidden Physics Models: Machine Learning of Nonlinear Partial Differential Equations
• Physics-Constrained Deep Learning for High-dimensional Surrogate Modeling and Uncertainty Quantification without Labeled Data
• Deep learning-based numerical methods for high-dimensional parabolic partial differential equations and backward stochastic differential equations
• Deep Hidden Physics Models: Deep Learning of Nonlinear Partial Differential Equations
• Simulator-free Solution of High-Dimensional Stochastic Elliptic Partial Differential Equations using Deep Neural Networks
• A Deep Neural Network Surrogate for High-Dimensional Random Partial Differential Equations
• Solving high-dimensional eigenvalue problems using deep neural networks: A diffusion Monte Carlo like approach
• Hidden Fluid Mechanics: A Navier-Stokes Informed Deep Learning Framework for Assimilating Flow Visualization Data
• A Derivative-Free Method for Solving Elliptic Partial Differential Equations with Deep Neural Networks
• On the Euler discretization error of Brownian motion about random times
• A Path Integral Monte Carlo Method based on Feynman-Kac Formula for Electrical Impedance Tomography

Cited by in corpus (1)

• Learning the random variables in Monte Carlo simulations with stochastic gradient descent: Machine learning for parametric PDEs and financial derivative pricing