Efficient almost-exact Lévy area sampling
From MaRDI portal
Publication:2453869
Abstract: We present a new method for sampling the Levy area for a two-dimensional Wiener process conditioned on its endpoints. An efficient sampler for the Levy area is required to implement a strong Milstein numerical scheme to approximate the solution of a stochastic differential equation driven by a two-dimensional Wiener process whose diffusion vector fields do not commute. Our method is simple and complementary to those of Gaines-Lyons and Wiktorsson, and amenable to quasi-Monte--Carlo implementation. It is based on representing the Levy area by an infinite weighted sum of independent Logistic random variables. We use Chebychev polynomials to approximate the inverse distribution function of sums of independent Logistic random variables in three characteristic regimes. The error is controlled by the degree of the polynomials, we set the error to be uniformly 10^(-12). We thus establish a strong almost-exact Levy area sampling method. The complexity of our method is square logarithmic. We indicate how our method can contribute to efficient sampling in higher dimensions.
Recommendations
- Complexity and effective dimension of discrete Lévy areas
- The multilevel Monte Carlo method used on a Lévy driven SDE
- On the approximation and simulation of iterated stochastic integrals and the corresponding Lévy areas in terms of a multidimensional Brownian motion
- Stratified sampling and quasi-Monte Carlo simulation of Lévy processes
- New methods of simulating Lévy processes
Cites work
- Aspects of Brownian motion
- Chi-square simulation of the CIR process and the Heston model
- scientific article; zbMATH DE number 6490089 (Why is no real title available?)
- scientific article; zbMATH DE number 1409619 (Why is no real title available?)
- scientific article; zbMATH DE number 236854 (Why is no real title available?)
- Joint characteristic function and simultaneous simulation of iterated Itô integrals for multiple independent Brownian motions
- On an infinite integral arising in the numerical integration of stochastic differential equations
- On the simulation of iterated Itô integrals.
- Random Generation of Stochastic Area Integrals
- The approximation of multiple stochastic integrals
Cited in
(16)- On the implementation of multilevel Monte Carlo simulation of the stochastic volatility and interest rate model using multi-GPU clusters
- Stability and error analysis of an implicit Milstein finite difference scheme for a two-dimensional Zakai SPDE
- Higher strong order methods for linear Itô SDEs on matrix Lie groups
- An analysis of approximation algorithms for iterated stochastic integrals and a Julia and \textsc{Matlab} simulation toolbox
- Splitting integrators for the stochastic Landau-Lifshitz equation
- Series Expansions and Direct Inversion for the Heston Model
- scientific article; zbMATH DE number 7627965 (Why is no real title available?)
- On the approximation and simulation of iterated stochastic integrals and the corresponding Lévy areas in terms of a multidimensional Brownian motion
- Algebraic structures and stochastic differential equations driven by Lévy processes
- scientific article; zbMATH DE number 7644507 (Why is no real title available?)
- Mean-square approximation of iterated Ito and Stratonovich stochastic integrals: method of generalized multiple Fourier series. Application to numerical integration of Ito SDEs and semilinear SPDEs
- A new proof of the expansion of iterated Itô stochastic integrals with respect to the components of a multidimensional Wiener process based on generalized multiple Fourier series and Hermite polynomials
- Strong convergence of an adaptive time-stepping Milstein method for SDEs with monotone coefficients
- A new approach to the series expansion of iterated Stratonovich stochastic integrals with respect to components of a multidimensional Wiener process. The case of arbitrary complete orthonormal systems in Hilbert space
- A new approach to the series expansion of iterated Stratonovich stochastic integrals with respect to components of a multidimensional Wiener process. The case of arbitrary complete orthonormal systems in Hilbert space. II
- Complexity and effective dimension of discrete Lévy areas
This page was built for publication: Efficient almost-exact Lévy area sampling
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2453869)
