CVT_MAX - Try to Prescribe CVT Cell Sizes

CVT_MAX was an 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. Each node was given a "capacity", corresponding to its desired size. Then all the sample points were generated, and "marked" by their nearest node. Then each node took as many of its nearest sample points as possible, up to its capacity. The remaining sample points were reapportioned to the next nearest nodes, and so on until all were allotted.

At the moment, this approach has not produced satisfactory results.

Files you may copy include:

cvt_max.f90, the source code.

cvt_max_prb.f90, a sample problem.

cvt_max_prb.out, sample problem output.

The list of routines includes:

CELL_VOLUME_COMPUTATION estimates the cell volumes by sampling.

CVT_MAX_ITERATION takes one step of the CVT Max iteration.

FIND_CLOSEST finds the Voronoi cell generator closest to a point X.

FIND_CLOSEST_AVAILABLE matches points to nearest available generators.

RANDOM_INITIALIZE initializes the FORTRAN 90 random number seed.

REGION_SAMPLER returns a sample point in the physical region.

TIMESTAMP prints the current YMDHMS date as a time stamp.

Return to the FORTRAN software page.

Last revised on 29 October 2001.