CVTITERATE Update germs of a CVT using random points with given density.
G2 = cvtIterate(G, FPTR, FARGS, N)
G: inital germs
FPTR: pointer to a function which accept a scalar M and return M random
points with a given distribution
FARGS: arguments to be given to the FPTR function (can be empty)
N: number of random points to generate
Example
P = randPointDiscUnif(50);
P2 = cvtIterate(P, @randPointDiscUnif, [], 1000);
P3 = cvtIterate(P2, @randPointDiscUnif, [], 1000);
See also
Rewritten from programs found in
http://people.scs.fsu.edu/~burkardt/m_src/cvt/cvt.html
Reference:
Qiang Du, Vance Faber, and Max Gunzburger,
Centroidal Voronoi Tessellations: Applications and Algorithms,
SIAM Review, Volume 41, 1999, pages 637-676.
Package: matgeom