I have created two components for poisson disk sampling. One for a standard 2D sampling and one for sampling in 3D space with spheres instead of discs.
The 3D version is still taking a lot of time which can be shortened by parallelization.
I am looking to release this once I have a geodesic version implemented.