A review on numerical schemes for solving a linear stochastic oscillator

In recent years several numerical methods to solve a linear stochastic oscillator with one additive noise have been proposed. The usual aim of these approaches was to preserve different long time properties of the oscillator solution. In this work we collect these properties, namely, symplecticity, linear growth of its second moment and asymptotic oscillation around zero. We show that these features can be studied in terms of the coefficients of the matrices that appear in the linear recurrence obtained when the schemes are applied to the oscillator. We use this study to compare the numerical schemes as well as to propose new schemes improving some properties of classical methods.