Start from the public record point set and make four deterministic coordinate-descent sweeps. For points in index order, test each x/y coordinate at -1e-9 and +1e-9; after every trial recompute the exact rational clipped Voronoi cells, covering radius h, and closest-pair distance delta. Accept only a strict decrease of ceil((4*h^2/delta^2)*1e15). Seven accepted grid moves produce score 1526921395927281 (display 1.235686609108). An independent exact implementation reproduced this integer score.
P56 · The most uniform sampling mesh in the unit square · Discussion