March 9 2003 1:56:55.675 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 = 1249707173 Enter command:> # romero_02.inp # Command echo has been enabled. Enter command:> generator_read = romero_02_lhs.txt FILE_ROW_COUNT: Number of records: 17 Number of data records: 13 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.4061000 0.1210000 2 0.6642000 0.9451000 3 0.2763000 0.8988000 4 0.5055000 0.7620000 5 0.6024000 0.1923000 6 0.3112000 3.7599999E-02 7 0.7617000 0.2914000 8 0.9583000 0.6686000 9 0.8960000 0.4869000 10 3.2900002E-02 0.8255000 11 0.1560000 0.6087000 12 8.6599998E-02 0.4604000 13 0.8368000 0.3476000 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: 13 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 1249707173 Enter command:> iterate STEP L2 Change 1 0.269785 2 0.113927 3 0.066758 4 0.048904 5 0.036774 6 0.028089 7 0.025629 8 0.028488 9 0.024767 10 0.024002 11 0.020987 12 0.018641 13 0.019695 14 0.013029 15 0.013803 16 0.011051 17 0.009105 18 0.009607 19 0.007950 20 0.009469 21 0.007378 22 0.007104 23 0.006127 24 0.007846 25 0.007496 26 0.008545 27 0.006023 28 0.007350 29 0.007936 30 0.007267 31 0.007207 32 0.007381 33 0.005415 34 0.006420 35 0.008991 36 0.008217 37 0.007816 38 0.009256 39 0.009042 40 0.006118 41 0.005798 42 0.006330 43 0.007105 44 0.006426 45 0.006840 46 0.006195 47 0.006815 48 0.005988 49 0.005669 50 0.007280 51 0.007429 52 0.006766 53 0.007302 54 0.005736 55 0.006529 56 0.007322 57 0.007096 58 0.007400 59 0.006831 60 0.006931 61 0.007015 62 0.006570 63 0.008896 64 0.007025 65 0.004215 66 0.007217 67 0.004252 68 0.007352 69 0.005669 70 0.005834 71 0.008032 72 0.007556 73 0.008233 74 0.006665 75 0.006712 76 0.006311 77 0.006663 78 0.007503 79 0.006772 80 0.008032 81 0.007427 82 0.005459 83 0.007570 84 0.005265 85 0.006956 86 0.005571 87 0.006238 88 0.006377 89 0.006489 90 0.007294 91 0.006926 92 0.007804 93 0.006248 94 0.006769 95 0.006609 96 0.005468 97 0.006717 98 0.007516 99 0.007747 100 0.006341 Enter command:> generator_print Current cell generators: 1 0.3819475 0.1659886 2 0.6162442 0.8370537 3 0.3824911 0.8344610 4 0.5005670 0.4966779 5 0.6193532 0.1672953 6 0.1249032 0.1253436 7 0.8758508 0.1251729 8 0.8744751 0.8736492 9 0.8292465 0.6204013 10 0.1268599 0.8730382 11 0.1671565 0.6186774 12 0.1646375 0.3804355 13 0.8365548 0.3816557 Enter command:> maxit = 5 Enter command:> ns_cvt = 20000 Enter command:> iterate STEP L2 Change 1 0.004597 2 0.004078 3 0.003554 4 0.003826 5 0.003974 Enter command:> maxit = 5 Enter command:> ns_cvt = 80000 Enter command:> iterate STEP L2 Change 1 0.002104 2 0.001483 3 0.001605 4 0.001472 5 0.001476 Enter command:> maxit = 5 Enter command:> ns_cvt = 320000 Enter command:> iterate STEP L2 Change 1 0.000847 2 0.000830 3 0.000827 4 0.000905 5 0.000871 Enter command:> generator_write = romero_02_cvt.txt Enter command:> quit QUIT requested. CVT_MAIN Normal end of execution. March 9 2003 1:57:23.715 PM