Many real world data, particularly in the natural sciences and computer vision, lie on known riemannian manifolds such as spheres, tori or the group of rotation matrices. The predominant approaches to learning a distribution on such a manifold require solving a differential equation in