Capacity-Constrained Delaunay Triangulation for point distributions
on centroidal voronoi tessellation - energy smoothness and fast computation .