line_cvt_lloyd, a FORTRAN90 code which carries out Lloyd's iteration for computing a Centroidal Voronoi Tesselation (CVT) of points over the interior of a line segment in 1D.
A constraint has been added to the computation, which forces the first and last points to be fixed to the endpoints of the line segment. This is not actually a requirement for Lloyd's method.
At least for the uniform density case, the exact solution of this problem is known in advance, and is simply a set of equally spaced points. For instance, for N = 4, the solution over the interval [0,1] would be 0, 1/3, 2/3, 1, if the endpoint constraint is imposed. If the endpoint constraint is NOT imposed, then the exact solution is: 1/8, 3/8, 5/8, 7/8.
The computer code and data files described and made available on this web page are distributed under the MIT license
line_cvt_lloyd is available in a C version and a C++ version and a FORTRAN90 version and a MATLAB version.
LINE_FELIPPA_RULE, a FORTRAN90 code which returns the points and weights of a Felippa quadrature rule over the interior of a line segment in 1D.
LINE_GRID, a FORTRAN90 code which computes a grid of points over the interior of a line segment in 1D.