March 10 2003 3:49:17.369 PM CVT_MAIN A sample problem for the probabilistic Centroidal Voronoi Tessellation algorithm. Given a region in ND, the problem is to determine GENERATORS, a set of points which define a division of the region into Voronoi cells, which are also CENTROIDS of the Voronoi cells. RANDOM_INITIALIZE Initialize RANDOM_NUMBER with arbitrary SEED = 33614256 Enter command:> # romero_10.inp # Command echo has been enabled. Enter command:> generator_read = romero_10_lhs.txt FILE_ROW_COUNT: Number of records: 16 Number of data records: 13 Number of comment records: 3 POINTS_READ: Read coordinate data from file. Enter command:> sample_random The region will be sampled using a random sequence. Enter command:> generator_print Current cell generators: 1 0.1502000 0.1630000 2 2.4000001E-03 0.9493000 3 0.6124000 0.8930000 4 0.9661000 0.8008000 5 0.1698000 0.3897000 6 0.3706000 0.2848000 7 0.4325000 0.5204000 8 0.8882000 3.3300001E-02 9 0.3022000 0.3788000 10 0.8002000 0.5971000 11 0.7291000 0.1469000 12 0.5073000 0.7029000 13 0.6278000 0.6411000 Enter command:> maxit = 100 Enter command:> ns_cvt = 5000 Enter command:> print Geometry parameters: ------------------- The spatial dimension is NDIM = 2 The minimum corner of the bounding box is: 0.00000 0.00000 The maximum corner of the bounding box is: 1.00000 1.00000 A simple user routine determines the region. CVT Algorithm parameters: ------------------------- The number of Voronoi cells to generate: 13 Number of iterations to determine CVT: 100 Number of sampling points per Voronoi cell: 5000 The initial generator positions will be assigned using a vector Halton sequence. The region will be sampled with a random sequence. Moment parameters: ------------------ Number of sampling points per Voronoi cell: 5000 The volume of the region is given. It is specified as REGION_VOLUME = 1.00000 Nearest Neighbor Search parameters: ----------------------------------- The nearest neighbor search is not speeded up. The nearest neighbor search is done by exhaustion. Miscellaneous parameters: ------------------------ Generator and moment output files will NOT be written. The random number seed is 33614256 Enter command:> iterate STEP L2 Change 1 0.283982 2 0.086320 3 0.044420 4 0.030797 5 0.023592 6 0.018756 7 0.015475 8 0.011147 9 0.011864 10 0.007736 11 0.009643 12 0.008207 13 0.010332 14 0.007211 15 0.010063 16 0.009431 17 0.011213 18 0.010789 19 0.012922 20 0.013479 21 0.015733 22 0.015617 23 0.013324 24 0.012815 25 0.011848 26 0.012071 27 0.010403 28 0.010506 29 0.011029 30 0.010881 31 0.009336 32 0.008144 33 0.008318 34 0.008437 35 0.008329 36 0.008746 37 0.009650 38 0.007795 39 0.008046 40 0.007450 41 0.005713 42 0.005513 43 0.005793 44 0.006988 45 0.006584 46 0.006883 47 0.007070 48 0.006435 49 0.006343 50 0.006743 51 0.005446 52 0.006183 53 0.005682 54 0.007723 55 0.006528 56 0.007037 57 0.005403 58 0.006616 59 0.007408 60 0.008041 61 0.007276 62 0.007706 63 0.005777 64 0.005794 65 0.006037 66 0.005964 67 0.006596 68 0.005894 69 0.008823 70 0.006530 71 0.005648 72 0.007811 73 0.007060 74 0.007277 75 0.009470 76 0.006965 77 0.007155 78 0.006229 79 0.005392 80 0.006675 81 0.006692 82 0.005739 83 0.007229 84 0.007689 85 0.004773 86 0.006246 87 0.006335 88 0.007404 89 0.006956 90 0.009027 91 0.006734 92 0.007232 93 0.006385 94 0.005789 95 0.006777 96 0.008297 97 0.007065 98 0.006936 99 0.008282 100 0.007695 Enter command:> generator_print Current cell generators: 1 0.1649961 0.1193649 2 0.1293808 0.8776249 3 0.6320948 0.8873075 4 0.8803160 0.8251315 5 0.1473928 0.6206595 6 0.4774870 0.1316629 7 0.3926906 0.4557660 8 0.8455980 0.1413418 9 0.1341238 0.3553221 10 0.8857743 0.4855550 11 0.6533870 0.3342225 12 0.3836165 0.8226193 13 0.6252835 0.6226684 Enter command:> maxit = 5 Enter command:> ns_cvt = 20000 Enter command:> iterate STEP L2 Change 1 0.004449 2 0.002880 3 0.003095 4 0.002988 5 0.003230 Enter command:> maxit = 5 Enter command:> ns_cvt = 80000 Enter command:> iterate STEP L2 Change 1 0.002482 2 0.001617 3 0.001680 4 0.001772 5 0.001827 Enter command:> maxit = 5 Enter command:> ns_cvt = 320000 Enter command:> iterate STEP L2 Change 1 0.001022 2 0.001170 3 0.000935 4 0.000661 5 0.001341 Enter command:> generator_write = romero_10_cvt.txt Enter command:> quit QUIT requested. CVT_MAIN Normal end of execution. March 10 2003 3:49:40.904 PM