Hostname: page-component-8448b6f56d-sxzjt Total loading time: 0 Render date: 2024-04-24T03:54:42.549Z Has data issue: false hasContentIssue false

Truncated Newton-Based Multigrid Algorithm for Centroidal Voronoi Diagram Calculation

Published online by Cambridge University Press:  28 May 2015

Zichao Di*
Affiliation:
Department of Mathematical Sciences, George Mason University, Fairfax, VA 22030, USA
Maria Emelianenko*
Affiliation:
Department of Mathematical Sciences, George Mason University, Fairfax, VA 22030, USA
Stephen Nash*
Affiliation:
Systems Engineering and Operations Research Department, George Mason University, Fairfax, VA 22030, USA
*
Corresponding author.Email address:zdi@gmu.edu
Corresponding author.Email address:memelian@gmu.edu
Corresponding author.Email address:snash@gmu.edu
Get access

Abstract

In a variety of modern applications there arises a need to tessellate the domain into representative regions, called Voronoi cells. A particular type of such tessellations, called centroidal Voronoi tessellations or CVTs, are in big demand due to their optimality properties important for many applications. The availability of fast and reliable algorithms for their construction is crucial for their successful use in practical settings. This paper introduces a new multigrid algorithm for constructing CVTs that is based on the MG/Opt algorithm that was originally designed to solve large nonlinear optimization problems. Uniform convergence of the new method and its speedup comparing to existing techniques are demonstrated for linear and nonlinear densities for several 1d and 2d problems, and O(k) complexity estimation is provided for a problem with k generators.

