Exact Calculation of Stochastic Time Integrals with an Application Towards High Order Monte Carlo of Non-Linear Stochastic Differential Equations

2015 ◽  
Author(s):  
Ahsan Amin
Author(s):  
Adrien Laurent ◽  
Gilles Vilmart

AbstractWe derive a new methodology for the construction of high-order integrators for sampling the invariant measure of ergodic stochastic differential equations with dynamics constrained on a manifold. We obtain the order conditions for sampling the invariant measure for a class of Runge–Kutta methods applied to the constrained overdamped Langevin equation. The analysis is valid for arbitrarily high order and relies on an extension of the exotic aromatic Butcher-series formalism. To illustrate the methodology, a method of order two is introduced, and numerical experiments on the sphere, the torus and the special linear group confirm the theoretical findings.


2020 ◽  
Vol 26 (3) ◽  
pp. 193-203
Author(s):  
Shady Ahmed Nagy ◽  
Mohamed A. El-Beltagy ◽  
Mohamed Wafa

AbstractMonte Carlo (MC) simulation depends on pseudo-random numbers. The generation of these numbers is examined in connection with the Brownian motion. We present the low discrepancy sequence known as Halton sequence that generates different stochastic samples in an equally distributed form. This will increase the convergence and accuracy using the generated different samples in the Multilevel Monte Carlo method (MLMC). We compare algorithms by using a pseudo-random generator and a random generator depending on a Halton sequence. The computational cost for different stochastic differential equations increases in a standard MC technique. It will be highly reduced using a Halton sequence, especially in multiplicative stochastic differential equations.


2019 ◽  
Vol 23 (Suppl. 1) ◽  
pp. 1-12 ◽  
Author(s):  
Burhaneddin Izgi ◽  
Coskun Cetin

We develop Milstein-type versions of semi-implicit split-step methods for numerical solutions of non-linear stochastic differential equations with locally Lipschitz coefficients. Under a one-sided linear growth condition on the drift term, we obtain some moment estimates and discuss convergence properties of these numerical methods. We compare the performance of multiple methods, including the backward Milstein, tamed Milstein, and truncated Milstein procedures on non-linear stochastic differential equations including generalized stochastic Ginzburg-Landau equations. In particular, we discuss their empirical rates of convergence.


Sign in / Sign up

Export Citation Format

Share Document