We consider the approximate Euler scheme for Lévy-driven stochastic differential equations. We study the rate of convergence in law of the paths. We show that when approximating the small jumps by Gaussian variables, the convergence is much faster than when simply neglecting them. For example, when the Lévy measure of the driving process behaves like |z|−1−αdz near 0, for some α ∈ (1,2), we obtain an error of order 1/√n with a computational cost of order nα. For a similar error when neglecting the...

We consider the approximate Euler scheme for Lévy-driven stochastic differential equations. We study the rate of convergence in law of the paths. We show that when approximating the small jumps by Gaussian variables, the convergence is much faster than when simply neglecting them. For example, when the Lévy measure of the driving process behaves like |z|−1−αdz near 0, for some α∈ (1,2), we obtain an error of order 1/√n with a computational cost of order nα. For a similar error when neglecting the...

The numerical approximation of parametric partial differential equations is a computational challenge, in particular when the number of involved parameter is large. This paper considers a model class of second order, linear, parametric, elliptic PDEs on a bounded domain D with diffusion coefficients depending on the parameters in an affine manner. For such models, it was shown in [9, 10] that under very weak assumptions on the diffusion coefficients, the entire family of solutions to such equations...

Let $AV\to V^{\text{'}}$ be a strongly elliptic operator on a $d$-dimensional manifold $D$ (polyhedra or boundaries of polyhedra are also allowed). An operator equation $Au=f$ with stochastic data $f$ is considered. The goal of the computation is the mean field and higher moments ${ℳ}^{1}u\in V$, ${ℳ}^{2}u\in V\otimes V$, $...$, ${ℳ}^{k}u\in V\otimes \cdots \otimes V$ of the solution. We discretize the mean field problem using a FEM with hierarchical basis and $N$ degrees of freedom. We present a Monte-Carlo algorithm and a deterministic algorithm for the approximation of the moment ${ℳ}^{k}u$ for $k\ge 1$. The key tool...

We use the scale of Besov spaces $B_{\tau ,\tau }^{\alpha }\left(\right)$, 1/τ = α/d + 1/p, α > 0, p fixed, to study the spatial regularity of solutions of linear parabolic stochastic partial differential equations on bounded Lipschitz domains ⊂ ℝ. The Besov smoothness determines the order of convergence that can be achieved by nonlinear approximation schemes. The proofs are based on a combination of weighted Sobolev estimates and characterizations of Besov spaces by wavelet expansions.

The aim of this paper is to take an in-depth look at the long time behaviour of some continuous time markovian dynamical systems and at its numerical analysis. We first propose a short overview of the main ergodicity properties of time continuous homogeneous Markov processes (stability, positive recurrence). The basic tool is a Lyapunov function. Then, we investigate if these properties still hold for the time discretization of these processes, either with constant or decreasing step (ODE method...

The aim of this paper is to take an in-depth look at the long time behaviour of some continuous time Markovian dynamical systems and at its numerical analysis. We first propose a short overview of the main ergodicity properties of time continuous homogeneous Markov processes (stability, positive recurrence). The basic tool is a Lyapunov function. Then, we investigate if these properties still hold for the time discretization of these processes, either with constant or decreasing step (ODE...