CVT_WEIGHT
Try to Prescribe CVT Cell Sizes

CVT_WEIGHT is a FORTRAN90 library which implemented an experimental effort to compute a centroidal Voronoi tessellation (CVT) in which there were cells of specified size.

Normally, given arbitrary initial data and a constant density, the CVT would comprise cells of roughly equal area. In this case, a modification was made to the CVT iteration, in an attempt to drive the algorithm towards cells of prescribed sizes. While we were convinced this idea ought to work...it didn't.

Licensing:

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

Languages:

CVT_WEIGHT is available in a FORTRAN90 version.

Source Code:

• cvt_weight.f90, the source code.
• cvt_weight.sh, commands to compile the source code.

Examples and Tests:

• cvt_weight_test.f90, a sample problem.
• cvt_weight_test.sh, commands to compile, link and run the sample problem.
• cvt_weight_test.txt, sample problem output.

List of Routines:

• CELL_VOLUME_COMPUTATION estimates the cell volumes by sampling.
• CVT_DENSE_ITERATION takes one step of the CVT density iteration.
• FIND_CLOSEST finds the Voronoi cell generator closest to a point.
• RANDOM_INITIALIZE initializes the FORTRAN 90 random number seed.
• REGION_SAMPLER returns sample points in the physical region.
• REGION_SAMPLER_WEIGHTED returns weighted sample points
• RVEC_BRACKET searches a sorted array for successive brackets of a value.
• RVEC_EVEN returns N real values, evenly spaced between ALO and AHI.
• SPLINE_CONSTANT_VAL evaluates a piecewise constant spline at a point.
• SPLINE_LINEAR_VAL evaluates a linear spline at a specific point.
• TEST_REGION determines if a point is within the physical region.
• TIMESTAMP prints the current YMDHMS date as a time stamp.

You can go up one level to the FORTRAN90 source codes.