March 9 2003 1:56:24.799 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 = 1249378849 Enter command:> # romero_01.inp # Command echo has been enabled. Enter command:> generator_read = romero_01_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.1250000 0.3842000 2 4.7499999E-02 0.8654000 3 0.3991000 0.4988000 4 0.8701000 0.6192000 5 0.5866000 0.7222000 6 0.4495000 0.2087000 7 0.7302000 0.9465000 8 0.2253000 0.1003000 9 0.9594000 0.2942000 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 1249378849 Enter command:> iterate STEP L2 Change 1 0.203803 2 0.072349 3 0.052477 4 0.038279 5 0.022959 6 0.018035 7 0.012241 8 0.011880 9 0.012476 10 0.010160 11 0.010073 12 0.007474 13 0.007412 14 0.007837 15 0.004586 16 0.007038 17 0.006239 18 0.007501 19 0.008137 20 0.005310 21 0.006911 22 0.006539 23 0.005384 24 0.005966 25 0.006964 26 0.007965 27 0.007854 28 0.007394 29 0.004123 30 0.005764 31 0.006852 32 0.007061 33 0.008258 34 0.008555 35 0.005750 36 0.007409 37 0.008127 38 0.007311 39 0.007309 40 0.008550 41 0.006878 42 0.006359 43 0.008184 44 0.006941 45 0.006617 46 0.006210 47 0.006138 48 0.008243 49 0.006466 50 0.008707 51 0.006029 52 0.009116 53 0.009300 54 0.005372 55 0.005797 56 0.005846 57 0.006081 58 0.005231 59 0.006949 60 0.006384 61 0.005559 62 0.008028 63 0.008221 64 0.007531 65 0.005888 66 0.005622 67 0.006423 68 0.007541 69 0.005264 70 0.006822 71 0.007806 72 0.006065 73 0.006431 74 0.004939 75 0.004989 76 0.007495 77 0.006979 78 0.009232 79 0.008014 80 0.008190 81 0.007441 82 0.006889 83 0.006653 84 0.006217 85 0.005965 86 0.006582 87 0.007624 88 0.007218 89 0.006510 90 0.007218 91 0.007349 92 0.007356 93 0.007408 94 0.007335 95 0.005382 96 0.005790 97 0.006377 98 0.005799 99 0.005614 100 0.009087 Enter command:> generator_print Current cell generators: 1 0.1668702 0.4945796 2 0.1677043 0.8327172 3 0.4953975 0.5091925 4 0.8322953 0.4974210 5 0.5060905 0.8386984 6 0.5046231 0.1727597 7 0.8323433 0.8296763 8 0.1688944 0.1632394 9 0.8332843 0.1633820 Enter command:> maxit = 5 Enter command:> ns_cvt = 20000 Enter command:> iterate STEP L2 Change 1 0.005400 2 0.003649 3 0.004098 4 0.003532 5 0.002631 Enter command:> maxit = 5 Enter command:> ns_cvt = 80000 Enter command:> iterate STEP L2 Change 1 0.002002 2 0.001610 3 0.001052 4 0.001466 5 0.001386 Enter command:> maxit = 5 Enter command:> ns_cvt = 320000 Enter command:> iterate STEP L2 Change 1 0.001019 2 0.001023 3 0.000929 4 0.001157 5 0.001101 Enter command:> generator_write = romero_01_cvt.txt Enter command:> quit QUIT requested. CVT_MAIN Normal end of execution. March 9 2003 1:56:43.192 PM