Abstract
We propose and study conformal integrators for linearly damped stochastic Poisson systems. We analyse the qualitative and quantitative properties of these numerical integrators: preservation of dynamics of certain Casimir and Hamiltonian functions, almost sure bounds of the numerical solutions, and strong and weak rates of convergence under appropriate conditions. These theoretical results are illustrated with several numerical experiments on, for example, the linearly damped free rigid body with random inertia tensor or the linearly damped stochastic Lotka-Volterra system.