19 November 2008 10:47:07.735 AM LCVT_PRB FORTRAN90 version Test the routines in the LCVT library. TEST01 CVT computes a Centroidal Voronoi Tessellation. R8MAT_LATINIZE makes it a Latin Hypersquare. In this test, we vary the sampling used during the CVT Latin iteration. Spatial dimension DIM_NUM = 2 Number of generators = 25 Initial random number seed = 123456789 Initialize using UNIFORM. Sample using RANDOM_NUMBER (Fortran90 intrinsic). Number of sample points = 100000 Number of sample steps = 50 STEP L2 Change 1 0.314183 2 0.105687 3 0.058406 4 0.040435 5 0.036713 6 0.034276 7 0.030940 8 0.029897 9 0.028522 10 0.028552 11 0.026604 12 0.027358 13 0.025306 14 0.021496 15 0.020033 16 0.018962 17 0.017060 18 0.015925 19 0.017171 20 0.012546 21 0.011617 22 0.010869 23 0.010034 24 0.010054 25 0.010407 26 0.009610 27 0.010565 28 0.009392 29 0.008519 30 0.009487 31 0.008889 32 0.007914 33 0.008887 34 0.008316 35 0.009005 36 0.008333 37 0.007442 38 0.008174 39 0.007662 40 0.008046 41 0.009087 42 0.008584 43 0.007057 44 0.008774 45 0.007046 46 0.008245 47 0.007891 48 0.007473 49 0.007740 50 0.007771 After CVT steps: Row 1 2 Col 1 0.116817 0.908190 2 0.716517 0.746416 3 0.695443 0.114387 4 0.493020 0.834200E-01 5 0.911007E-01 0.694562 6 0.102220 0.488087 7 0.482681 0.723949 8 0.898285 0.882196E-01 9 0.917118 0.461623 10 0.975366E-01 0.911553E-01 11 0.881600 0.903418 12 0.294584 0.115630 13 0.386759 0.917646 14 0.285132 0.344748 15 0.887795 0.271131 16 0.598750 0.565622 17 0.767830 0.537991 18 0.912444 0.690863 19 0.274163 0.774139 20 0.486936 0.266599 21 0.458063 0.463130 22 0.682387 0.343377 23 0.296318 0.576262 24 0.996464E-01 0.281999 25 0.625310 0.911064 After Latin step: Row 1 2 Col 1 0.180000 0.900000 2 0.740000 0.780000 3 0.700000 0.140000 4 0.540000 0.200000E-01 5 0.200000E-01 0.700000 6 0.140000 0.500000 7 0.460000 0.740000 8 0.900000 0.600000E-01 9 0.980000 0.420000 10 0.600000E-01 0.100000 11 0.820000 0.860000 12 0.300000 0.180000 13 0.380000 0.980000 14 0.260000 0.380000 15 0.860000 0.260000 16 0.580000 0.580000 17 0.780000 0.540000 18 0.940000 0.660000 19 0.220000 0.820000 20 0.500000 0.220000 21 0.420000 0.460000 22 0.660000 0.340000 23 0.340000 0.620000 24 0.100000 0.300000 25 0.620000 0.940000 STEP L2 Change 1 0.337151 2 0.111009 3 0.073269 4 0.063906 5 0.062759 6 0.046974 7 0.036446 8 0.029860 9 0.025250 10 0.020119 11 0.019682 12 0.015881 13 0.014749 14 0.011907 15 0.013257 16 0.012845 17 0.011425 18 0.011250 19 0.011182 20 0.010909 21 0.011512 22 0.010854 23 0.008470 24 0.009438 25 0.009243 26 0.009021 27 0.009507 28 0.010373 29 0.008335 30 0.008316 31 0.008685 32 0.007039 33 0.008060 34 0.008188 35 0.007264 36 0.007575 37 0.009578 38 0.008643 39 0.006669 40 0.006597 41 0.007275 42 0.006869 43 0.007996 44 0.009404 45 0.007502 46 0.007763 47 0.008826 48 0.007537 49 0.008445 50 0.008775 After CVT steps: Row 1 2 Col 1 0.122602 0.702328 2 0.915037E-01 0.901543 3 0.895116 0.785327E-01 4 0.856899 0.743512 5 0.893510 0.243990 6 0.894786 0.418448 7 0.937399E-01 0.495058 8 0.103309 0.293907 9 0.286559 0.892491 10 0.953390E-01 0.959417E-01 11 0.285711 0.113534 12 0.478399 0.972857E-01 13 0.306777 0.332392 14 0.494038 0.899929 15 0.698167 0.321139 16 0.614764 0.705713 17 0.382958 0.711206 18 0.908320 0.910409 19 0.684573 0.107430 20 0.285045 0.536982 21 0.907785 0.592333 22 0.506237 0.292682 23 0.700107 0.902033 24 0.708125 0.529094 25 0.496685 0.504723 After Latin step: Row 1 2 Col 1 0.180000 0.660000 2 0.200000E-01 0.900000 3 0.900000 0.200000E-01 4 0.780000 0.780000 5 0.820000 0.220000 6 0.860000 0.420000 7 0.600000E-01 0.460000 8 0.140000 0.300000 9 0.300000 0.820000 10 0.100000 0.600000E-01 11 0.260000 0.180000 12 0.420000 0.100000 13 0.340000 0.380000 14 0.460000 0.860000 15 0.660000 0.340000 16 0.580000 0.700000 17 0.380000 0.740000 18 0.980000 0.980000 19 0.620000 0.140000 20 0.220000 0.580000 21 0.940000 0.620000 22 0.540000 0.260000 23 0.700000 0.940000 24 0.740000 0.540000 25 0.500000 0.500000 STEP L2 Change 1 0.274362 2 0.127961 3 0.084423 4 0.062268 5 0.050495 6 0.044570 7 0.040085 8 0.035842 9 0.032338 10 0.028954 11 0.029467 12 0.026143 13 0.023117 14 0.021518 15 0.018857 16 0.017298 17 0.017169 18 0.016243 19 0.015355 20 0.015219 21 0.015479 22 0.015159 23 0.017514 24 0.016777 25 0.017201 26 0.015579 27 0.016556 28 0.016335 29 0.015329 30 0.013315 31 0.013414 32 0.016153 33 0.012872 34 0.012844 35 0.013842 36 0.012674 37 0.012628 38 0.011948 39 0.011983 40 0.013368 41 0.012866 42 0.013743 43 0.013259 44 0.014654 45 0.013312 46 0.015289 47 0.014782 48 0.014391 49 0.014299 50 0.011353 After CVT steps: Row 1 2 Col 1 0.565964 0.906223 2 0.301598 0.361894 3 0.724182 0.310289 4 0.497612 0.280804 5 0.355491 0.926814 6 0.918176 0.885516 7 0.101663 0.901664E-01 8 0.253577 0.595193 9 0.430385 0.740202 10 0.753248 0.869185 11 0.903095 0.125109 12 0.101423 0.910902 13 0.629102 0.693293 14 0.905197 0.654478 15 0.104766 0.277666 16 0.967852E-01 0.475063 17 0.763355 0.544240 18 0.302185 0.122542 19 0.426212 0.536296 20 0.823116E-01 0.694394 21 0.240939 0.800035 22 0.502938 0.891662E-01 23 0.913026 0.384756 24 0.703045 0.103929 25 0.586213 0.468147 After Latin step: Row 1 2 Col 1 0.580000 0.900000 2 0.300000 0.340000 3 0.740000 0.300000 4 0.500000 0.260000 5 0.380000 0.980000 6 0.980000 0.860000 7 0.140000 0.600000E-01 8 0.260000 0.580000 9 0.460000 0.740000 10 0.780000 0.820000 11 0.860000 0.180000 12 0.100000 0.940000 13 0.660000 0.660000 14 0.900000 0.620000 15 0.180000 0.220000 16 0.600000E-01 0.460000 17 0.820000 0.540000 18 0.340000 0.140000 19 0.420000 0.500000 20 0.200000E-01 0.700000 21 0.220000 0.780000 22 0.540000 0.200000E-01 23 0.940000 0.380000 24 0.700000 0.100000 25 0.620000 0.420000 TEST01 CVT computes a Centroidal Voronoi Tessellation. R8MAT_LATINIZE makes it a Latin Hypersquare. In this test, we vary the sampling used during the CVT Latin iteration. Spatial dimension DIM_NUM = 2 Number of generators = 25 Initial random number seed = 123456789 Initialize using UNIFORM. Sample using UNIFORM. Number of sample points = 100000 Number of sample steps = 50 STEP L2 Change 1 0.312474 2 0.105229 3 0.057077 4 0.042855 5 0.035959 6 0.036197 7 0.036203 8 0.033884 9 0.031400 10 0.027388 11 0.025810 12 0.024076 13 0.021438 14 0.019048 15 0.017143 16 0.016397 17 0.015198 18 0.012561 19 0.010360 20 0.011249 21 0.010410 22 0.008952 23 0.008881 24 0.011511 25 0.010030 26 0.010919 27 0.010776 28 0.010629 29 0.010289 30 0.009319 31 0.011271 32 0.008357 33 0.008795 34 0.009263 35 0.008482 36 0.008563 37 0.008436 38 0.009120 39 0.009266 40 0.007812 41 0.008304 42 0.008500 43 0.007636 44 0.007843 45 0.008545 46 0.007899 47 0.007240 48 0.008271 49 0.007754 50 0.006911 After CVT steps: Row 1 2 Col 1 0.116933 0.907766 2 0.720610 0.776328 3 0.667343 0.813520E-01 4 0.465368 0.914709E-01 5 0.905720E-01 0.696740 6 0.102859 0.488649 7 0.494212 0.719049 8 0.895026 0.978081E-01 9 0.905019 0.305356 10 0.962695E-01 0.932955E-01 11 0.882181 0.908216 12 0.282569 0.122934 13 0.385606 0.915650 14 0.301766 0.352592 15 0.721570 0.235893 16 0.685709 0.596080 17 0.894163 0.502839 18 0.906633 0.701316 19 0.278585 0.769737 20 0.502472 0.279015 21 0.487803 0.495267 22 0.687004 0.411335 23 0.298109 0.574324 24 0.102286 0.287968 25 0.616493 0.915848 After Latin step: Row 1 2 Col 1 0.180000 0.860000 2 0.740000 0.820000 3 0.620000 0.200000E-01 4 0.420000 0.600000E-01 5 0.200000E-01 0.660000 6 0.140000 0.460000 7 0.500000 0.740000 8 0.900000 0.140000 9 0.940000 0.340000 10 0.600000E-01 0.100000 11 0.820000 0.900000 12 0.260000 0.180000 13 0.380000 0.940000 14 0.340000 0.380000 15 0.780000 0.220000 16 0.660000 0.620000 17 0.860000 0.540000 18 0.980000 0.700000 19 0.220000 0.780000 20 0.540000 0.260000 21 0.460000 0.500000 22 0.700000 0.420000 23 0.300000 0.580000 24 0.100000 0.300000 25 0.580000 0.980000 STEP L2 Change 1 0.272960 2 0.126839 3 0.076874 4 0.051343 5 0.035603 6 0.031433 7 0.029451 8 0.027206 9 0.027264 10 0.023962 11 0.023047 12 0.023884 13 0.024100 14 0.021692 15 0.019882 16 0.019437 17 0.017826 18 0.014613 19 0.014302 20 0.012283 21 0.011673 22 0.012704 23 0.012243 24 0.012476 25 0.012810 26 0.011923 27 0.010377 28 0.010558 29 0.010651 30 0.010645 31 0.010683 32 0.010321 33 0.009609 34 0.010294 35 0.011446 36 0.015424 37 0.015094 38 0.015264 39 0.015964 40 0.015499 41 0.011782 42 0.013059 43 0.011504 44 0.010421 45 0.011016 46 0.009460 47 0.010944 48 0.009618 49 0.009125 50 0.010022 After CVT steps: Row 1 2 Col 1 0.285544 0.278939 2 0.895245 0.912542 3 0.341895 0.684733 4 0.544049 0.542154 5 0.215774 0.517126 6 0.742361 0.251790 7 0.471365 0.909350 8 0.890686E-01 0.131367 9 0.914929 0.328355 10 0.900220 0.736261 11 0.887913E-01 0.889404 12 0.672194 0.892615E-01 13 0.271700 0.859692E-01 14 0.691970 0.896142 15 0.391401 0.446291 16 0.544097 0.739400 17 0.904779 0.536171 18 0.718872 0.446211 19 0.273236 0.885989 20 0.898767 0.104032 21 0.875564E-01 0.379518 22 0.981319E-01 0.668781 23 0.537617 0.301166 24 0.723253 0.666643 25 0.466457 0.108691 After Latin step: Row 1 2 Col 1 0.300000 0.260000 2 0.820000 0.980000 3 0.340000 0.700000 4 0.540000 0.580000 5 0.180000 0.500000 6 0.780000 0.220000 7 0.460000 0.940000 8 0.100000 0.180000 9 0.980000 0.340000 10 0.900000 0.740000 11 0.600000E-01 0.860000 12 0.620000 0.600000E-01 13 0.220000 0.200000E-01 14 0.660000 0.900000 15 0.380000 0.460000 16 0.580000 0.780000 17 0.940000 0.540000 18 0.700000 0.420000 19 0.260000 0.820000 20 0.860000 0.100000 21 0.200000E-01 0.380000 22 0.140000 0.660000 23 0.500000 0.300000 24 0.740000 0.620000 25 0.420000 0.140000 STEP L2 Change 1 0.272935 2 0.098665 3 0.064802 4 0.050033 5 0.039388 6 0.036052 7 0.031663 8 0.029825 9 0.026627 10 0.025127 11 0.020570 12 0.019641 13 0.018287 14 0.018153 15 0.015787 16 0.016466 17 0.016924 18 0.018527 19 0.017511 20 0.019138 21 0.019389 22 0.018982 23 0.017347 24 0.015818 25 0.014955 26 0.014053 27 0.014128 28 0.012366 29 0.012696 30 0.010807 31 0.011695 32 0.010395 33 0.010895 34 0.008767 35 0.010165 36 0.010192 37 0.008469 38 0.008233 39 0.008463 40 0.008500 41 0.007880 42 0.008425 43 0.009133 44 0.008529 45 0.009242 46 0.008489 47 0.008514 48 0.008559 49 0.006566 50 0.008250 After CVT steps: Row 1 2 Col 1 0.995147E-01 0.728372 2 0.100363 0.530052 3 0.724188 0.840686E-01 4 0.905822 0.703991 5 0.723838 0.810222 6 0.120211 0.909443 7 0.490135 0.528139 8 0.585502 0.922395 9 0.909339 0.102561 10 0.898478 0.501413 11 0.273181 0.262524 12 0.917902E-01 0.322945 13 0.345828 0.898368 14 0.504468 0.735563 15 0.298466 0.677859 16 0.890688 0.301609 17 0.316780 0.880483E-01 18 0.695661 0.615722 19 0.686815 0.428018 20 0.473327 0.311439 21 0.674163 0.252330 22 0.518982 0.992726E-01 23 0.293752 0.461662 24 0.885504 0.913514 25 0.101102 0.104949 After Latin step: Row 1 2 Col 1 0.600000E-01 0.740000 2 0.100000 0.580000 3 0.780000 0.200000E-01 4 0.940000 0.700000 5 0.740000 0.820000 6 0.180000 0.900000 7 0.460000 0.540000 8 0.580000 0.980000 9 0.980000 0.140000 10 0.900000 0.500000 11 0.220000 0.260000 12 0.200000E-01 0.380000 13 0.380000 0.860000 14 0.500000 0.780000 15 0.300000 0.660000 16 0.860000 0.300000 17 0.340000 0.600000E-01 18 0.700000 0.620000 19 0.660000 0.420000 20 0.420000 0.340000 21 0.620000 0.220000 22 0.540000 0.100000 23 0.260000 0.460000 24 0.820000 0.940000 25 0.140000 0.180000 TEST01 CVT computes a Centroidal Voronoi Tessellation. R8MAT_LATINIZE makes it a Latin Hypersquare. In this test, we vary the sampling used during the CVT Latin iteration. Spatial dimension DIM_NUM = 2 Number of generators = 25 Initial random number seed = 123456789 Initialize using UNIFORM. Sample using HALTON. Number of sample points = 100000 Number of sample steps = 50 STEP L2 Change 1 0.313709 2 0.103967 3 0.055809 4 0.040693 5 0.035751 6 0.035253 7 0.035543 8 0.032250 9 0.028541 10 0.026049 11 0.024934 12 0.022688 13 0.020879 14 0.018426 15 0.015865 16 0.013686 17 0.012200 18 0.010608 19 0.009798 20 0.008726 21 0.008004 22 0.007694 23 0.007017 24 0.006917 25 0.006937 26 0.006641 27 0.006681 28 0.006831 29 0.006402 30 0.006128 31 0.005799 32 0.005589 33 0.005190 34 0.004990 35 0.004845 36 0.004266 37 0.004225 38 0.004257 39 0.003994 40 0.003773 41 0.003626 42 0.003426 43 0.003309 44 0.002984 45 0.002830 46 0.002730 47 0.002524 48 0.002385 49 0.002135 50 0.002354 After CVT steps: Row 1 2 Col 1 0.117622 0.907556 2 0.722278 0.775014 3 0.666090 0.825889E-01 4 0.465020 0.919154E-01 5 0.928176E-01 0.695466 6 0.102324 0.487503 7 0.495350 0.717915 8 0.894006 0.963251E-01 9 0.907085 0.304569 10 0.947520E-01 0.928697E-01 11 0.883292 0.908789 12 0.282002 0.122723 13 0.383260 0.916528 14 0.302535 0.351410 15 0.723881 0.237409 16 0.685904 0.595619 17 0.893656 0.501826 18 0.906395 0.699860 19 0.279633 0.767809 20 0.503557 0.280666 21 0.487020 0.494807 22 0.688613 0.410508 23 0.297071 0.571729 24 0.102600 0.287418 25 0.616686 0.916086 After Latin step: Row 1 2 Col 1 0.180000 0.860000 2 0.740000 0.820000 3 0.620000 0.200000E-01 4 0.420000 0.600000E-01 5 0.200000E-01 0.660000 6 0.100000 0.460000 7 0.500000 0.740000 8 0.900000 0.140000 9 0.980000 0.340000 10 0.600000E-01 0.100000 11 0.820000 0.900000 12 0.260000 0.180000 13 0.380000 0.980000 14 0.340000 0.380000 15 0.780000 0.220000 16 0.660000 0.620000 17 0.860000 0.540000 18 0.940000 0.700000 19 0.220000 0.780000 20 0.540000 0.260000 21 0.460000 0.500000 22 0.700000 0.420000 23 0.300000 0.580000 24 0.140000 0.300000 25 0.580000 0.940000 STEP L2 Change 1 0.338451 2 0.111042 3 0.071632 4 0.062695 5 0.062697 6 0.047747 7 0.035373 8 0.027995 9 0.023427 10 0.019701 11 0.016764 12 0.014582 13 0.012791 14 0.011507 15 0.010297 16 0.009470 17 0.008881 18 0.008281 19 0.008058 20 0.007547 21 0.007097 22 0.006733 23 0.006493 24 0.005991 25 0.005650 26 0.005302 27 0.005048 28 0.004708 29 0.004559 30 0.004657 31 0.004276 32 0.004224 33 0.003948 34 0.003698 35 0.003409 36 0.003363 37 0.003340 38 0.003220 39 0.003159 40 0.002961 41 0.003000 42 0.002950 43 0.003138 44 0.003177 45 0.002974 46 0.003147 47 0.003378 48 0.003280 49 0.003646 50 0.003249 After CVT steps: Row 1 2 Col 1 0.122179 0.702410 2 0.931445E-01 0.899895 3 0.893460 0.796305E-01 4 0.857334 0.744741 5 0.893268 0.244376 6 0.896464 0.417668 7 0.939085E-01 0.494810 8 0.104933 0.292763 9 0.287457 0.893535 10 0.952161E-01 0.959487E-01 11 0.287731 0.114034 12 0.481452 0.951924E-01 13 0.307638 0.334347 14 0.495084 0.899305 15 0.696365 0.323631 16 0.616235 0.707412 17 0.385510 0.709564 18 0.906766 0.911068 19 0.684934 0.107623 20 0.285426 0.540002 21 0.909007 0.592897 22 0.504162 0.291479 23 0.700564 0.901106 24 0.708839 0.530095 25 0.497403 0.505024 After Latin step: Row 1 2 Col 1 0.180000 0.660000 2 0.200000E-01 0.900000 3 0.860000 0.200000E-01 4 0.780000 0.780000 5 0.820000 0.220000 6 0.900000 0.420000 7 0.600000E-01 0.460000 8 0.140000 0.300000 9 0.260000 0.820000 10 0.100000 0.100000 11 0.300000 0.180000 12 0.420000 0.600000E-01 13 0.340000 0.380000 14 0.460000 0.860000 15 0.660000 0.340000 16 0.580000 0.700000 17 0.380000 0.740000 18 0.940000 0.980000 19 0.620000 0.140000 20 0.220000 0.580000 21 0.980000 0.620000 22 0.540000 0.260000 23 0.700000 0.940000 24 0.740000 0.540000 25 0.500000 0.500000 STEP L2 Change 1 0.275100 2 0.125962 3 0.084490 4 0.061369 5 0.049675 6 0.043655 7 0.038759 8 0.035638 9 0.032584 10 0.030297 11 0.027832 12 0.025839 13 0.023771 14 0.022716 15 0.020979 16 0.019370 17 0.018606 18 0.017239 19 0.017526 20 0.017219 21 0.016785 22 0.015958 23 0.015021 24 0.014201 25 0.012860 26 0.011464 27 0.010576 28 0.009568 29 0.008941 30 0.008369 31 0.007910 32 0.007864 33 0.007210 34 0.007098 35 0.006875 36 0.006574 37 0.006066 38 0.005920 39 0.005397 40 0.005482 41 0.005180 42 0.004752 43 0.004448 44 0.004257 45 0.003773 46 0.003733 47 0.003163 48 0.003315 49 0.003470 50 0.002980 After CVT steps: Row 1 2 Col 1 0.600629 0.908951 2 0.330944 0.432433 3 0.734943 0.470829 4 0.379644 0.245260 5 0.418435 0.903613 6 0.925064 0.856260 7 0.106918 0.954984E-01 8 0.234344 0.597594 9 0.559145 0.731769 10 0.772310 0.907061 11 0.880381 0.989345E-01 12 0.834157E-01 0.904869 13 0.529296 0.539650 14 0.914645 0.594343 15 0.132777 0.280827 16 0.910098E-01 0.471195 17 0.745980 0.694047 18 0.334005 0.840030E-01 19 0.372356 0.692983 20 0.913593E-01 0.701655 21 0.247246 0.875084 22 0.589981 0.934230E-01 23 0.912363 0.343108 24 0.733458 0.239361 25 0.554582 0.335274 After Latin step: Row 1 2 Col 1 0.660000 0.980000 2 0.300000 0.380000 3 0.740000 0.420000 4 0.420000 0.220000 5 0.460000 0.860000 6 0.980000 0.780000 7 0.140000 0.100000 8 0.220000 0.580000 9 0.580000 0.740000 10 0.820000 0.940000 11 0.860000 0.140000 12 0.200000E-01 0.900000 13 0.500000 0.500000 14 0.940000 0.540000 15 0.180000 0.260000 16 0.600000E-01 0.460000 17 0.780000 0.660000 18 0.340000 0.200000E-01 19 0.380000 0.620000 20 0.100000 0.700000 21 0.260000 0.820000 22 0.620000 0.600000E-01 23 0.900000 0.340000 24 0.700000 0.180000 25 0.540000 0.300000 TEST01 CVT computes a Centroidal Voronoi Tessellation. R8MAT_LATINIZE makes it a Latin Hypersquare. In this test, we vary the sampling used during the CVT Latin iteration. Spatial dimension DIM_NUM = 2 Number of generators = 25 Initial random number seed = 123456789 Initialize using UNIFORM. Sample using GRID. Number of sample points = 100000 Number of sample steps = 50 STEP L2 Change 1 0.312960 2 0.104298 3 0.055621 4 0.040594 5 0.034914 6 0.032824 7 0.032171 8 0.031990 9 0.030718 10 0.028242 11 0.026117 12 0.023893 13 0.021650 14 0.018758 15 0.016188 16 0.013999 17 0.012129 18 0.010746 19 0.009725 20 0.008827 21 0.008075 22 0.007499 23 0.006967 24 0.006654 25 0.006487 26 0.006322 27 0.006361 28 0.006502 29 0.006233 30 0.006346 31 0.006021 32 0.005711 33 0.005192 34 0.005642 35 0.004846 36 0.004396 37 0.004265 38 0.004314 39 0.003901 40 0.003521 41 0.003080 42 0.003614 43 0.003306 44 0.002706 45 0.003068 46 0.002599 47 0.002207 48 0.002826 49 0.002486 50 0.002286 After CVT steps: Row 1 2 Col 1 0.117378 0.907721 2 0.716593 0.774318 3 0.664889 0.824716E-01 4 0.464305 0.915203E-01 5 0.923139E-01 0.696472 6 0.102332 0.490041 7 0.492804 0.717175 8 0.891355 0.966168E-01 9 0.904400 0.304509 10 0.945985E-01 0.934519E-01 11 0.877690 0.908594 12 0.281417 0.122941 13 0.381538 0.916179 14 0.301524 0.351511 15 0.721275 0.237280 16 0.682725 0.594358 17 0.888936 0.501227 18 0.900574 0.699775 19 0.278289 0.767238 20 0.501315 0.280026 21 0.484752 0.494163 22 0.685754 0.409354 23 0.296169 0.571788 24 0.102177 0.289789 25 0.613459 0.916256 After Latin step: Row 1 2 Col 1 0.180000 0.860000 2 0.740000 0.820000 3 0.620000 0.200000E-01 4 0.420000 0.600000E-01 5 0.200000E-01 0.660000 6 0.140000 0.460000 7 0.500000 0.740000 8 0.900000 0.140000 9 0.980000 0.340000 10 0.600000E-01 0.100000 11 0.820000 0.900000 12 0.260000 0.180000 13 0.380000 0.940000 14 0.340000 0.380000 15 0.780000 0.220000 16 0.660000 0.620000 17 0.860000 0.540000 18 0.940000 0.700000 19 0.220000 0.780000 20 0.540000 0.260000 21 0.460000 0.500000 22 0.700000 0.420000 23 0.300000 0.580000 24 0.100000 0.300000 25 0.580000 0.980000 STEP L2 Change 1 0.337894 2 0.111462 3 0.072592 4 0.062935 5 0.063232 6 0.047543 7 0.034878 8 0.028003 9 0.023101 10 0.019849 11 0.017135 12 0.014358 13 0.012781 14 0.011310 15 0.010804 16 0.009665 17 0.009454 18 0.008590 19 0.008148 20 0.007622 21 0.007083 22 0.006589 23 0.006321 24 0.005901 25 0.005801 26 0.005462 27 0.005540 28 0.004632 29 0.004703 30 0.004197 31 0.004944 32 0.003534 33 0.003270 34 0.003388 35 0.003354 36 0.002872 37 0.002667 38 0.002566 39 0.002422 40 0.002340 41 0.002156 42 0.001860 43 0.001482 44 0.001452 45 0.001262 46 0.001219 47 0.001317 48 0.001311 49 0.001310 50 0.001210 After CVT steps: Row 1 2 Col 1 0.122075 0.703468 2 0.930657E-01 0.900179 3 0.888235 0.800818E-01 4 0.862021 0.749471 5 0.897071 0.244681 6 0.894083 0.419968 7 0.920191E-01 0.495608 8 0.106279 0.294782 9 0.287415 0.894629 10 0.934350E-01 0.971009E-01 11 0.282783 0.111626 12 0.474567 0.100947 13 0.312718 0.328976 14 0.495081 0.897688 15 0.706713 0.307564 16 0.622998 0.702342 17 0.384262 0.710845 18 0.899488 0.914574 19 0.678491 0.103577 20 0.279234 0.538023 21 0.889518 0.593849 22 0.514808 0.302611 23 0.697008 0.899857 24 0.692992 0.510392 25 0.485879 0.511709 After Latin step: Row 1 2 Col 1 0.180000 0.700000 2 0.600000E-01 0.940000 3 0.820000 0.200000E-01 4 0.780000 0.780000 5 0.940000 0.220000 6 0.900000 0.420000 7 0.200000E-01 0.460000 8 0.140000 0.260000 9 0.300000 0.820000 10 0.100000 0.600000E-01 11 0.260000 0.180000 12 0.420000 0.100000 13 0.340000 0.380000 14 0.500000 0.860000 15 0.740000 0.340000 16 0.580000 0.660000 17 0.380000 0.740000 18 0.980000 0.980000 19 0.620000 0.140000 20 0.220000 0.580000 21 0.860000 0.620000 22 0.540000 0.300000 23 0.700000 0.900000 24 0.660000 0.500000 25 0.460000 0.540000 STEP L2 Change 1 0.274582 2 0.125347 3 0.084359 4 0.061432 5 0.050010 6 0.043372 7 0.038854 8 0.035330 9 0.032662 10 0.029907 11 0.027778 12 0.026349 13 0.024574 14 0.022722 15 0.020957 16 0.019368 17 0.018329 18 0.017544 19 0.017549 20 0.017527 21 0.016893 22 0.015997 23 0.014815 24 0.013648 25 0.012314 26 0.010988 27 0.010225 28 0.009523 29 0.008798 30 0.007962 31 0.007651 32 0.007360 33 0.006952 34 0.006874 35 0.006825 36 0.006473 37 0.006149 38 0.006036 39 0.005611 40 0.005354 41 0.004849 42 0.004866 43 0.004555 44 0.004354 45 0.004083 46 0.003974 47 0.003932 48 0.003746 49 0.003614 50 0.003649 After CVT steps: Row 1 2 Col 1 0.596934 0.907657 2 0.329658 0.425128 3 0.729610 0.472162 4 0.381340 0.241921 5 0.417100 0.903910 6 0.919375 0.856502 7 0.106154 0.953542E-01 8 0.238339 0.583877 9 0.552540 0.729923 10 0.767377 0.907234 11 0.878171 0.994474E-01 12 0.822801E-01 0.903837 13 0.520696 0.540445 14 0.908780 0.595419 15 0.132116 0.280078 16 0.891819E-01 0.470711 17 0.740222 0.694383 18 0.331740 0.836677E-01 19 0.360906 0.699381 20 0.956686E-01 0.700970 21 0.245724 0.880254 22 0.589355 0.930419E-01 23 0.908868 0.344603 24 0.731338 0.240013 25 0.553146 0.335697 After Latin step: Row 1 2 Col 1 0.660000 0.980000 2 0.300000 0.380000 3 0.700000 0.460000 4 0.420000 0.220000 5 0.460000 0.900000 6 0.980000 0.780000 7 0.140000 0.100000 8 0.220000 0.540000 9 0.540000 0.740000 10 0.820000 0.940000 11 0.860000 0.140000 12 0.200000E-01 0.860000 13 0.500000 0.500000 14 0.900000 0.580000 15 0.180000 0.260000 16 0.600000E-01 0.420000 17 0.780000 0.620000 18 0.340000 0.200000E-01 19 0.380000 0.660000 20 0.100000 0.700000 21 0.260000 0.820000 22 0.620000 0.600000E-01 23 0.940000 0.340000 24 0.740000 0.180000 25 0.580000 0.300000 TEST02 CVT computes a Centroidal Voronoi Tessellation. R8MAT_LATINIZE makes it a Latin Hypersquare. In this test, we initialize the generators to grid points; this is an unstable CVT solution. Spatial dimension DIM_NUM = 2 Number of generators = 25 Initial random number seed = 123456789 USER will initialize data. Sample using UNIFORM. Number of sample points = 100000 Number of sample steps = 50 Initial generators (rows): Row 1 2 Col 1 0.100000 0.100000 2 0.100000 0.300000 3 0.100000 0.500000 4 0.100000 0.700000 5 0.100000 0.900000 6 0.300000 0.100000 7 0.300000 0.300000 8 0.300000 0.500000 9 0.300000 0.700000 10 0.300000 0.900000 11 0.500000 0.100000 12 0.500000 0.300000 13 0.500000 0.500000 14 0.500000 0.700000 15 0.500000 0.900000 16 0.700000 0.100000 17 0.700000 0.300000 18 0.700000 0.500000 19 0.700000 0.700000 20 0.700000 0.900000 21 0.900000 0.100000 22 0.900000 0.300000 23 0.900000 0.500000 24 0.900000 0.700000 25 0.900000 0.900000 STEP L2 Change 1 0.006640 2 0.009399 3 0.008002 4 0.008833 5 0.007833 6 0.006963 7 0.007560 8 0.007965 9 0.008468 10 0.006457 11 0.007293 12 0.006824 13 0.008566 14 0.007811 15 0.008052 16 0.009991 17 0.008560 18 0.008213 19 0.009217 20 0.007873 21 0.007564 22 0.008067 23 0.009780 24 0.008688 25 0.008935 26 0.009017 27 0.009542 28 0.009288 29 0.009209 30 0.007959 31 0.008316 32 0.007808 33 0.009434 34 0.009190 35 0.008586 36 0.008384 37 0.008312 38 0.008791 39 0.008651 40 0.007653 41 0.008112 42 0.006694 43 0.006798 44 0.007459 45 0.007248 46 0.006619 47 0.007639 48 0.007852 49 0.006549 50 0.007713 After CVT steps: Row 1 2 Col 1 0.102250 0.916803E-01 2 0.101744 0.279880 3 0.100881 0.477441 4 0.994338E-01 0.684091 5 0.990910E-01 0.893837 6 0.299120 0.112164 7 0.300636 0.330854 8 0.300272 0.542215 9 0.303680 0.738593 10 0.302822 0.916638 11 0.499432 0.821377E-01 12 0.502101 0.260051 13 0.498682 0.458127 14 0.500214 0.665765 15 0.500996 0.886632 16 0.700731 0.111268 17 0.698087 0.330361 18 0.697310 0.537542 19 0.694169 0.735810 20 0.703735 0.916048 21 0.900549 0.929568E-01 22 0.896097 0.281587 23 0.898170 0.480598 24 0.899424 0.684301 25 0.901173 0.891040 After Latin step: Row 1 2 Col 1 0.180000 0.600000E-01 2 0.140000 0.260000 3 0.100000 0.460000 4 0.600000E-01 0.660000 5 0.200000E-01 0.900000 6 0.220000 0.180000 7 0.300000 0.380000 8 0.260000 0.580000 9 0.380000 0.780000 10 0.340000 0.980000 11 0.460000 0.200000E-01 12 0.580000 0.220000 13 0.420000 0.420000 14 0.500000 0.620000 15 0.540000 0.820000 16 0.740000 0.140000 17 0.700000 0.340000 18 0.660000 0.540000 19 0.620000 0.740000 20 0.780000 0.940000 21 0.940000 0.100000 22 0.820000 0.300000 23 0.860000 0.500000 24 0.900000 0.700000 25 0.980000 0.860000 STEP L2 Change 1 0.179457 2 0.048693 3 0.027182 4 0.022060 5 0.019375 6 0.019138 7 0.017794 8 0.016536 9 0.015006 10 0.013615 11 0.011072 12 0.011383 13 0.010875 14 0.009622 15 0.010147 16 0.008503 17 0.008411 18 0.007133 19 0.007809 20 0.008277 21 0.007458 22 0.007988 23 0.008177 24 0.007856 25 0.009038 26 0.007996 27 0.007334 28 0.008354 29 0.008601 30 0.008025 31 0.007787 32 0.006872 33 0.007558 34 0.008464 35 0.007768 36 0.008107 37 0.008020 38 0.007339 39 0.006743 40 0.007461 41 0.007451 42 0.006998 43 0.007824 44 0.007274 45 0.007628 46 0.007074 47 0.008236 48 0.007175 49 0.007413 50 0.008998 After CVT steps: Row 1 2 Col 1 0.999069E-01 0.909814E-01 2 0.103280 0.277088 3 0.102898 0.476563 4 0.100763 0.685865 5 0.981655E-01 0.893570 6 0.297349 0.113582 7 0.298684 0.336336 8 0.298170 0.546245 9 0.302323 0.741549 10 0.302317 0.918583 11 0.495389 0.842624E-01 12 0.496858 0.262048 13 0.496403 0.459087 14 0.499857 0.667491 15 0.502631 0.887626 16 0.697715 0.112836 17 0.696082 0.331579 18 0.695568 0.539783 19 0.698104 0.735292 20 0.701825 0.915449 21 0.898914 0.924848E-01 22 0.897278 0.278430 23 0.898221 0.479599 24 0.899695 0.683245 25 0.899564 0.891528 After Latin step: Row 1 2 Col 1 0.600000E-01 0.600000E-01 2 0.180000 0.260000 3 0.140000 0.460000 4 0.100000 0.700000 5 0.200000E-01 0.900000 6 0.220000 0.180000 7 0.300000 0.380000 8 0.260000 0.580000 9 0.380000 0.780000 10 0.340000 0.980000 11 0.420000 0.200000E-01 12 0.500000 0.220000 13 0.460000 0.420000 14 0.540000 0.620000 15 0.580000 0.820000 16 0.700000 0.140000 17 0.660000 0.340000 18 0.620000 0.540000 19 0.740000 0.740000 20 0.780000 0.940000 21 0.900000 0.100000 22 0.820000 0.300000 23 0.860000 0.500000 24 0.980000 0.660000 25 0.940000 0.860000 STEP L2 Change 1 0.155244 2 0.054002 3 0.033023 4 0.024398 5 0.022571 6 0.016987 7 0.015616 8 0.012999 9 0.012353 10 0.011060 11 0.010695 12 0.009098 13 0.008696 14 0.008236 15 0.007908 16 0.007994 17 0.008635 18 0.008526 19 0.008765 20 0.008763 21 0.007843 22 0.007138 23 0.007547 24 0.007710 25 0.008041 26 0.009164 27 0.006964 28 0.007545 29 0.008365 30 0.007492 31 0.007475 32 0.008441 33 0.007863 34 0.006969 35 0.008090 36 0.009244 37 0.008215 38 0.007044 39 0.008055 40 0.006840 41 0.007261 42 0.005878 43 0.007555 44 0.007496 45 0.007376 46 0.007397 47 0.007899 48 0.006724 49 0.007871 50 0.006585 After CVT steps: Row 1 2 Col 1 0.996875E-01 0.920661E-01 2 0.104555 0.280829 3 0.100363 0.478226 4 0.100379 0.682634 5 0.999222E-01 0.893572 6 0.299359 0.112905 7 0.302327 0.333920 8 0.299105 0.543555 9 0.303829 0.739018 10 0.301184 0.916063 11 0.500746 0.817582E-01 12 0.502075 0.261845 13 0.497662 0.460535 14 0.501193 0.666995 15 0.501060 0.888249 16 0.700532 0.112159 17 0.695906 0.335459 18 0.696993 0.540627 19 0.698546 0.738937 20 0.700210 0.915012 21 0.898977 0.905614E-01 22 0.895290 0.277769 23 0.898986 0.475041 24 0.900330 0.679836 25 0.900599 0.894167 After Latin step: Row 1 2 Col 1 0.200000E-01 0.100000 2 0.180000 0.300000 3 0.100000 0.500000 4 0.140000 0.700000 5 0.600000E-01 0.860000 6 0.260000 0.180000 7 0.340000 0.340000 8 0.220000 0.580000 9 0.380000 0.780000 10 0.300000 0.980000 11 0.460000 0.200000E-01 12 0.580000 0.220000 13 0.420000 0.420000 14 0.540000 0.620000 15 0.500000 0.820000 16 0.780000 0.140000 17 0.620000 0.380000 18 0.660000 0.540000 19 0.700000 0.740000 20 0.740000 0.940000 21 0.860000 0.600000E-01 22 0.820000 0.260000 23 0.900000 0.460000 24 0.940000 0.660000 25 0.980000 0.900000 LCVT_PRB Normal end of execution. 19 November 2008 10:52:29.357 AM