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