An Analysis of Approximation Algorithms for Iterated Stochastic Integrals and a Julia and MATLAB Simulation Toolbox

Abstract

For the approximation and simulation of twofold iterated stochastic integrals and the corresponding Lévy areas w.r.t. a multi-dimensional Wiener process, we review four algorithms based on a Fourier series approach. Especially, the very efficient algorithm due to Wiktorsson and a newly proposed algorithm due to Mrongowius and Rössler are considered. To put recent advances into context, we analyse the four Fourier-based algorithms in a unified framework to highlight differences and similarities in their derivation. A comparison of theoretical properties is complemented by a numerical simulation that reveals the order of convergence for each algorithm. Further, concrete instructions for the choice of the optimal algorithm and parameters for the simulation of solutions for stochastic (partial) differential equations are given. Additionally, we provide advice for an efficient implementation of the considered algorithms and incorporated these insights into an open source toolbox that is freely available for both Julia and MATLAB programming languages. The performance of this toolbox is analysed by comparing it to some existing implementations, where we observe a significant speed-up.

References in corpus (6)

• Efficiently and easily integrating differential equations with JiTCODE, JiTCDDE, and JiTCSDE
• Development and Application of the Fourier Method to the Mean-Square Approximation of Iterated Ito and Stratonovich Stochastic Integrals
• Stochastic Runge-Kutta Software Package for Stochastic Differential Equations
• Comparative Analysis of the Efficiency of Application of Legendre Polynomials and Trigonometric Functions to the Numerical Integration of Ito Stochastic Differential Equations
• Brownian bridge expansions for Lévy area approximations and particular values of the Riemann zeta function
• Issues in the software implementation of stochastic numerical Runge-Kutta

Cited by in corpus (1)

• Generative Modelling of Lévy Area for High Order SDE Simulation