cvt_1d_lumping_test

cvt_1d_lumping_test, a MATLAB code which calls cvt_1d_lumping(), which carries out a lumped version of Lloyd's iterative algorithm for the computation of a centroidal Voronoi Tessellation (CVT) in the interval [-1,+1], for a general density function rho(x).

Licensing:

The computer code and data files described and made available on this web page are distributed under the MIT license

Related Data and Programs:

cvt_1d_lumping, a MATLAB code which allows the user to carry out a lumped version of Lloyd's iterative algorithm for a centroidal Voronoi Tessellation (CVT() in the interval [-1,+1], and is applied to investigate a relationship between CVT's and the Chebyshev Zero nodes.

Source Code:

cvt_1d_lumping_test.m, calls all the tests.
cvt_1d_lumping_test.sh, runs all the tests.
cvt_1d_lumping_test.txt, the output file.

density.png, a plot of densities rho(x) and mu(x), truncating values above 10.

The code was run with the Chebyshev Zero Node density function, N = 3, 5, 9, 17, 33, 65, 129 and 255, with the generator values initialized with equaly spaced values (U).

The output files:

output_003_u.txt, using 250 steps and 100,000 samples and uniform initialization.
output_005_u.txt, using 250 steps and 100,000 samples and uniform initialization.
output_009_u.txt, using 250 steps and 100,000 samples and uniform initialization.
output_017_u.txt, using 250 steps and 100,000 samples and uniform initialization.
output_033_u.txt, using 750 steps and 100,000 samples and uniform initialization.
output_065_u.txt, using 2000 steps and 100,000 samples and uniform initialization.
output_129_u.txt, using 6,000 steps and 200,000 samples and uniform initialization.
output_257_u.txt , using 15,000 steps and 200,000 samples and uniform initialization.

Data Files:

Files containing the coordinates of the Chebyshev Zero nodes, and the CVT generators initialized by Chebyshev Zero nodes, by random values, or by uniform values:

The plots of energy as a function of iteration step:

The plots of the norm of the change in position of the generators as a function of iteration step:

The plots of the position of the generators as a function of iteration step:

The code was run with the Chebyshev Zero Node density function, N = 3, 5, 9, 17, 33, 65, 129 and 255, with the generator values initialized with the Chebyshev Zero Nodes (C). It was expected that, in these runs, the CVT generators would stay near their initial values, so that the motion plots, in particular, would be near vertical.

The output files:

output_003_c.txt, using 250 steps and 100,000 samples and uniform initialization.
output_005_c.txt, using 250 steps and 100,000 samples and uniform initialization.
output_009_c.txt, using 250 steps and 100,000 samples and uniform initialization.
output_017_c.txt, using 250 steps and 100,000 samples and uniform initialization.
output_033_c.txt, using 750 steps and 100,000 samples and uniform initialization.
output_065_c.txt, using 2000 steps and 100,000 samples and uniform initialization.
output_129_c.txt, using 6,000 steps and 200,000 samples and uniform initialization.
output_257_c.txt , using 15,000 steps and 200,000 samples and uniform initialization.

The plots of energy as a function of iteration step:

The plots of the norm of the change in position of the generators as a function of iteration step:

The plots of the position of the generators as a function of iteration step:

Last revised on 14 December 2018.