added a random point when centroids run out
refactoring
got one r/fold parallel function in, 2.5x speedup for large ?