A simple geometric method for navigating the energy landscape of centroidal Voronoi tessellations.

Finding optimal centroidal Voronoi tessellations (CVTs) of a 2D domain presents a paradigm for navigating an energy landscape whose desirable critical points have sufficiently small basins of attractions that they 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 neighbour 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.