March 9 2003 2:18:03.878 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 = 1262365823 Enter command:> # romero_08.inp # Command echo has been enabled. Enter command:> generator_read = romero_08_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.3444000 0.5558000 2 0.7286000 0.4015000 3 0.6305000 0.7345000 4 0.8347000 9.3199998E-02 5 0.5885000 0.4287000 6 0.9816000 0.6696000 7 0.2283000 0.7024000 8 8.0399998E-02 0.8508000 9 0.3784000 0.7892000 10 0.4524000 0.2814000 11 0.6910000 6.2799998E-02 12 0.4333000 0.9135000 13 0.9482000 0.5751000 14 0.9106000 4.5299999E-02 15 3.9900001E-02 0.1448000 16 0.6669000 0.8749000 17 0.4000000 0.3688000 18 0.3076000 0.7192000 19 0.1619000 0.6264000 20 0.5796000 0.8228000 21 0.5100000 0.5078000 22 0.9012000 0.7756000 23 0.9608000 0.3313000 24 0.8118000 0.9936000 25 0.1102000 7.1000000E-03 26 0.3340000 0.5152000 27 0.2835000 0.9494000 28 0.6383000 0.1683000 29 0.2065000 0.3554000 30 0.8027000 0.4652000 31 0.1296000 0.2208000 32 2.2100000E-02 0.4531000 33 0.4814000 0.8975000 34 0.2673000 0.1970000 35 6.1500002E-02 0.1747000 36 0.1732000 0.9733000 37 0.7699000 0.3141000 38 0.5359000 0.2620000 39 0.7555000 0.1214000 40 0.5504000 0.6058000 41 0.8640000 0.6431000 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 1262365823 Enter command:> iterate STEP L2 Change 1 0.210525 2 0.099142 3 0.067759 4 0.054542 5 0.042837 6 0.032883 7 0.025631 8 0.023001 9 0.019900 10 0.018673 11 0.017917 12 0.017467 13 0.016127 14 0.013993 15 0.015543 16 0.013003 17 0.014130 18 0.012016 19 0.012126 20 0.011274 21 0.011633 22 0.011271 23 0.010652 24 0.010747 25 0.010146 26 0.009564 27 0.009402 28 0.009200 29 0.009691 30 0.008143 31 0.008275 32 0.007882 33 0.008653 34 0.009487 35 0.008590 36 0.008691 37 0.009330 38 0.008251 39 0.008841 40 0.007711 41 0.008963 42 0.007795 43 0.007600 44 0.007990 45 0.008406 46 0.007192 47 0.007760 48 0.007286 49 0.007327 50 0.006926 51 0.008159 52 0.007471 53 0.007335 54 0.006701 55 0.006673 56 0.006119 57 0.006592 58 0.006030 59 0.005685 60 0.006859 61 0.006958 62 0.005861 63 0.006694 64 0.006738 65 0.006622 66 0.007069 67 0.006952 68 0.007605 69 0.007021 70 0.007059 71 0.006915 72 0.006776 73 0.007411 74 0.007005 75 0.006314 76 0.007138 77 0.006847 78 0.006947 79 0.007307 80 0.005963 81 0.006318 82 0.007091 83 0.006049 84 0.007715 85 0.006862 86 0.006332 87 0.006134 88 0.006678 89 0.006838 90 0.006252 91 0.008145 92 0.007249 93 0.007510 94 0.006169 95 0.006063 96 0.006987 97 0.006640 98 0.006606 99 0.007167 100 0.007457 Enter command:> generator_print Current cell generators: 1 0.4630535 0.6444063 2 0.6502456 0.3984759 3 0.7175637 0.7749521 4 0.7901528 0.1024039 5 0.5174735 0.4774783 6 0.9303898 0.5792974 7 0.2093864 0.8029273 8 7.1296662E-02 0.7510043 9 0.3728679 0.7783104 10 0.3768899 0.2267369 11 0.6504958 6.5307997E-02 12 0.4502512 0.9226047 13 0.9235168 0.4061159 14 0.9274930 7.4242920E-02 15 7.0187449E-02 8.8695399E-02 16 0.7788751 0.9181674 17 0.3785138 0.3867769 18 0.2483643 0.6567819 19 8.1266604E-02 0.5940009 20 0.5535886 0.7729222 21 0.3611849 0.5422465 22 0.9124699 0.7572802 23 0.9255747 0.2359052 24 0.9281847 0.9191272 25 0.2139150 6.3379720E-02 26 0.2169461 0.4991458 27 0.2740776 0.9327514 28 0.5082952 9.8641261E-02 29 0.2233484 0.3453193 30 0.7775891 0.4851415 31 0.2234242 0.1994024 32 7.3194399E-02 0.4251967 33 0.6166053 0.9241942 34 0.3619415 7.2675318E-02 35 7.3809125E-02 0.2582902 36 8.6538702E-02 0.9237027 37 0.7860638 0.3024726 38 0.5188829 0.2923788 39 0.6531212 0.2178069 40 0.6320236 0.6023210 41 0.7914342 0.6459860 Enter command:> maxit = 5 Enter command:> ns_cvt = 20000 Enter command:> iterate STEP L2 Change 1 0.004982 2 0.003903 3 0.003974 4 0.003925 5 0.003609 Enter command:> maxit = 5 Enter command:> ns_cvt = 80000 Enter command:> iterate STEP L2 Change 1 0.002487 2 0.002458 3 0.002138 4 0.001948 5 0.002477 Enter command:> maxit = 5 Enter command:> ns_cvt = 320000 Enter command:> iterate STEP L2 Change 1 0.002288 2 0.002104 3 0.002199 4 0.002319 5 0.002198 Enter command:> generator_write = romero_08_cvt.txt Enter command:> quit QUIT requested. CVT_MAIN Normal end of execution. March 9 2003 2:20:25.238 PM