CVT_WEIGHT - Try to Prescribe CVT Cell Sizes

CVT_WEIGHT was yet another other 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.

Files you may copy include:

• cvt_weight.f90, the source code.
• cvt_weight_prb.f90, a sample problem.
• cvt_weight_prb.out, sample problem output.

The list of routines includes:

• 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.

Return to the FORTRAN software page.