Global Optimization of Centroidal Voronoi Tessellation with Monte Carlo Approach

published_or_final_versio

Probabilistic and parallel algorithms for centroidal Voronoi tessellations with application to meshless computing and numerical analysis on surfaces

Centroidal Voronoi tessellations (CVT) are Voronoi tessellations of a region such that the generating points of the tessellations are also the centroids of the corresponding Voronoi regions. Such tessellations are of use in very diverse applications, including data compression, clustering analysis, cell biology, territorial behavior of animals, optimal allocation of resources, and grid generation. A detailed review is given in chapter 1. In chapter 2, some probabilistic methods for determining centroidal Voronoi tessellations and their parallel implementation on distributed memory systems are presented. The results of computational experiments performed on a CRAY T3E-600 system are given for each algorithm. These demonstrate the superior sequential and parallel performance of a new algorithm we introduce. Then, new algorithms are presented in chapter 3 for the determination of point sets and associated support regions that can then be used in meshless computing methods. The algorithms are probabilistic in nature so that they are totally meshfree, i.e., they do not require, at any stage, the use of any coarse or fine boundary conforming or superimposed meshes. Computational examples are provided that show, for both uniform and non-uniform point distributions that the algorithms result in high-quality point sets and high-quality support regions. The extensions of centroidal Voronoi tessellations to general spaces and sets are also available. For example, tessellations of surfaces in a Euclidean space may be considered. In chapter 4, a precise definition of such constrained centroidal Voronoi tessellations (CCVT\u27s) is given and a number of their properties are derived, including their characterization as minimizers of a kind of energy. Deterministic and probabilistic algorithms for the construction of CCVT\u27s are presented and some analytical results for one of the algorithms are given. Some computational examples are provided which serve to illustrate the high quality of CCVT point sets. CCVT point sets are also applied to polynomial interpolation and numerical integration on the sphere. Finally, some conclusions are given in chapter 5

A Hierarchical Approach for Regular Centroidal Voronoi Tessellations

International audienceIn this paper we consider Centroidal Voronoi Tessellations (CVTs) and study their regularity. CVTs are geometric structures that enable regular tessellations of geometric objects and are widely used in shape modeling and analysis. While several efficient iterative schemes, with defined local convergence properties, have been proposed to compute CVTs, little attention has been paid to the evaluation of the resulting cell decompositions. In this paper, we propose a regularity criterion that allows us to evaluate and compare CVTs independently of their sizes and of their cell numbers. This criterion allows us to compare CVTs on a common basis. It builds on earlier theoretical work showing that second moments of cells converge to a lower bound when optimising CVTs. In addition to proposing a regularity criterion, this paper also considers computational strategies to determine regular CVTs. We introduce a hierarchical framework that propagates regularity over decomposition levels and hence provides CVTs with provably better regularities than existing methods. We illustrate these principles with a wide range of experiments on synthetic and real models

Decentralized Resource Allocation through Constrained Centroidal Voronoi Tessellations

The advancements in the fields of microelectronics facilitate incorporating team elements like coordination into engineering systems through advanced computing power. Such incorporation is useful since many engineering systems can be characterized as a collection of interacting subsystems each having access to local information, making local decisions, interacting with neighbors, and seeking to optimize local objectives that may well conflict with other subsystems, while also trying to optimize certain global objective. In this dissertation, we take advantage of such technological advancements to explore the problem of resource allocation through different aspects of the decentralized architecture like information structure in a team. Introduced in 1968 as a toy example in the field of team decision theory to demonstrate the significance of information structure within a team, the Witsenhausen counterexample remained unsolved until the analytical person-by-person optimal solution was developed within the past decade. We develop a numerical method to implement the optimal laws and show that our laws coincide with the optimal affine laws. For the region where the optimal laws are non-linear, we show that our laws result in the lowest costs when compared with previously reported costs. Recognizing that, in the framework of team decision theory, the difficulties arising from the non-classical information structure within a team currently limit its applicability in real-world applications, we move on to investigating Centroidal Voronoi Tessellations (CVTs) to solve the resource allocation problem. In one-dimensional spaces, a line communication network is sufficient to obtain CVTs in a decentralized manner, while being scalable to any number of agents in the team. We first solve the static resource allocation problem where the amount of resource is fixed. Using such static allocation solution as an initialization step, we solve the dynamic resource allocation problem in a truly decentralized manner. Furthermore, we allow for flexibility in agents\u27 embedding their local preferences through what we call a civility model. We end the dissertation by revisiting the application of Demand-response in smart grids and demonstrate the developed decentralized dynamic resource allocation method to solve the problem of power allocation in a group of building loads

Pointwise convergence of the Lloyd algorithm in higher dimension

We establish the pointwise convergence of the iterative Lloyd algorithm, also known as $k$-means algorithm, when the quadratic quantization error of the starting grid (with size $N\ge 2$) is lower than the minimal quantization error with respect to the input distribution is lower at level $N-1$. Such a protocol is known as the splitting method and allows for convergence even when the input distribution has an unbounded support. We also show under very light assumption that the resulting limiting grid still has full size $N$. These results are obtained without continuity assumption on the input distribution. A variant of the procedure taking advantage of the asymptotic of the optimal quantizer radius is proposed which always guarantees the boundedness of the iterated grids