22 December 2011 03:59:08 PM LCVT_PRB C++ version Test 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 M = 2 Number of generators = 25 Initial random number seed = 123456789 Initialize using UNIFORM. Sample using RANDOM_NUMBER (C++ STDLIB intrinsic). Number of sample points = 100000 Number of sample steps = 50 STEP L2 Change 1 0.314024 2 0.104015 3 0.0556271 4 0.0422177 5 0.0391032 6 0.0380846 7 0.0346962 8 0.0297201 9 0.0253493 10 0.0243533 11 0.0220516 12 0.0211193 13 0.0203718 14 0.0187989 15 0.0166751 16 0.0154369 17 0.0143212 18 0.0131501 19 0.0119631 20 0.0106695 21 0.00941972 22 0.00892315 23 0.00775858 24 0.00662621 25 0.00627641 26 0.00531321 27 0.00518186 28 0.00478021 29 0.00388683 30 0.00360003 31 0.00326982 32 0.00324899 33 0.00314693 34 0.00260767 35 0.00243986 36 0.00253931 37 0.00245512 38 0.00282205 39 0.00297869 40 0.00333447 41 0.00339234 42 0.00353608 43 0.00337712 44 0.00382045 45 0.00433563 46 0.00458005 47 0.00414493 48 0.00467738 49 0.00466608 50 0.00472917 After CVT steps: Row: 1 2 Col 1 0.114474 0.913193 2 0.722008 0.753118 3 0.667585 0.0890559 4 0.464982 0.0856054 5 0.0927853 0.704644 6 0.098652 0.496371 7 0.492229 0.716044 8 0.896326 0.0966554 9 0.912122 0.310878 10 0.0957155 0.0919195 11 0.885337 0.905274 12 0.28031 0.124556 13 0.395064 0.918826 14 0.285555 0.36128 15 0.737512 0.243865 16 0.668518 0.577124 17 0.875828 0.500078 18 0.910782 0.691908 19 0.277595 0.786783 20 0.483598 0.270271 21 0.46827 0.491574 22 0.656306 0.393699 23 0.290432 0.587407 24 0.098498 0.288998 25 0.624247 0.91134 After Latin step: Row: 1 2 Col 1 0.18 0.94 2 0.74 0.78 3 0.66 0.06 4 0.42 0.02 5 0.02 0.7 6 0.14 0.5 7 0.54 0.74 8 0.9 0.14 9 0.98 0.34 10 0.06 0.1 11 0.86 0.86 12 0.26 0.18 13 0.38 0.98 14 0.3 0.38 15 0.78 0.22 16 0.7 0.58 17 0.82 0.54 18 0.94 0.66 19 0.22 0.82 20 0.5 0.26 21 0.46 0.46 22 0.62 0.42 23 0.34 0.62 24 0.1 0.3 25 0.58 0.9 STEP L2 Change 1 0.339509 2 0.111239 3 0.0710747 4 0.0627591 5 0.0643316 6 0.0466494 7 0.0340442 8 0.0273164 9 0.023101 10 0.0194435 11 0.0169677 12 0.0155466 13 0.0138986 14 0.0131601 15 0.0115912 16 0.0102181 17 0.00952871 18 0.00843081 19 0.0083268 20 0.00768556 21 0.00749585 22 0.00713794 23 0.00692038 24 0.00674538 25 0.00632723 26 0.00588079 27 0.00533538 28 0.00509968 29 0.00483437 30 0.00471372 31 0.00418014 32 0.0039421 33 0.00405923 34 0.00368943 35 0.00377881 36 0.00318975 37 0.00311377 38 0.00325536 39 0.00319766 40 0.00274708 41 0.00316549 42 0.00336729 43 0.00322782 44 0.00330001 45 0.00328479 46 0.00317922 47 0.00342192 48 0.00337002 49 0.00321616 50 0.00345958 After CVT steps: Row: 1 2 Col 1 0.123732 0.709349 2 0.0941693 0.902837 3 0.885624 0.0783171 4 0.851621 0.741092 5 0.904877 0.239181 6 0.898045 0.414985 7 0.0887427 0.507177 8 0.111946 0.30848 9 0.289145 0.894101 10 0.0912471 0.101402 11 0.277898 0.109537 12 0.470752 0.110103 13 0.329303 0.320763 14 0.499249 0.899001 15 0.719835 0.30326 16 0.616856 0.706506 17 0.378886 0.710431 18 0.912002 0.905562 19 0.673381 0.105113 20 0.278012 0.528873 21 0.910461 0.591343 22 0.532944 0.316541 23 0.705517 0.90205 24 0.704416 0.51538 25 0.487387 0.518311 After Latin step: Row: 1 2 Col 1 0.18 0.7 2 0.1 0.94 3 0.82 0.02 4 0.78 0.78 5 0.9 0.22 6 0.86 0.42 7 0.02 0.46 8 0.14 0.3 9 0.3 0.82 10 0.06 0.06 11 0.22 0.14 12 0.42 0.18 13 0.34 0.38 14 0.5 0.86 15 0.74 0.26 16 0.58 0.66 17 0.38 0.74 18 0.98 0.98 19 0.62 0.1 20 0.26 0.58 21 0.94 0.62 22 0.54 0.34 23 0.7 0.9 24 0.66 0.5 25 0.46 0.54 STEP L2 Change 1 0.274452 2 0.126967 3 0.0851309 4 0.0615766 5 0.0491243 6 0.0424887 7 0.0373852 8 0.0332278 9 0.0305374 10 0.0279289 11 0.0254533 12 0.0237418 13 0.0231961 14 0.022798 15 0.0217271 16 0.0218066 17 0.0198424 18 0.0181839 19 0.0165055 20 0.0155802 21 0.0152034 22 0.0149855 23 0.0137563 24 0.0128151 25 0.0123969 26 0.0123686 27 0.011598 28 0.0110797 29 0.0102432 30 0.00966663 31 0.00900827 32 0.0083191 33 0.00829738 34 0.0078154 35 0.00718193 36 0.0067152 37 0.0061264 38 0.00537035 39 0.00550782 40 0.00510884 41 0.00518116 42 0.00480782 43 0.00471767 44 0.00468792 45 0.00455238 46 0.0047909 47 0.00496679 48 0.00504358 49 0.00557483 50 0.00587104 After CVT steps: Row: 1 2 Col 1 0.609972 0.900033 2 0.297042 0.413628 3 0.747992 0.460833 4 0.418406 0.243023 5 0.417133 0.909784 6 0.928946 0.850889 7 0.10945 0.0926097 8 0.237431 0.639115 9 0.512918 0.745736 10 0.784029 0.917547 11 0.875162 0.101076 12 0.0921883 0.907604 13 0.597105 0.579852 14 0.914615 0.585138 15 0.125394 0.27651 16 0.0889561 0.480364 17 0.751875 0.714646 18 0.329439 0.0929158 19 0.397801 0.604206 20 0.0810445 0.705362 21 0.259521 0.860407 22 0.593205 0.0857605 23 0.90872 0.329221 24 0.696606 0.249765 25 0.526129 0.399621 After Latin step: Row: 1 2 Col 1 0.66 0.86 2 0.3 0.38 3 0.74 0.42 4 0.46 0.18 5 0.42 0.94 6 0.98 0.78 7 0.14 0.06 8 0.22 0.62 9 0.5 0.74 10 0.82 0.98 11 0.86 0.14 12 0.1 0.9 13 0.62 0.5 14 0.94 0.54 15 0.18 0.26 16 0.06 0.46 17 0.78 0.7 18 0.34 0.1 19 0.38 0.58 20 0.02 0.66 21 0.26 0.82 22 0.58 0.02 23 0.9 0.3 24 0.7 0.22 25 0.54 0.34 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 M = 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.0570768 4 0.0428552 5 0.0359586 6 0.0361972 7 0.0362033 8 0.0338843 9 0.0314002 10 0.0273877 11 0.0258102 12 0.0240761 13 0.0214378 14 0.0190485 15 0.0171434 16 0.0163974 17 0.015198 18 0.0125611 19 0.0103603 20 0.0112487 21 0.0104096 22 0.00895235 23 0.00888119 24 0.011511 25 0.0100304 26 0.010919 27 0.0107757 28 0.0106287 29 0.0102891 30 0.00931901 31 0.0112714 32 0.00835721 33 0.00879489 34 0.0092627 35 0.00848223 36 0.00856298 37 0.00843568 38 0.00912026 39 0.00926614 40 0.00781166 41 0.00830438 42 0.00850008 43 0.00763605 44 0.0078426 45 0.00854515 46 0.00789935 47 0.0072398 48 0.00827137 49 0.00775375 50 0.00691107 After CVT steps: Row: 1 2 Col 1 0.116933 0.907766 2 0.72061 0.776328 3 0.667343 0.081352 4 0.465368 0.0914709 5 0.090572 0.69674 6 0.102859 0.488649 7 0.494212 0.719049 8 0.895026 0.0978081 9 0.905019 0.305356 10 0.0962695 0.0932955 11 0.882181 0.908216 12 0.282569 0.122934 13 0.385606 0.91565 14 0.301766 0.352592 15 0.72157 0.235893 16 0.685709 0.59608 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.18 0.86 2 0.74 0.82 3 0.62 0.02 4 0.42 0.06 5 0.02 0.66 6 0.14 0.46 7 0.5 0.74 8 0.9 0.14 9 0.94 0.34 10 0.06 0.1 11 0.82 0.9 12 0.26 0.18 13 0.38 0.94 14 0.34 0.38 15 0.78 0.22 16 0.66 0.62 17 0.86 0.54 18 0.98 0.7 19 0.22 0.78 20 0.54 0.26 21 0.46 0.5 22 0.7 0.42 23 0.3 0.58 24 0.1 0.3 25 0.58 0.98 STEP L2 Change 1 0.27296 2 0.126839 3 0.076874 4 0.0513427 5 0.0356027 6 0.0314327 7 0.0294513 8 0.0272059 9 0.0272645 10 0.0239617 11 0.0230466 12 0.0238836 13 0.0241003 14 0.0216922 15 0.0198821 16 0.0194373 17 0.0178256 18 0.0146131 19 0.0143019 20 0.0122827 21 0.0116732 22 0.0127043 23 0.0122433 24 0.0124757 25 0.01281 26 0.0119234 27 0.0103774 28 0.0105575 29 0.010651 30 0.0106451 31 0.0106835 32 0.0103209 33 0.00960925 34 0.0102941 35 0.0114456 36 0.0154237 37 0.0150937 38 0.015264 39 0.015964 40 0.0154994 41 0.0117824 42 0.0130586 43 0.0115037 44 0.010421 45 0.0110161 46 0.00945979 47 0.0109435 48 0.00961757 49 0.00912484 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.25179 7 0.471365 0.90935 8 0.0890686 0.131367 9 0.914929 0.328355 10 0.90022 0.736261 11 0.0887913 0.889404 12 0.672194 0.0892615 13 0.2717 0.0859692 14 0.69197 0.896142 15 0.391401 0.446291 16 0.544097 0.7394 17 0.904779 0.536171 18 0.718872 0.446211 19 0.273236 0.885989 20 0.898767 0.104032 21 0.0875564 0.379518 22 0.0981319 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.3 0.26 2 0.82 0.98 3 0.34 0.7 4 0.54 0.58 5 0.18 0.5 6 0.78 0.22 7 0.46 0.94 8 0.1 0.18 9 0.98 0.34 10 0.9 0.74 11 0.06 0.86 12 0.62 0.06 13 0.22 0.02 14 0.66 0.9 15 0.38 0.46 16 0.58 0.78 17 0.94 0.54 18 0.7 0.42 19 0.26 0.82 20 0.86 0.1 21 0.02 0.38 22 0.14 0.66 23 0.5 0.3 24 0.74 0.62 25 0.42 0.14 STEP L2 Change 1 0.272935 2 0.0986654 3 0.0648024 4 0.0500333 5 0.0393882 6 0.0360516 7 0.0316631 8 0.0298251 9 0.0266269 10 0.0251273 11 0.0205702 12 0.0196412 13 0.0182869 14 0.0181531 15 0.0157867 16 0.0164663 17 0.0169245 18 0.0185267 19 0.0175109 20 0.0191375 21 0.0193889 22 0.0189825 23 0.0173466 24 0.0158185 25 0.0149552 26 0.014053 27 0.0141278 28 0.0123657 29 0.012696 30 0.0108069 31 0.0116952 32 0.0103947 33 0.010895 34 0.00876674 35 0.0101645 36 0.0101923 37 0.00846925 38 0.00823281 39 0.00846332 40 0.00849974 41 0.00788009 42 0.00842539 43 0.00913292 44 0.00852903 45 0.00924226 46 0.00848935 47 0.00851411 48 0.00855905 49 0.00656634 50 0.00825004 After CVT steps: Row: 1 2 Col 1 0.0995147 0.728372 2 0.100363 0.530052 3 0.724188 0.0840686 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.0917902 0.322945 13 0.345828 0.898368 14 0.504468 0.735563 15 0.298466 0.677859 16 0.890688 0.301609 17 0.31678 0.0880483 18 0.695661 0.615722 19 0.686815 0.428018 20 0.473327 0.311439 21 0.674163 0.25233 22 0.518982 0.0992726 23 0.293752 0.461662 24 0.885504 0.913514 25 0.101102 0.104949 After Latin step: Row: 1 2 Col 1 0.06 0.74 2 0.1 0.58 3 0.78 0.02 4 0.94 0.7 5 0.74 0.82 6 0.18 0.9 7 0.46 0.54 8 0.58 0.98 9 0.98 0.14 10 0.9 0.5 11 0.22 0.26 12 0.02 0.38 13 0.38 0.86 14 0.5 0.78 15 0.3 0.66 16 0.86 0.3 17 0.34 0.06 18 0.7 0.62 19 0.66 0.42 20 0.42 0.34 21 0.62 0.22 22 0.54 0.1 23 0.26 0.46 24 0.82 0.94 25 0.14 0.18 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 M = 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.0406927 5 0.0357508 6 0.035253 7 0.0355435 8 0.03225 9 0.028541 10 0.0260489 11 0.0249343 12 0.0226881 13 0.0208791 14 0.0184261 15 0.0158646 16 0.0136864 17 0.0121998 18 0.0106075 19 0.00979826 20 0.0087259 21 0.00800402 22 0.00769355 23 0.0070173 24 0.00691705 25 0.00693696 26 0.0066406 27 0.0066807 28 0.00683059 29 0.00640193 30 0.00612838 31 0.00579879 32 0.00558938 33 0.00518958 34 0.00498987 35 0.00484544 36 0.00426551 37 0.00422506 38 0.00425657 39 0.00399379 40 0.00377305 41 0.0036263 42 0.00342603 43 0.00330895 44 0.00298417 45 0.00283019 46 0.00272989 47 0.00252374 48 0.00238476 49 0.00213491 50 0.0023535 After CVT steps: Row: 1 2 Col 1 0.117622 0.907556 2 0.722278 0.775014 3 0.66609 0.0825889 4 0.46502 0.0919154 5 0.0928176 0.695466 6 0.102324 0.487503 7 0.49535 0.717915 8 0.894006 0.0963251 9 0.907085 0.304569 10 0.094752 0.0928697 11 0.883292 0.908789 12 0.282002 0.122723 13 0.38326 0.916528 14 0.302535 0.35141 15 0.723881 0.237409 16 0.685904 0.595619 17 0.893656 0.501826 18 0.906395 0.69986 19 0.279633 0.767809 20 0.503557 0.280666 21 0.48702 0.494807 22 0.688613 0.410508 23 0.297071 0.571729 24 0.1026 0.287418 25 0.616686 0.916086 After Latin step: Row: 1 2 Col 1 0.18 0.86 2 0.74 0.82 3 0.62 0.02 4 0.42 0.06 5 0.02 0.66 6 0.1 0.46 7 0.5 0.74 8 0.9 0.14 9 0.98 0.34 10 0.06 0.1 11 0.82 0.9 12 0.26 0.18 13 0.38 0.98 14 0.34 0.38 15 0.78 0.22 16 0.66 0.62 17 0.86 0.54 18 0.94 0.7 19 0.22 0.78 20 0.54 0.26 21 0.46 0.5 22 0.7 0.42 23 0.3 0.58 24 0.14 0.3 25 0.58 0.94 STEP L2 Change 1 0.338451 2 0.111042 3 0.071632 4 0.0626951 5 0.0626966 6 0.0477467 7 0.0353726 8 0.0279948 9 0.023427 10 0.0197014 11 0.0167637 12 0.0145824 13 0.0127907 14 0.0115074 15 0.0102967 16 0.00947034 17 0.0088811 18 0.00828136 19 0.00805758 20 0.0075472 21 0.00709667 22 0.00673319 23 0.00649308 24 0.005991 25 0.00564957 26 0.00530154 27 0.00504769 28 0.00470819 29 0.00455862 30 0.0046571 31 0.00427569 32 0.00422415 33 0.00394813 34 0.00369806 35 0.00340872 36 0.00336278 37 0.0033397 38 0.00322033 39 0.00315892 40 0.00296059 41 0.00300028 42 0.00295037 43 0.00313776 44 0.00317702 45 0.00297438 46 0.00314692 47 0.00337754 48 0.00327974 49 0.0036463 50 0.00324877 After CVT steps: Row: 1 2 Col 1 0.122179 0.70241 2 0.0931445 0.899895 3 0.89346 0.0796305 4 0.857334 0.744741 5 0.893268 0.244376 6 0.896464 0.417668 7 0.0939085 0.49481 8 0.104933 0.292763 9 0.287457 0.893535 10 0.0952161 0.0959487 11 0.287731 0.114034 12 0.481452 0.0951924 13 0.307638 0.334347 14 0.495084 0.899305 15 0.696365 0.323631 16 0.616235 0.707412 17 0.38551 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.18 0.66 2 0.02 0.9 3 0.86 0.02 4 0.78 0.78 5 0.82 0.22 6 0.9 0.42 7 0.06 0.46 8 0.14 0.3 9 0.26 0.82 10 0.1 0.1 11 0.3 0.18 12 0.42 0.06 13 0.34 0.38 14 0.46 0.86 15 0.66 0.34 16 0.58 0.7 17 0.38 0.74 18 0.94 0.98 19 0.62 0.14 20 0.22 0.58 21 0.98 0.62 22 0.54 0.26 23 0.7 0.94 24 0.74 0.54 25 0.5 0.5 STEP L2 Change 1 0.2751 2 0.125962 3 0.0844898 4 0.0613687 5 0.0496751 6 0.0436553 7 0.038759 8 0.0356382 9 0.0325845 10 0.0302966 11 0.0278321 12 0.0258386 13 0.0237712 14 0.0227156 15 0.0209787 16 0.0193704 17 0.0186055 18 0.0172389 19 0.0175257 20 0.0172194 21 0.0167848 22 0.0159576 23 0.0150214 24 0.0142012 25 0.0128603 26 0.0114645 27 0.0105762 28 0.00956813 29 0.00894074 30 0.00836945 31 0.00791028 32 0.00786404 33 0.00720992 34 0.00709843 35 0.00687498 36 0.00657366 37 0.0060662 38 0.00592004 39 0.00539722 40 0.0054816 41 0.00518 42 0.00475212 43 0.00444837 44 0.00425722 45 0.00377342 46 0.00373325 47 0.00316305 48 0.00331493 49 0.00347011 50 0.00297985 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.24526 5 0.418435 0.903613 6 0.925064 0.85626 7 0.106918 0.0954984 8 0.234344 0.597594 9 0.559145 0.731769 10 0.77231 0.907061 11 0.880381 0.0989345 12 0.0834157 0.904869 13 0.529296 0.53965 14 0.914645 0.594343 15 0.132777 0.280827 16 0.0910098 0.471195 17 0.74598 0.694047 18 0.334005 0.084003 19 0.372356 0.692983 20 0.0913593 0.701655 21 0.247246 0.875084 22 0.589981 0.093423 23 0.912363 0.343108 24 0.733458 0.239361 25 0.554582 0.335274 After Latin step: Row: 1 2 Col 1 0.66 0.98 2 0.3 0.38 3 0.74 0.42 4 0.42 0.22 5 0.46 0.86 6 0.98 0.78 7 0.14 0.1 8 0.22 0.58 9 0.58 0.74 10 0.82 0.94 11 0.86 0.14 12 0.02 0.9 13 0.5 0.5 14 0.94 0.54 15 0.18 0.26 16 0.06 0.46 17 0.78 0.66 18 0.34 0.02 19 0.38 0.62 20 0.1 0.7 21 0.26 0.82 22 0.62 0.06 23 0.9 0.34 24 0.7 0.18 25 0.54 0.3 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 M = 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.31296 2 0.104298 3 0.0556211 4 0.0405943 5 0.0349141 6 0.032824 7 0.032171 8 0.0319904 9 0.0307178 10 0.0282423 11 0.0261171 12 0.0238931 13 0.0216498 14 0.018758 15 0.0161877 16 0.0139986 17 0.0121291 18 0.0107462 19 0.00972515 20 0.00882685 21 0.00807502 22 0.00749861 23 0.00696677 24 0.00665431 25 0.00648741 26 0.00632216 27 0.00636079 28 0.00650193 29 0.00623285 30 0.00634604 31 0.00602058 32 0.00571068 33 0.00519151 34 0.00564244 35 0.00484621 36 0.00439562 37 0.00426498 38 0.00431384 39 0.00390123 40 0.00352097 41 0.00307958 42 0.00361401 43 0.00330628 44 0.00270564 45 0.00306803 46 0.00259868 47 0.00220686 48 0.00282577 49 0.00248577 50 0.00228596 After CVT steps: Row: 1 2 Col 1 0.117378 0.907721 2 0.716593 0.774318 3 0.664889 0.0824716 4 0.464305 0.0915203 5 0.0923139 0.696472 6 0.102332 0.490041 7 0.492804 0.717175 8 0.891355 0.0966168 9 0.9044 0.304509 10 0.0945985 0.0934519 11 0.87769 0.908594 12 0.281417 0.122941 13 0.381538 0.916179 14 0.301524 0.351511 15 0.721275 0.23728 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.18 0.86 2 0.74 0.82 3 0.62 0.02 4 0.42 0.06 5 0.02 0.66 6 0.14 0.46 7 0.5 0.74 8 0.9 0.14 9 0.98 0.34 10 0.06 0.1 11 0.82 0.9 12 0.26 0.18 13 0.38 0.94 14 0.34 0.38 15 0.78 0.22 16 0.66 0.62 17 0.86 0.54 18 0.94 0.7 19 0.22 0.78 20 0.54 0.26 21 0.46 0.5 22 0.7 0.42 23 0.3 0.58 24 0.1 0.3 25 0.58 0.98 STEP L2 Change 1 0.337894 2 0.111462 3 0.0725916 4 0.0629352 5 0.0632318 6 0.0475428 7 0.0348775 8 0.0280028 9 0.0231014 10 0.0198489 11 0.0171354 12 0.0143585 13 0.0127814 14 0.0113104 15 0.0108035 16 0.00966506 17 0.009454 18 0.00858987 19 0.00814846 20 0.00762192 21 0.00708333 22 0.00658943 23 0.00632145 24 0.00590059 25 0.00580139 26 0.0054616 27 0.00553953 28 0.00463159 29 0.00470327 30 0.0041975 31 0.00494364 32 0.00353426 33 0.00327026 34 0.00338754 35 0.00335437 36 0.00287206 37 0.00266701 38 0.00256649 39 0.00242225 40 0.00234019 41 0.00215599 42 0.00186042 43 0.00148245 44 0.00145194 45 0.00126241 46 0.00121873 47 0.00131654 48 0.00131143 49 0.00130953 50 0.00121048 After CVT steps: Row: 1 2 Col 1 0.122075 0.703468 2 0.0930657 0.900179 3 0.888235 0.0800818 4 0.862021 0.749471 5 0.897071 0.244681 6 0.894083 0.419968 7 0.0920191 0.495608 8 0.106279 0.294782 9 0.287415 0.894629 10 0.093435 0.0971009 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.18 0.7 2 0.06 0.94 3 0.82 0.02 4 0.78 0.78 5 0.94 0.22 6 0.9 0.42 7 0.02 0.46 8 0.14 0.26 9 0.3 0.82 10 0.1 0.06 11 0.26 0.18 12 0.42 0.1 13 0.34 0.38 14 0.5 0.86 15 0.74 0.34 16 0.58 0.66 17 0.38 0.74 18 0.98 0.98 19 0.62 0.14 20 0.22 0.58 21 0.86 0.62 22 0.54 0.3 23 0.7 0.9 24 0.66 0.5 25 0.46 0.54 STEP L2 Change 1 0.274582 2 0.125347 3 0.0843592 4 0.0614323 5 0.0500101 6 0.043372 7 0.0388543 8 0.0353305 9 0.0326621 10 0.0299066 11 0.0277784 12 0.0263485 13 0.0245741 14 0.0227222 15 0.0209569 16 0.0193682 17 0.0183288 18 0.0175441 19 0.017549 20 0.0175268 21 0.0168927 22 0.0159968 23 0.0148146 24 0.0136479 25 0.0123142 26 0.0109877 27 0.0102251 28 0.00952287 29 0.00879783 30 0.00796207 31 0.00765112 32 0.00735978 33 0.00695192 34 0.00687419 35 0.00682489 36 0.00647327 37 0.00614859 38 0.00603553 39 0.00561091 40 0.00535425 41 0.00484898 42 0.00486607 43 0.00455513 44 0.00435351 45 0.00408314 46 0.00397443 47 0.00393231 48 0.00374559 49 0.00361385 50 0.00364924 After CVT steps: Row: 1 2 Col 1 0.596934 0.907657 2 0.329658 0.425128 3 0.72961 0.472162 4 0.38134 0.241921 5 0.4171 0.90391 6 0.919375 0.856502 7 0.106154 0.0953542 8 0.238339 0.583877 9 0.55254 0.729923 10 0.767377 0.907234 11 0.878171 0.0994474 12 0.0822801 0.903837 13 0.520696 0.540445 14 0.90878 0.595419 15 0.132116 0.280078 16 0.0891819 0.470711 17 0.740222 0.694383 18 0.33174 0.0836677 19 0.360906 0.699381 20 0.0956686 0.70097 21 0.245724 0.880254 22 0.589355 0.0930419 23 0.908868 0.344603 24 0.731338 0.240013 25 0.553146 0.335697 After Latin step: Row: 1 2 Col 1 0.66 0.98 2 0.3 0.38 3 0.7 0.46 4 0.42 0.22 5 0.46 0.9 6 0.98 0.78 7 0.14 0.1 8 0.22 0.54 9 0.54 0.74 10 0.82 0.94 11 0.86 0.14 12 0.02 0.86 13 0.5 0.5 14 0.9 0.58 15 0.18 0.26 16 0.06 0.42 17 0.78 0.62 18 0.34 0.02 19 0.38 0.66 20 0.1 0.7 21 0.26 0.82 22 0.62 0.06 23 0.94 0.34 24 0.74 0.18 25 0.58 0.3 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 M = 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.1 0.1 2 0.1 0.3 3 0.1 0.5 4 0.1 0.7 5 0.1 0.9 6 0.3 0.1 7 0.3 0.3 8 0.3 0.5 9 0.3 0.7 10 0.3 0.9 11 0.5 0.1 12 0.5 0.3 13 0.5 0.5 14 0.5 0.7 15 0.5 0.9 16 0.7 0.1 17 0.7 0.3 18 0.7 0.5 19 0.7 0.7 20 0.7 0.9 21 0.9 0.1 22 0.9 0.3 23 0.9 0.5 24 0.9 0.7 25 0.9 0.9 STEP L2 Change 1 0.00664017 2 0.00939851 3 0.00800209 4 0.00883319 5 0.00783317 6 0.00696342 7 0.00756035 8 0.00796526 9 0.0084684 10 0.00645734 11 0.0072929 12 0.00682437 13 0.00856602 14 0.00781099 15 0.00805214 16 0.00999112 17 0.00855996 18 0.00821275 19 0.00921709 20 0.007873 21 0.00756445 22 0.00806682 23 0.00977954 24 0.00868791 25 0.00893513 26 0.00901741 27 0.00954189 28 0.00928846 29 0.00920926 30 0.00795912 31 0.00831646 32 0.00780843 33 0.00943376 34 0.00919017 35 0.00858552 36 0.00838388 37 0.00831205 38 0.00879097 39 0.00865118 40 0.00765336 41 0.00811211 42 0.0066943 43 0.00679794 44 0.00745861 45 0.00724754 46 0.00661947 47 0.00763879 48 0.00785157 49 0.00654922 50 0.00771335 After CVT steps: Row: 1 2 Col 1 0.10225 0.0916803 2 0.101744 0.27988 3 0.100881 0.477441 4 0.0994338 0.684091 5 0.099091 0.893837 6 0.29912 0.112164 7 0.300636 0.330854 8 0.300272 0.542215 9 0.30368 0.738593 10 0.302822 0.916638 11 0.499432 0.0821377 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.69731 0.537542 19 0.694169 0.73581 20 0.703735 0.916048 21 0.900549 0.0929568 22 0.896097 0.281587 23 0.89817 0.480598 24 0.899424 0.684301 25 0.901173 0.89104 After Latin step: Row: 1 2 Col 1 0.18 0.06 2 0.14 0.26 3 0.1 0.46 4 0.06 0.66 5 0.02 0.9 6 0.22 0.18 7 0.3 0.38 8 0.26 0.58 9 0.38 0.78 10 0.34 0.98 11 0.46 0.02 12 0.58 0.22 13 0.42 0.42 14 0.5 0.62 15 0.54 0.82 16 0.74 0.14 17 0.7 0.34 18 0.66 0.54 19 0.62 0.74 20 0.78 0.94 21 0.94 0.1 22 0.82 0.3 23 0.86 0.5 24 0.9 0.7 25 0.98 0.86 STEP L2 Change 1 0.179457 2 0.0486928 3 0.0271824 4 0.0220595 5 0.0193753 6 0.0191377 7 0.0177937 8 0.0165359 9 0.0150058 10 0.0136154 11 0.0110715 12 0.0113829 13 0.0108745 14 0.00962225 15 0.010147 16 0.00850323 17 0.0084107 18 0.00713258 19 0.00780919 20 0.00827701 21 0.0074581 22 0.00798757 23 0.0081771 24 0.00785564 25 0.00903792 26 0.00799633 27 0.00733387 28 0.00835392 29 0.00860054 30 0.00802456 31 0.00778743 32 0.00687197 33 0.00755769 34 0.00846418 35 0.0077681 36 0.00810695 37 0.00802014 38 0.00733861 39 0.00674251 40 0.0074609 41 0.00745142 42 0.00699838 43 0.00782354 44 0.00727416 45 0.00762762 46 0.00707362 47 0.00823569 48 0.00717548 49 0.00741275 50 0.00899789 After CVT steps: Row: 1 2 Col 1 0.0999069 0.0909814 2 0.10328 0.277088 3 0.102898 0.476563 4 0.100763 0.685865 5 0.0981655 0.89357 6 0.297349 0.113582 7 0.298684 0.336336 8 0.29817 0.546245 9 0.302323 0.741549 10 0.302317 0.918583 11 0.495389 0.0842624 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.0924848 22 0.897278 0.27843 23 0.898221 0.479599 24 0.899695 0.683245 25 0.899564 0.891528 After Latin step: Row: 1 2 Col 1 0.06 0.06 2 0.18 0.26 3 0.14 0.46 4 0.1 0.7 5 0.02 0.9 6 0.22 0.18 7 0.3 0.38 8 0.26 0.58 9 0.38 0.78 10 0.34 0.98 11 0.42 0.02 12 0.5 0.22 13 0.46 0.42 14 0.54 0.62 15 0.58 0.82 16 0.7 0.14 17 0.66 0.34 18 0.62 0.54 19 0.74 0.74 20 0.78 0.94 21 0.9 0.1 22 0.82 0.3 23 0.86 0.5 24 0.98 0.66 25 0.94 0.86 STEP L2 Change 1 0.155244 2 0.0540019 3 0.0330229 4 0.0243982 5 0.0225713 6 0.0169874 7 0.0156163 8 0.0129992 9 0.0123533 10 0.0110597 11 0.0106954 12 0.00909802 13 0.00869555 14 0.00823592 15 0.00790784 16 0.00799396 17 0.0086346 18 0.0085257 19 0.00876549 20 0.00876285 21 0.00784343 22 0.00713787 23 0.00754707 24 0.00770984 25 0.00804082 26 0.00916429 27 0.00696382 28 0.00754476 29 0.00836523 30 0.00749182 31 0.00747464 32 0.00844097 33 0.00786318 34 0.0069687 35 0.00808985 36 0.0092441 37 0.00821491 38 0.00704352 39 0.00805498 40 0.00683954 41 0.00726091 42 0.00587844 43 0.00755522 44 0.00749553 45 0.00737641 46 0.00739666 47 0.00789929 48 0.0067235 49 0.00787061 50 0.00658525 After CVT steps: Row: 1 2 Col 1 0.0996875 0.0920661 2 0.104555 0.280829 3 0.100363 0.478226 4 0.100379 0.682634 5 0.0999222 0.893572 6 0.299359 0.112905 7 0.302327 0.33392 8 0.299105 0.543555 9 0.303829 0.739018 10 0.301184 0.916063 11 0.500746 0.0817582 12 0.502075 0.261845 13 0.497662 0.460535 14 0.501193 0.666995 15 0.50106 0.888249 16 0.700532 0.112159 17 0.695906 0.335459 18 0.696993 0.540627 19 0.698546 0.738937 20 0.70021 0.915012 21 0.898977 0.0905614 22 0.89529 0.277769 23 0.898986 0.475041 24 0.90033 0.679836 25 0.900599 0.894167 After Latin step: Row: 1 2 Col 1 0.02 0.1 2 0.18 0.3 3 0.1 0.5 4 0.14 0.7 5 0.06 0.86 6 0.26 0.18 7 0.34 0.34 8 0.22 0.58 9 0.38 0.78 10 0.3 0.98 11 0.46 0.02 12 0.58 0.22 13 0.42 0.42 14 0.54 0.62 15 0.5 0.82 16 0.78 0.14 17 0.62 0.38 18 0.66 0.54 19 0.7 0.74 20 0.74 0.94 21 0.86 0.06 22 0.82 0.26 23 0.9 0.46 24 0.94 0.66 25 0.98 0.9 LCVT_PRB Normal end of execution. 22 December 2011 04:00:05 PM