Generative modelling of Lévy area for high order SDE simulation
File(s)
Accepted version
Author(s)
Type
Journal Article
Abstract
It is well known that, when numerically simulating solutions to SDEs, achieving a strong convergence rate better than O(√h) (where h is the step-size) usually requires the use of
certain iterated integrals of Brownian motion, commonly referred to as its “Lévy areas”. However, these stochastic integrals are difficult to simulate due to their non-Gaussian nature and for a d-dimensional Brownian motion with d > 2, no fast almost-exact sampling algorithm is known. In this paper, we propose LévyGAN, a deep-learning-based model for generating approximate samples of Lévy area conditional on a Brownian increment. Due to our “Bridge-flipping” operation, the output samples match all joint and conditional odd moments exactly. Our generator employs a tailored GNN-inspired architecture, which enforces the correct dependency structure between the output distribution and the conditioning variable. Furthermore, we incorporate a mathematically principled characteristic-function based discriminator. Lastly, we introduce a novel training mechanism termed “Chen-training”, which circumvents the need for expensive-to-generate training data-sets. This new training procedure is underpinned by our two main theoretical results. For 4-dimensional Brownian motion, we show that LévyGAN exhibits state-of-the-art performance across several metrics which measure both the joint and marginal distributions. We conclude with a numerical experiment on the log-Heston model, a popular SDE in mathematical finance, demonstrating that high-quality synthetic Lévy area can lead to high order weak convergence and variance reduction when using multilevel
Monte Carlo (MLMC).
certain iterated integrals of Brownian motion, commonly referred to as its “Lévy areas”. However, these stochastic integrals are difficult to simulate due to their non-Gaussian nature and for a d-dimensional Brownian motion with d > 2, no fast almost-exact sampling algorithm is known. In this paper, we propose LévyGAN, a deep-learning-based model for generating approximate samples of Lévy area conditional on a Brownian increment. Due to our “Bridge-flipping” operation, the output samples match all joint and conditional odd moments exactly. Our generator employs a tailored GNN-inspired architecture, which enforces the correct dependency structure between the output distribution and the conditioning variable. Furthermore, we incorporate a mathematically principled characteristic-function based discriminator. Lastly, we introduce a novel training mechanism termed “Chen-training”, which circumvents the need for expensive-to-generate training data-sets. This new training procedure is underpinned by our two main theoretical results. For 4-dimensional Brownian motion, we show that LévyGAN exhibits state-of-the-art performance across several metrics which measure both the joint and marginal distributions. We conclude with a numerical experiment on the log-Heston model, a popular SDE in mathematical finance, demonstrating that high-quality synthetic Lévy area can lead to high order weak convergence and variance reduction when using multilevel
Monte Carlo (MLMC).
Date Issued
2025-12-01
Date Acceptance
2025-04-11
Citation
SIAM Journal on Mathematics of Data Science, 2025, 7 (4), pp.1541-1567
ISSN
2577-0187
Publisher
Society for Industrial and Applied Mathematics
Start Page
1541
End Page
1567
Journal / Book Title
SIAM Journal on Mathematics of Data Science
Volume
7
Issue
4
Copyright Statement
© 2025 Society for Industrial and Applied Mathematics. This is the author’s accepted manuscript made available under a CC-BY licence in accordance with Imperial’s Research Publications Open Access policy (www.imperial.ac.uk/oa-policy)
License URL
Publication Status
Published
Date Publish Online
2025-10-03