Type
Research Article
Copyright
Copyright © Global Science Press Limited 2012

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Aurenhammer, F., Voronoi diagrams. a survey of a fundamental geometric data structure, ACM Computing Surveys, 23 (1990), pp. 345405.CrossRefGoogle Scholar
[2]Cortes, J., Martinez, S., Karata, T., and Bullo, F., Coverage control for mobile sensing networks, IEEE Transactions on Robotics and Automation, 20 (2004), pp. 243255.CrossRefGoogle Scholar
[3]Du, Q. and Emelianenko, M., Acceleration of algorithms for the computation of centroidal Voronoi tessellations, Numerical Linear Algebra and its Applications, 13 (2006), pp. 173192.CrossRefGoogle Scholar
[4]Du, Q., Uniform convergence of a multilevel energy-based quantization scheme, in Lecture Notes in Computer Science and Engineering, Widlund, O. B. and Keyes, D. E., eds., vol. 55, Springer, Berlin, 2007, pp. 533541.Google Scholar
[5]Du, Q.Uniform convergence of a nonlinear energy-based multilevel quantization scheme via centroidal Voronoi tessellations, SIAM Journal on Numerical Analysis, 46 (2008), pp. 14831502.CrossRefGoogle Scholar
[6]Du, Q., Emelianenko, M., and Ju, L., Convergence properties of the Lloyd algorithm for computing the centroidal Voronoi tessellations, SIAM Journal on Numerical Analysis, 44 (2006), pp. 102119.CrossRefGoogle Scholar
[7]Du, Q., Faber, V., and Gunzburger, M., Centroidal Voronoi tessellations: applications and algorithms, SIAM Review, 41 (1999), pp. 637676.CrossRefGoogle Scholar
[8]Du, Q. and Gunzburger, M., Grid generation and optimization based on centroidal Voronoi tessellations, Applied Mathematics and Computation, 133 (2002), pp. 591607.CrossRefGoogle Scholar
[9]Du, Q., Gunzburger, M., and Ju, L., Meshfree probabilistic determination of points, support spheres, and connectivities for meshless computing, Computational Methods in Applied Me chanics and Engineering, 191 (2002), pp. 13491366.CrossRefGoogle Scholar
[10]Du, Q., Constrained centroidal Voronoi tessellations on general surfaces, SIAM Journal on Sci entific Computing, 24 (2003), pp. 14881506.CrossRefGoogle Scholar
[11]Du, Q., Gunzburger, M., and Ju, L., Advances in studies and applications of centroidal voronoi tessellations, NMTMA, 3 (2010), pp. 119142.CrossRefGoogle Scholar
[12]Du, Q. and Wang, D., Tetrahedral mesh generation and optimization based on centroidal Voronoi tessellations, International Journal for Numerical Methods in Engineering, 56 (2002), pp. 13551373.CrossRefGoogle Scholar
[13]Du, Q.Recent progress in robust and quality Delaunay mesh generation, Journal of Computa tional and Applied Mathematics, 195 (2006), pp. 823.CrossRefGoogle Scholar
[14]Dwyer, R., Higher-dimensional Voronoi diagrams in linear expected time, Discrete and Computational Geometry, 6 (1991), pp. 343367.CrossRefGoogle Scholar
[15]Emelianenko, M., Fast multilevel CVT-based adaptive data visualization algorithm, Numerical Mathematics: Theory, Methods and Applications, 3 (2010), pp. 195211.Google Scholar
[16]Emelianenko, M., Ju, L., and Rand, A., Weak global convergence of the Lloyd method for computing centroidal Voronoi tessellations in Rd, SIAM Journal on Numerical Analysis, 46 (2008), pp. 14231441.CrossRefGoogle Scholar
[17]Fortune, S., A sweepline algorithm for Voronoi diagrams, Algorithmica, 2 (1987), pp. 153174.CrossRefGoogle Scholar
[18]Gersho, A., Asymptotically optimal block quantization, IEEE Transactions on Information Theory, 25 (1979), pp. 373380.CrossRefGoogle Scholar
[19]Gray, R. and Neuhoff, D., Quantization, IEEE Transactions on Information Theory, 44 (1998), pp. 23252383.CrossRefGoogle Scholar
[20]Griva, I., Nash, S., and Sofer, A., Linear and Nonlinear Optimization, SIAM, Philadelphia, 2008.Google Scholar
[21]Hiller, S., Hellwig, H., and Deussen, O., Beyond stippling methods for distributing objects on the plane, Computer Graphics Forum, 22 (2003), pp. 515522.CrossRefGoogle Scholar
[22]Ju, L., Du, Q., and Gunzburger, M., Probabilistic methods for centroidal Voronoi tessellations and their parallel implementations, Parallel Computing, 28 (2002), pp. 14771500.CrossRefGoogle Scholar
[23]Koren, Y. and Yavneh, I., Adaptive multiscale redistributionfor vector quantization, SIAM Journal on Scientific Computing, 27 (2006), pp. 15731593.CrossRefGoogle Scholar
[24]Koren, Y., Yavneh, I., and Spira, A., A multigrid approach to the 1-D quantization problem, tech. report, 2003.Google Scholar
[25]Lévy, B. and Liu, Y, Lp centroidal Voronoi tessellation and its applications, in ACM Transactions on Graphics, SIGGRAPH conference proceedings, vol. 24, 2010. Article 119.Google Scholar
[26]Lewis, R. and Nash, S., Model problems for the multigrid optimization of systems governed by differential equations, SIAM Journal on Scientific Computing, 26 (2005), pp. 18111837.CrossRefGoogle Scholar
[27]Linde, Y., Buzo, A., and Gray, R., An algorithm for vector quantizer design, IEEE Transactions on Communications, 28 (1980), pp. 8495.CrossRefGoogle Scholar
[28]Lloyd, S., Least square quantization in PCM, IEEE Transactions on Information Theory, 28 (1982), pp. 129137.CrossRefGoogle Scholar
[29]Mccormick, S., Multilevel Projection Methods for Partial Differential Equations, Society for Industrial and Applied Mathematics, 1992.Google Scholar
[30]Mendes, A. and Themido, I., Multi-outlet retail site location assessment, International Transactions in Operational Research, 11 (2004), pp. 118.CrossRefGoogle Scholar
[32]Nash, S., A multigrid approach to discretized optimization problems, Journal of Computational and Applied Mathematics, 14 (2000), pp. 99116.Google Scholar
[33]Nash, S., Convergence and descent properties for a class of multilevel optimization algorithms, tech. report, Systems Engineering and Operations Research Department, George Mason University, Fairfax, VA, 2010.Google Scholar
[34]Nash, S. and Nocedal, J., A numerical study of the limited memory BFGS method and the truncated-Newton method for large scale optimization, 1 (1991), pp. 358372.Google Scholar
[35]Okabe, A., Boots, B., and Sugihara, K., Spatial Tessellations; Concepts and Applications of Voronoi Diagrams, Wiley, Chichester, 1992.Google Scholar
[36]Rong, G., Liu, Y., Wang, W., Yin, X., Gu, X., and Guo, X., GPU-assisted computation of centroidal Voronoi tessellation, Visualization and Computer Graphics, IEEE Transactions on, 17 (2011), pp. 345–356.Google ScholarPubMed
[37]Sabin, J. and Gray, R., Global convergence and empirical consistency of the generalized Lloyd algorithm, IEEE Transactions on Information Theory, IT–32 (1986), pp. 148–155.Google Scholar
[38]Trushkin, A., On the design of an optimal quantizer, IEEE Transactions on Information Theory, 39 (1993), pp. 1180–1194.CrossRefGoogle Scholar
[39]Trottenberg, U., Oosterlee, C.W., and Schuller, A., Multigrid, (2001), pp. 61–62.Google Scholar
[40]Valette, S. and Chassery, J., Approximated centroidal Voronoi diagrams for uniform polygonal mesh coarsening, Computer Graphics Forum, 23 (2004), pp. 381–390.CrossRefGoogle Scholar
[41]Vasconcelos, C. N., , A., Carvalho, P. C., and Gattass, M., Lloyd’s algorithm on GPU, in ISVC ‘08: Proceedings of the 4th International Symposium on Advances in Visual Computing, Berlin, Heidelberg, 2008, Springer-Verlag, pp. 953–964.Google Scholar
[42]Wager, C., Coull, B., and Lange, N., Modelling spatial intensity for replicated inhomogeneous point patterns in brain imaging, Journal of Royal Statistical Society B.Google Scholar
[43]Liu, Y., Wang, W., Lévy, B., Sun, F., Yan, D., Lu, L., and Yang, Ch., On centroidal voronoi tessellation-energy smoothness and fast computation, ACM Trans. Graph. 28, 4, Article 101, 2009.CrossRefGoogle Scholar