As a useful addition to my post, I would like to show a simple way of generating a set of hyperuniformly distributed points with help of the LloydAlgorithm user contributed function available from Wolfram Function Repository. The set of points can be obtained after Lloyd optimization applied to a set of randomly distributed points, as detailed in the following
pts = RandomReal[{0, 10}, {(*number of points *) 100, 2}];
newpts =
ResourceFunction["LloydAlgorithm"][pts, {{0, 10}, {0, 10}},(* nr.
iterations *) 100];
{Show[{VoronoiMesh[pts], Graphics[Point[pts]]}],
Show[{VoronoiMesh[newpts], Graphics[Point[newpts]]}]}
