We place the numerical method of Cruzeiro, Malliavin and Thalmeier for simulation of elliptic diffusions in the context of Riemannian geometry and discuss possible extensions to the hypoelliptic case.