A SIMPLE GEOMETRIC METHOD FOR NAVIGATING THE ENERGY LANDSCAPE OF CENTROIDAL VORONOI TESSELLATIONS

Finding optimal (or low energy) centroidal Voronoi tessellations (CVTs) on a two-dimensional domain is a challenging problem. One must navigate an energy landscape whose desirable critical points have sufficiently small basins of attractions that are inaccessible with Monte Carlo initialized gradient descent methods. We present a simple deterministic method for efficiently navigating the energy landscape in order to access these low energy CVTs. The method has two parameters and is based upon each generator moving away from the closest neighbor by a certain distance. We give a statistical analysis of the performance of this hybrid method comparing with the results of a large number of runs for both Lloyd's method and state-of-the-art quasi-Newton methods. Stochastic alternatives are also considered.