Efficient almost-exact Lévy area sampling

Simon J A Malham*, Anke Wiese

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

15 Citations (Scopus)


We present a new method for sampling the Lévy area for a two-dimensional Wiener process conditioned on its endpoints. An efficient sampler for the Lévy 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 Lévy area by an infinite weighted sum of independent Logistic random variables. We use Chebyshev 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 Lévy area sampling method. The complexity of our method is square logarithmic. We indicate how it can contribute to efficient sampling in higher dimensions.

Original languageEnglish
Pages (from-to)50-55
Number of pages6
JournalStatistics and Probability Letters
Issue number1
Publication statusPublished - 2014


  • Chebyshev approximation
  • Lévy area
  • Logistic expansion
  • Milstein method
  • Strong simulation

ASJC Scopus subject areas

  • Statistics, Probability and Uncertainty
  • Statistics and Probability


Dive into the research topics of 'Efficient almost-exact Lévy area sampling'. Together they form a unique fingerprint.

Cite this