March 9 2003 2:02:08.845 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 = 1252809912 Enter command:> # romero_05.inp # Command echo has been enabled. Enter command:> generator_read = romero_05_lhs.txt FILE_ROW_COUNT: Number of records: 13 Number of data records: 9 Number of comment records: 4 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.9406000 0.5214000 2 0.1817000 0.9738000 3 0.1006000 0.4012000 4 0.3664000 0.7240000 5 0.2801000 0.2913000 6 0.7888000 0.8222000 7 0.7377000 0.5689000 8 0.5904000 4.0199999E-02 9 0.5022000 0.2176000 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: 9 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 1252809912 Enter command:> iterate STEP L2 Change 1 0.226420 2 0.104483 3 0.060029 4 0.046580 5 0.042343 6 0.035203 7 0.026918 8 0.025863 9 0.019138 10 0.015712 11 0.011710 12 0.012988 13 0.011507 14 0.009647 15 0.008798 16 0.008927 17 0.007453 18 0.007560 19 0.008889 20 0.008386 21 0.008969 22 0.006796 23 0.005678 24 0.005561 25 0.008541 26 0.007389 27 0.007695 28 0.008095 29 0.008717 30 0.008107 31 0.008020 32 0.008466 33 0.007710 34 0.008212 35 0.005351 36 0.006861 37 0.006469 38 0.007938 39 0.008246 40 0.005894 41 0.007534 42 0.005753 43 0.010265 44 0.005992 45 0.008901 46 0.004858 47 0.004756 48 0.006881 49 0.007848 50 0.005398 51 0.007131 52 0.007792 53 0.005968 54 0.004913 55 0.005870 56 0.008139 57 0.007059 58 0.006931 59 0.005313 60 0.007998 61 0.005244 62 0.005538 63 0.009129 64 0.005793 65 0.004768 66 0.006261 67 0.008027 68 0.007183 69 0.006816 70 0.008248 71 0.006637 72 0.005695 73 0.007643 74 0.007764 75 0.005377 76 0.006612 77 0.006965 78 0.008031 79 0.007457 80 0.007994 81 0.006368 82 0.006137 83 0.006865 84 0.006265 85 0.007093 86 0.005340 87 0.006096 88 0.006504 89 0.006490 90 0.005454 91 0.004539 92 0.006828 93 0.004853 94 0.007295 95 0.008212 96 0.005647 97 0.006308 98 0.006752 99 0.007451 100 0.007619 Enter command:> generator_print Current cell generators: 1 0.8362098 0.5031891 2 0.1654213 0.8385784 3 0.1694011 0.5017648 4 0.4971613 0.8322296 5 0.1691584 0.1675971 6 0.8311545 0.8333253 7 0.5039968 0.5002521 8 0.8296767 0.1690338 9 0.4986358 0.1671344 Enter command:> maxit = 5 Enter command:> ns_cvt = 20000 Enter command:> iterate STEP L2 Change 1 0.005306 2 0.002811 3 0.003059 4 0.003112 5 0.003725 Enter command:> maxit = 5 Enter command:> ns_cvt = 80000 Enter command:> iterate STEP L2 Change 1 0.002978 2 0.001617 3 0.001438 4 0.001524 5 0.001806 Enter command:> maxit = 5 Enter command:> ns_cvt = 320000 Enter command:> iterate STEP L2 Change 1 0.001409 2 0.000953 3 0.000653 4 0.000777 5 0.000803 Enter command:> generator_write = romero_05_cvt.txt Enter command:> quit QUIT requested. CVT_MAIN Normal end of execution. March 9 2003 2:02:26.469 PM