Related papers: Implementable coupling of L\'evy process and Brown…
Mixed level orthogonal arrays are basic structures in experimental design. We develop three algorithms that compute Rao and Gilbert-Varshamov type bounds for mixed level orthogonal arrays. The computational complexity of the terms involved…
We develop a novel approximate simulation algorithm for the joint law of the position, the running supremum and the time of the supremum of a general L\'evy process at an arbitrary finite time. We identify the law of the error in simple…
One of the outstanding problems in the numerical discretization of the Feynman-Kac formula calls for the design of arbitrary-order short-time approximations that are constructed in a stable way, yet only require knowledge of the potential…
We present a hybrid Brownian dynamics / Monte Carlo algorithm for simulating solutions of highly entangled semiflexible polymers or filaments. The algorithm combines a Brownian dynamics time-stepping approach with an efficient scheme for…
Given an n-dimensional stochastic process X driven by P-Brownian motions and Poisson random measures, we seek the probability measure Q, with minimal relative entropy to P, such that the Q-expectations of some terminal and running costs are…
The first passage time process of a L\'evy subordinator with heavy-tailed L\'evy measure has long-range dependent paths. The random fluctuations that appear under two natural schemes of summation and time scaling of such stochastic…
For a broad class of planar Markov processes, viz. L\'evy processes satisfying certain conditions (valid \textit{eg} in the case of Brownian motion and L\'evy flights), we establish an exact, universal formula describing the shape of the…
Rough path analysis can be developed using the concept of controlled paths, and with respect to a topology in which L\'evy's area plays a role. For vectors of irregular paths we investigate the relationship between the property of being…
We analyze a specific class of random systems that are driven by a symmetric L\'{e}vy stable noise. In view of the L\'{e}vy noise sensitivity to the confining "potential landscape" where jumps take place (in other words, to environmental…
We demonstrate effectiveness of the first-order algorithm from [Milstein, Tretyakov. Theory Prob. Appl. 47 (2002), 53-68] in application to barrier option pricing. The algorithm uses the weak Euler approximation far from barriers and a…
In this paper, we address rare-event simulation for heavy-tailed L\'evy processes with infinite activities. The presence of infinite activities poses a critical challenge, making it impractical to simulate or store the precise sample path…
We investigate the properties of the path integral over reparametrizations (= the boundary value of the Liouville field in open string theory). Discretizing the path integral, we apply the Metropolis-Hastings algorithm to numerical…
The joint distribution of maximum increase and decrease for Brownian motion up to an independent exponential time is computed. This is achieved by decomposing the Brownian path at the hitting times of the infimum and the supremum before the…
The signature is a collection of iterated integrals describing the "shape" of a path. It appears naturally in the Taylor expansions of controlled differential equations and, as a consequence, is arguably the central object within rough path…
The approach to the theory of a relativistic random process is considered by the path integral method as Brownian motion taking into account the boundedness of speed. An attempt was made to build a relativistic analogue of the Wiener…
Assume a L\'evy process $X$ on the time interval $[0,1]$ that is an $L_2$-martingale and let $Y$ be either its stochastic exponential or $X$ itself. We consider Riemann-approximations of certain stochastic integrals driven by $Y$ and relate…
This article provides an overview of recent work on descriptions and properties of the convex minorant of random walks and L\'evy processes which summarize and extend the literature on these subjects. The results surveyed include point…
Estimation methods for the L\'{e}vy density of a L\'{e}vy process are developed under mild qualitative assumptions. A classical model selection approach made up of two steps is studied. The first step consists in the selection of a good…
The linear fractional stable motion generalizes two prominent classes of stochastic processes, namely stable L\'evy processes, and fractional Brownian motion. For this reason it may be regarded as a basic building block for continuous time…
We describe a novel algorithm for rounding packing integer programs based on multidimensional Brownian motion in $\mathbb{R}^n$. Starting from an optimal fractional feasible solution $\bar{x}$, the procedure converges in polynomial time to…