The Warp::squareToUniformDiskPdf(const Point2f& pt) function passes points all over the uniform square and not just within the disk. Hence you are returning a pdf value of 1/2pi for all points (and hence your integrated density is uniform over the test region). You will have to do a rejection test to check if the point is within the disk (basically x^2 + y^2 <= 1) and return 1/2pi for only those points and zero everywhere else.
On Wednesday, 7 June 2017 04:58:43 UTC-4, 卡卡西 wrote:

