ENH Implement mh.disk in C++ for 2D cases

10x faster (even if this function is generally called for small radii).