Citation:
Bernal, F. & Acebrón, J. A. (2016). A Comparison of Higher-Order Weak Numerical Schemes for Stopped Stochastic Differential Equations. Communications in Computational Physics, 20(3), pp. 703–732.
Sponsor:
This work was supported by Portuguese national funds through FCT under grants UID/CEC/50021/2013 and PTDC/EIA-CCO/098910/2008. FB also acknowledges FCT funding under grant SFRH/BPD/79986/2011.
We review, implement, and compare numerical integration schemes for spatially bounded diffusions stopped at the boundary which possess a convergence rate of the discretization error with respect to the timestep h higher than . We address specific implementatioWe review, implement, and compare numerical integration schemes for spatially bounded diffusions stopped at the boundary which possess a convergence rate of the discretization error with respect to the timestep h higher than . We address specific implementation issues of the most general-purpose of such schemes. They have been coded into a single Matlab program and compared, according to their accuracy and computational cost, on a wide range of problems in up to ℝ48. The paper is self-contained and the code will be made freely downloadable.[+][-]