Related papers: Implementable coupling of L\'evy process and Brown…
In this paper, we show an approximation in law of the complex Brownian motion by processes constructed from a stochastic process with independent increments. We give sufficient conditions for the characteristic function of the process with…
We construct optimal Markov couplings of L\'{e}vy processes, whose L\'evy (jump) measure has an absolutely continuous component. The construction is based on properties of subordinate Brownian motions and the coupling of Brownian motions by…
We introduce two general non-parametric methods for recovering paths of the Brownian and jump components from high-frequency observations of a L\'evy process. The first procedure relies on reordering of independently sampled normal…
The pricing of options in exponential Levy models amounts to the computation of expectations of functionals of Levy processes. In many situations, Monte-Carlo methods are used. However, the simulation of a Levy process with infinite Levy…
We present an iterative sampling method which delivers upper and lower bounding processes for the Brownian path. We develop such processes with particular emphasis on being able to unbiasedly simulate them on a personal computer. The…
We study the error in approximating the minimum of a Brownian motion on the unit interval based on finitely many point evaluations. We construct an algorithm that adaptively chooses the points at which to evaluate the Brownian path. In…
It is known that the Brownian bridge or L\'evy-Ciesielski construction of Brownian paths almost surely converges uniformly to the true Brownian path. In the present article the focus is on the uniform error. In particular, we show…
We characterize the small-time asymptotic behavior of the exit probability of a L\'evy process out of a two-sided interval and of the law of its overshoot, conditionally on the terminal value of the process. The asymptotic expansions are…
It is shown how to construct a successful co-adapted coupling of two copies of an $n$-dimensional Brownian motion $(B_1,...,B_n)$ while simultaneously coupling all corresponding copies of L\'{e}vy stochastic areas $\int B_i dB_j-\int B_j…
A careful look at rough path topology applied to Brownian motion reveals new possible properties of the well-known L\'evy area, in particular the presence of an intrinsic drift of this area. Using renormalization limit of Markov chains on…
We consider a general class of high order weak approximation schemes for stochastic differential equations driven by L\'evy processes with infinite activity. These schemes combine a compound Poisson approximation for the jump part of the…
We apply multilevel Monte Carlo for option pricing problems using exponential L\'{e}vy models with a uniform timestep discretisation to monitor the running maximum required for lookback and barrier options. The numerical results demonstrate…
L\'{e}vy processes with completely monotone jumps appear frequently in various applications of probability. For example, all popular stock price models based on L\'{e}vy processes (such as the Variance Gamma, CGMY/KoBoL and Normal Inverse…
We develop a computational method for expected functionals of the drawdown and its duration in exponential L\'evy models. It is based on a novel simulation algorithm for the joint law of the state, supremum and time the supremum is attained…
We introduce modifications to Monte Carlo simulations of the Feynman path integral that improve sampling of localised interactions. The new algorithms generate trajectories in simple background potentials designed to concentrate them about…
We introduce an algorithm for the pricing of finite expiry American options driven by L\'evy processes. The idea is to tweak Carr's `Canadisation' method, cf. Carr [9] (see also Bouchard et al [5]), in such a way that the adjusted algorithm…
We introduce and analyze multilevel Monte Carlo algorithms for the computation of $\mathbb {E}f(Y)$, where $Y=(Y_t)_{t\in[0,1]}$ is the solution of a multidimensional L\'{e}vy-driven stochastic differential equation and $f$ is a real-valued…
We construct a coupling between the random walk composed of L\'evy area increments from a $d$-dimensional Brownian motion and a random walk composed of quadratic polynomials of Gaussian random variables. This coupling construction is used…
An American option grants the holder the right to select the time at which to exercise the option, so pricing an American option entails solving an optimal stopping problem. Difficulties in applying standard numerical methods to complex…
We consider the path approximation of Bessel processes and develop a new and efficient algorithm. This study is based on a recent work by the authors, on the path approximation of the Brownian motion, and on the construction of specific own…