March 9 2003 2:26:27.982 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 = 1267401348 Enter command:> # romero_12.inp # Command echo has been enabled. Enter command:> generator_read = romero_12_lhs.txt FILE_ROW_COUNT: Number of records: 45 Number of data records: 41 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.1175000 0.4289000 2 0.8223000 0.2364000 3 0.7774000 0.2882000 4 0.1371000 0.1266000 5 0.9718000 0.5882000 6 0.5064000 0.2109000 7 0.7300000 0.4498000 8 5.3900000E-02 0.3191000 9 0.6505000 0.8944000 10 0.3687000 3.8800001E-02 11 0.9772000 0.4707000 12 0.4821000 0.8566000 13 0.8758000 0.6240000 14 4.7600001E-02 0.1826000 15 0.9062000 0.9264000 16 0.3276000 0.2680000 17 0.8435000 0.7763000 18 0.2816000 0.3022000 19 0.9510000 5.4200001E-02 20 0.4406000 7.5499997E-02 21 0.1609000 0.5602000 22 0.3063000 0.4003000 23 0.6718000 6.8000001E-03 24 0.3636000 0.1591000 25 0.1991000 0.5825000 26 0.5710000 0.3444000 27 0.2537000 0.8264000 28 0.5203000 0.1190000 29 0.4074000 0.8486000 30 0.8978000 0.4978000 31 9.5799997E-02 0.9937000 32 0.4372000 0.3804000 33 0.1942000 0.6958000 34 0.7403000 0.9668000 35 0.8013000 0.6343000 36 0.5399000 0.6620000 37 7.9999998E-04 0.9461000 38 0.6181000 0.7268000 39 0.6002000 0.5316000 40 0.2312000 0.7325000 41 0.7031000 0.7933000 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: 41 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 1267401348 Enter command:> iterate STEP L2 Change 1 0.249005 2 0.080649 3 0.046231 4 0.035346 5 0.029963 6 0.026143 7 0.022867 8 0.022239 9 0.019844 10 0.018210 11 0.018414 12 0.016130 13 0.014951 14 0.014340 15 0.014078 16 0.012848 17 0.013737 18 0.014310 19 0.013373 20 0.013981 21 0.013829 22 0.012388 23 0.012129 24 0.011463 25 0.010558 26 0.010068 27 0.009461 28 0.009910 29 0.008851 30 0.008635 31 0.008403 32 0.009001 33 0.008114 34 0.009254 35 0.008257 36 0.008655 37 0.007950 38 0.008636 39 0.009019 40 0.009004 41 0.009127 42 0.009323 43 0.007710 44 0.007809 45 0.007780 46 0.008401 47 0.008531 48 0.008554 49 0.008094 50 0.008347 51 0.008729 52 0.007584 53 0.008193 54 0.008574 55 0.008283 56 0.008155 57 0.009117 58 0.008152 59 0.007915 60 0.007944 61 0.007482 62 0.007786 63 0.008571 64 0.008059 65 0.007712 66 0.008309 67 0.007355 68 0.007419 69 0.007641 70 0.007732 71 0.007310 72 0.007345 73 0.006840 74 0.006027 75 0.007413 76 0.007069 77 0.007356 78 0.006926 79 0.006411 80 0.006733 81 0.006488 82 0.006550 83 0.006923 84 0.006273 85 0.007156 86 0.007384 87 0.007986 88 0.007568 89 0.007289 90 0.006996 91 0.007820 92 0.006330 93 0.007315 94 0.007317 95 0.007984 96 0.006810 97 0.006520 98 0.006208 99 0.006729 100 0.006764 Enter command:> generator_print Current cell generators: 1 7.7579521E-02 0.3959518 2 0.9035776 0.2204787 3 0.7184536 0.2266971 4 0.2086171 0.1067864 5 0.9351041 0.6492948 6 0.5402927 0.2301500 7 0.7831853 0.3668121 8 7.9897434E-02 0.2336455 9 0.6032556 0.9259511 10 0.3405432 6.4101197E-02 11 0.9334112 0.3644068 12 0.5290233 0.7847731 13 0.7848019 0.6480762 14 7.1809940E-02 7.4290290E-02 15 0.9258428 0.9285131 16 0.4171818 0.3628013 17 0.8999201 0.7883816 18 0.2359511 0.2930532 19 0.9310983 7.5923361E-02 20 0.4880656 7.7490754E-02 21 8.2066514E-02 0.5627913 22 0.2264645 0.4539263 23 0.7866991 8.1360318E-02 24 0.3704045 0.2123951 25 0.2451092 0.6324946 26 0.6108940 0.3691836 27 0.2585283 0.9302479 28 0.6367715 7.9727277E-02 29 0.4308691 0.9263958 30 0.9019615 0.5081557 31 7.9377674E-02 0.9170970 32 0.3654315 0.5106044 33 0.1993869 0.7887648 34 0.7687488 0.9252868 35 0.7121904 0.5064530 36 0.4375764 0.6460083 37 6.9782712E-02 0.7301021 38 0.6145105 0.6437070 39 0.5318750 0.5020297 40 0.3556452 0.7784207 41 0.7081181 0.7864454 Enter command:> maxit = 5 Enter command:> ns_cvt = 20000 Enter command:> iterate STEP L2 Change 1 0.005110 2 0.003679 3 0.003201 4 0.003808 5 0.003250 Enter command:> maxit = 5 Enter command:> ns_cvt = 80000 Enter command:> iterate STEP L2 Change 1 0.002202 2 0.001696 3 0.001711 4 0.002022 5 0.001810 Enter command:> maxit = 5 Enter command:> ns_cvt = 320000 Enter command:> iterate STEP L2 Change 1 0.001416 2 0.001019 3 0.001015 4 0.000875 5 0.000963 Enter command:> generator_write = romero_12_cvt.txt Enter command:> quit QUIT requested. CVT_MAIN Normal end of execution. March 9 2003 2:28:50.496 PM