March 9 2003 1:58:58.334 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 = 1250914285 Enter command:> # romero_04.inp # Command echo has been enabled. Enter command:> generator_read = romero_04_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.9848000 0.1912000 2 9.8700002E-02 0.7289000 3 0.6021000 0.9578000 4 0.4273000 8.0600001E-02 5 1.0400000E-02 0.4662000 6 0.3039000 0.1424000 7 2.7400000E-02 0.7955000 8 0.8571000 0.8227000 9 0.2106000 0.8901000 10 0.7093000 0.5072000 11 0.9364000 0.5697000 12 0.7627000 0.3508000 13 8.7600000E-02 0.6214000 14 0.1603000 0.9861000 15 0.5587000 0.1641000 16 0.9513000 0.6427000 17 0.8228000 0.6633000 18 0.6881000 0.4537000 19 0.2415000 0.2215000 20 0.6620000 0.2172000 21 0.3554000 0.4343000 22 0.6541000 3.4899998E-02 23 0.2841000 0.8505000 24 0.5291000 0.2507000 25 0.7371000 0.3193000 26 0.5788000 1.2700000E-02 27 0.4631000 0.7595000 28 0.3290000 0.8662000 29 0.4095000 0.3834000 30 0.7831000 0.9406000 31 0.3780000 0.4092000 32 0.1288000 0.2855000 33 0.1910000 0.1184000 34 0.6167000 0.6052000 35 4.9500000E-02 0.6946000 36 0.8398000 0.7354000 37 0.4715000 0.5150000 38 0.2653000 5.2800000E-02 39 0.8839000 0.3055000 40 0.5101000 0.5517000 41 0.9220000 0.9090000 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 1250914285 Enter command:> iterate STEP L2 Change 1 0.280452 2 0.096389 3 0.060058 4 0.044555 5 0.033636 6 0.024944 7 0.022040 8 0.018562 9 0.017200 10 0.016105 11 0.015539 12 0.014169 13 0.014107 14 0.014550 15 0.014500 16 0.013568 17 0.012675 18 0.011765 19 0.012914 20 0.012778 21 0.011952 22 0.011169 23 0.011099 24 0.009923 25 0.009803 26 0.009898 27 0.010405 28 0.009792 29 0.009298 30 0.010046 31 0.008576 32 0.008309 33 0.007704 34 0.008483 35 0.007788 36 0.007668 37 0.008974 38 0.008247 39 0.009392 40 0.007569 41 0.007384 42 0.008294 43 0.009195 44 0.008185 45 0.008020 46 0.008041 47 0.008610 48 0.007393 49 0.008407 50 0.007845 51 0.008405 52 0.007961 53 0.007210 54 0.006835 55 0.006373 56 0.007814 57 0.007160 58 0.006479 59 0.007578 60 0.008237 61 0.005703 62 0.005899 63 0.006454 64 0.007000 65 0.007220 66 0.006426 67 0.007122 68 0.006658 69 0.006576 70 0.007329 71 0.007339 72 0.007536 73 0.007142 74 0.005821 75 0.007386 76 0.006333 77 0.006874 78 0.007320 79 0.007381 80 0.006521 81 0.007009 82 0.007009 83 0.005775 84 0.006422 85 0.007393 86 0.007711 87 0.007517 88 0.006963 89 0.007167 90 0.006914 91 0.006304 92 0.006014 93 0.006943 94 0.006435 95 0.007491 96 0.006950 97 0.006966 98 0.006573 99 0.007100 100 0.006723 Enter command:> generator_print Current cell generators: 1 0.9239476 8.6596042E-02 2 0.1929251 0.8117155 3 0.6028308 0.9073046 4 0.4311132 6.7698233E-02 5 8.9753173E-02 0.4095405 6 0.3492571 0.2139063 7 6.7281947E-02 0.7404039 8 0.9270596 0.7530699 9 0.2670364 0.9347167 10 0.7346300 0.5262281 11 0.9176489 0.4291252 12 0.7672203 0.3682548 13 0.2225480 0.6620209 14 7.9481684E-02 0.9216610 15 0.5065969 0.1955093 16 0.9233457 0.5808728 17 0.7939791 0.6591073 18 0.5933503 0.4300533 19 0.1973395 0.2170863 20 0.7772641 0.2065729 21 0.2412262 0.5112888 22 0.7663006 6.6371284E-02 23 0.3472205 0.7886376 24 0.4536647 0.3385247 25 0.6346498 0.2690635 26 0.6106258 8.2139045E-02 27 0.5027242 0.7623701 28 0.4424039 0.9246820 29 0.4116275 0.4799834 30 0.7563273 0.9405376 31 0.2689581 0.3607478 32 6.4095721E-02 0.2519854 33 8.3380468E-02 7.9606250E-02 34 0.6485437 0.7156944 35 7.5959459E-02 0.5715011 36 0.7722163 0.8087168 37 0.3872174 0.6318194 38 0.2564354 7.0748664E-02 39 0.9249568 0.2617043 40 0.5555047 0.5846819 41 0.9206126 0.9187292 Enter command:> maxit = 5 Enter command:> ns_cvt = 20000 Enter command:> iterate STEP L2 Change 1 0.004964 2 0.003790 3 0.003555 4 0.003382 5 0.003344 Enter command:> maxit = 5 Enter command:> ns_cvt = 80000 Enter command:> iterate STEP L2 Change 1 0.002387 2 0.002141 3 0.001891 4 0.002129 5 0.001884 Enter command:> maxit = 5 Enter command:> ns_cvt = 320000 Enter command:> iterate STEP L2 Change 1 0.001305 2 0.001177 3 0.001056 4 0.001139 5 0.001188 Enter command:> generator_write = romero_04_cvt.txt Enter command:> quit QUIT requested. CVT_MAIN Normal end of execution. March 9 2003 2:01:17.247 PM