Implemented random point generation within cuboids, spheres and
cylinders. The method always finds a point in these shapes unlike the
generic one and is sometimes faster (depends on how the bounding box is
filled). This should make MonteCarloAbsorption more robust.

Re #24636