cvtTorus

Centroidal Voronoi tessellation on a 2D torus

Developer documentation: this item describes internal implementation details.

---

Discussion

[x,y] = cvtTorus(N, Lx, Ly) [x,y] = cvtTorus(N, Lx, Ly, nIter)

N – number of generators Lx,Ly – domain size (periodic in both x and y) nIter – number of Lloyd iterations (default: 20)

x,y – N×1 vectors of final seed positions in [0,Lx)×[0,Ly)

Example: [x,y] = cvtTorus(100, 4, 2, 30); scatter(x,y,20,’filled’); axis equal; xlim([0 4]); ylim([0 2]);