Volume 17, Issue 3

Fast Fourier transforms are used to develop algorithms for the fast generation of correlated Gaussian random fields on rectangular regions of . The complexities of the algorithms are derived, simulation results and error analysis are presented.

A wavelet-based method for simulation of strictly sub-Gaussian processes with given accuracy and reliability in C ([0, T ]) is proposed.

This work develops novel error expansions with computable leading order terms for the global weak error in the tau-leap discretization of pure jump processes arising in kinetic Monte Carlo models. Accurate computable a posteriori error approximations are the basis for adaptive algorithms, a fundamental tool for numerical simulation of both deterministic and stochastic dynamical systems. These pure jump processes are simulated either by the tau-leap method, or by exact simulation, also referred to as dynamic Monte Carlo, the Gillespie Algorithm or the Stochastic Simulation Slgorithm. Two types of estimates are presented: an a priori estimate for the relative error that gives a comparison between the work for the two methods depending on the propensity regime, and an a posteriori estimate with computable leading order term.

Exact yet simple simulation algorithms are developed for a wide class of Ornstein–Uhlenbeck processes with tempered stable stationary distribution of finite variation with the help of their exact transition probability between consecutive time points. Random elements involved can be divided into independent tempered stable and compound Poisson distributions, each of which can be simulated in the exact sense through acceptance-rejection sampling, respectively, with stable and gamma proposal distributions. We discuss various alternative simulation methods within our algorithms on the basis of acceptance rate in acceptance-rejection sampling for both high- and low-frequency sampling. Numerical results illustrate their advantage relative to the existing approximative simulation method based on infinite shot noise series representation.