Set of C++ classes for integrating the Langevin equation using Monte Carlo integrators with different time stepping methods. Both Standard Monte Carlo and more advanced Multilevel Monte Carlo methods are supported. The code can be used to solve several SDEs: Simple test cases, such as propagation in a quadratic and quartic oscillator potential Propagation of particles in a turbulent velocity field in the atmospheric boundary layer For more details see https://bitbucket.org/em459/mlmclangevin The code provided here can be used to reproduce the results in the paper "Multilevel Monte Carlo and Improved Time Stepping Methods in Atmospheric Modelling" by Grigoris Katsiolides, Eike Mueller, Robert Scheichl, Tony Shardlow and David Thomson. (tag 'JCPpaper_v2', SHA da54f9a in bitbucket repository)

{"references": ["preprint version of paper https://arxiv.org/abs/1612.07717"]}