On Wed, 5 Nov 2014, Ed Smith-Rowland wrote: > Like the uniform_on_sphere_distribution which is used inside, the > 2-dimensional case uses rejection Could you point out where in the code you are special-casing dimension 2? Somehow I can't see it in the patch. -- Marc Glisse