September 8 2010 5:28:49.841 PM CVT_PRB FORTRAN90 version Test the CVT library. TEST01 CVT computes a Centroidal Voronoi Tessellation. Step SEED L2-Change Energy 0 123456789 0.341992E-01 1 396524255 0.881218 0.339084E-01 2 1804077568 0.311239 0.232023E-01 3 58279807 0.191678 0.204726E-01 4 890554932 0.151504 0.195617E-01 5 214536894 0.106324 0.189207E-01 6 1592773222 0.117191 0.189294E-01 7 1874550474 0.957184E-01 0.184396E-01 8 2041615702 0.977054E-01 0.183635E-01 9 2056573688 0.910568E-01 0.182822E-01 10 1522461057 0.839276E-01 0.177141E-01 11 1338654111 0.821802E-01 0.178577E-01 12 758290347 0.918929E-01 0.176410E-01 13 1895558656 0.753929E-01 0.176485E-01 14 553115540 0.706886E-01 0.171265E-01 15 1150573110 0.665998E-01 0.174222E-01 16 1130930562 0.696999E-01 0.173089E-01 17 1398936087 0.706177E-01 0.172119E-01 18 1351174075 0.581978E-01 0.170374E-01 19 1271270653 0.480012E-01 0.173369E-01 20 383674773 0.556872E-01 0.172128E-01 21 533903734 0.486140E-01 0.173451E-01 22 412527942 0.574532E-01 0.169264E-01 23 498761874 0.554653E-01 0.170907E-01 24 1767322623 0.594193E-01 0.172429E-01 25 8802113 0.421582E-01 0.171463E-01 26 337656644 0.357121E-01 0.169909E-01 27 982491412 0.473630E-01 0.170611E-01 28 1661195312 0.410137E-01 0.171046E-01 29 1727223671 0.378915E-01 0.170059E-01 30 1934076831 0.424840E-01 0.171542E-01 31 286968437 0.555760E-01 0.170516E-01 32 1651892980 0.339753E-01 0.173517E-01 33 1093660108 0.453935E-01 0.172445E-01 34 1312296485 0.495113E-01 0.170867E-01 35 1505118054 0.449818E-01 0.172817E-01 36 1225827697 0.469830E-01 0.173189E-01 37 2103031061 0.460460E-01 0.170275E-01 38 237994259 0.451147E-01 0.170597E-01 39 308344071 0.470322E-01 0.171993E-01 40 878651545 0.402539E-01 0.169572E-01 Dimension DIM_NUM = 2 Number of points N = 10 Initial SEED = 123456789 Current SEED = 336554041 INIT = "uniform". Max iterations IT_MAX = 40 IT_FIXED (fixed samples) 1 Iterations IT_NUM = 40 Difference IT_DIFF = 0.402539E-01 CVT ENERGY = 0.169572E-01 SAMPLE = "uniform". Samples SAMPLE_NUM = 10000 Sampling BATCH size = 1000 EPSILON (unit roundoff) = 0.222045E-15 Generators (rows): Row 1 2 Col 1 0.173177 0.883037 2 0.835817 0.827856 3 0.517580 0.498432 4 0.509540 0.165024 5 0.182392 0.622352 6 0.179416 0.374999 7 0.513340 0.833788 8 0.834205 0.165368 9 0.834444 0.494964 10 0.173130 0.122778 TEST02 CVT computes a Centroidal Voronoi Tessellation. Repeat test 1, but with twice the number of iterations. Step SEED L2-Change Energy 0 123456789 0.341992E-01 1 396524255 0.881218 0.339084E-01 2 1804077568 0.311239 0.232023E-01 3 58279807 0.191678 0.204726E-01 4 890554932 0.151504 0.195617E-01 5 214536894 0.106324 0.189207E-01 6 1592773222 0.117191 0.189294E-01 7 1874550474 0.957184E-01 0.184396E-01 8 2041615702 0.977054E-01 0.183635E-01 9 2056573688 0.910568E-01 0.182822E-01 10 1522461057 0.839276E-01 0.177141E-01 11 1338654111 0.821802E-01 0.178577E-01 12 758290347 0.918929E-01 0.176410E-01 13 1895558656 0.753929E-01 0.176485E-01 14 553115540 0.706886E-01 0.171265E-01 15 1150573110 0.665998E-01 0.174222E-01 16 1130930562 0.696999E-01 0.173089E-01 17 1398936087 0.706177E-01 0.172119E-01 18 1351174075 0.581978E-01 0.170374E-01 19 1271270653 0.480012E-01 0.173369E-01 20 383674773 0.556872E-01 0.172128E-01 21 533903734 0.486140E-01 0.173451E-01 22 412527942 0.574532E-01 0.169264E-01 23 498761874 0.554653E-01 0.170907E-01 24 1767322623 0.594193E-01 0.172429E-01 25 8802113 0.421582E-01 0.171463E-01 26 337656644 0.357121E-01 0.169909E-01 27 982491412 0.473630E-01 0.170611E-01 28 1661195312 0.410137E-01 0.171046E-01 29 1727223671 0.378915E-01 0.170059E-01 30 1934076831 0.424840E-01 0.171542E-01 31 286968437 0.555760E-01 0.170516E-01 32 1651892980 0.339753E-01 0.173517E-01 33 1093660108 0.453935E-01 0.172445E-01 34 1312296485 0.495113E-01 0.170867E-01 35 1505118054 0.449818E-01 0.172817E-01 36 1225827697 0.469830E-01 0.173189E-01 37 2103031061 0.460460E-01 0.170275E-01 38 237994259 0.451147E-01 0.170597E-01 39 308344071 0.470322E-01 0.171993E-01 40 878651545 0.402539E-01 0.169572E-01 41 336554041 0.444687E-01 0.169739E-01 42 69693939 0.417197E-01 0.172409E-01 43 107597223 0.444121E-01 0.170524E-01 44 284388055 0.449483E-01 0.170296E-01 45 821532605 0.336557E-01 0.170280E-01 46 1912444717 0.438442E-01 0.168697E-01 47 2139238658 0.640732E-01 0.171048E-01 48 1514249906 0.440515E-01 0.172451E-01 49 765384573 0.472810E-01 0.171670E-01 50 599511796 0.477820E-01 0.170414E-01 51 1929718247 0.526363E-01 0.169852E-01 52 1495267469 0.577634E-01 0.171637E-01 53 346287018 0.357526E-01 0.171425E-01 54 972591208 0.403661E-01 0.172913E-01 55 1170908305 0.345025E-01 0.169270E-01 56 2101087755 0.471512E-01 0.169423E-01 57 1585775448 0.427815E-01 0.169886E-01 58 1719007417 0.491513E-01 0.171720E-01 59 47014704 0.341510E-01 0.169112E-01 60 325311895 0.488737E-01 0.172634E-01 61 2114713435 0.425484E-01 0.169471E-01 62 469955390 0.426950E-01 0.172777E-01 63 2044044375 0.426745E-01 0.169180E-01 64 1902665122 0.429220E-01 0.171021E-01 65 999386148 0.488367E-01 0.170877E-01 66 1762727540 0.441203E-01 0.168713E-01 67 886787757 0.413277E-01 0.168756E-01 68 1442441010 0.446965E-01 0.170550E-01 69 766630007 0.452238E-01 0.169835E-01 70 473489672 0.286167E-01 0.170666E-01 71 1247561645 0.220248E-01 0.172088E-01 72 1093781229 0.437517E-01 0.171667E-01 73 1723735430 0.416874E-01 0.169444E-01 74 1100723884 0.540078E-01 0.170263E-01 75 2145312587 0.442257E-01 0.169726E-01 76 1556984304 0.393601E-01 0.170680E-01 77 1035291506 0.376604E-01 0.168979E-01 78 158213521 0.434274E-01 0.171770E-01 79 1863192335 0.410352E-01 0.170552E-01 80 343175222 0.408352E-01 0.170907E-01 Dimension DIM_NUM = 2 Number of points N = 10 Initial SEED = 123456789 Current SEED = 834195697 INIT = "uniform". Max iterations IT_MAX = 80 IT_FIXED (fixed samples) 1 Iterations IT_NUM = 80 Difference IT_DIFF = 0.408352E-01 CVT ENERGY = 0.170907E-01 SAMPLE = "uniform". Samples SAMPLE_NUM = 10000 Sampling BATCH size = 1000 EPSILON (unit roundoff) = 0.222045E-15 Generators (rows): Row 1 2 Col 1 0.171691 0.880159 2 0.836688 0.823236 3 0.517862 0.508216 4 0.509477 0.172695 5 0.178214 0.631626 6 0.183649 0.374228 7 0.509519 0.845054 8 0.839871 0.164440 9 0.838305 0.493313 10 0.168388 0.125100 TEST03 CVT computes a Centroidal Voronoi Tessellation. Repeat test 1, but with 100 times the sample points. Step SEED L2-Change Energy 0 123456789 0.342772E-01 1 720688506 0.893257 0.343633E-01 2 681172 0.290172 0.228573E-01 3 539098006 0.196318 0.205969E-01 4 1544551628 0.147817 0.196212E-01 5 1039137261 0.116655 0.190707E-01 6 1856120581 0.946072E-01 0.187394E-01 7 1398888834 0.854150E-01 0.184761E-01 8 777275849 0.820861E-01 0.182815E-01 9 193908656 0.750730E-01 0.180899E-01 10 594766309 0.669630E-01 0.179151E-01 11 1528604881 0.614130E-01 0.177392E-01 12 2135424592 0.530776E-01 0.176312E-01 13 843278689 0.496226E-01 0.175196E-01 14 1298034490 0.427557E-01 0.174754E-01 15 1774985135 0.412219E-01 0.174116E-01 16 58892340 0.392985E-01 0.173803E-01 17 497310354 0.364582E-01 0.173483E-01 18 1555340604 0.359526E-01 0.172887E-01 19 573219264 0.358032E-01 0.172589E-01 20 769258892 0.354920E-01 0.172147E-01 21 1733894875 0.341897E-01 0.171860E-01 22 1708504563 0.342741E-01 0.171507E-01 23 332981687 0.275611E-01 0.171264E-01 24 871833470 0.237202E-01 0.171083E-01 25 491107991 0.209864E-01 0.170851E-01 26 64836647 0.170961E-01 0.170871E-01 27 954805664 0.151398E-01 0.170704E-01 28 1033236242 0.130872E-01 0.170698E-01 29 1875753496 0.109179E-01 0.170671E-01 30 932903206 0.868216E-02 0.170777E-01 31 900419333 0.877443E-02 0.170689E-01 32 898406847 0.765317E-02 0.170628E-01 33 957321340 0.601815E-02 0.170555E-01 34 1879036371 0.610091E-02 0.170541E-01 35 550745514 0.600529E-02 0.170619E-01 36 1822412817 0.437686E-02 0.170696E-01 37 1926341691 0.460332E-02 0.170756E-01 38 578466603 0.462387E-02 0.170466E-01 39 1767869167 0.476128E-02 0.170515E-01 40 1282315487 0.489911E-02 0.170678E-01 Dimension DIM_NUM = 2 Number of points N = 10 Initial SEED = 123456789 Current SEED = 1300455975 INIT = "uniform". Max iterations IT_MAX = 40 IT_FIXED (fixed samples) 1 Iterations IT_NUM = 40 Difference IT_DIFF = 0.489911E-02 CVT ENERGY = 0.170678E-01 SAMPLE = "uniform". Samples SAMPLE_NUM = 1000000 Sampling BATCH size = 1000 EPSILON (unit roundoff) = 0.222045E-15 Generators (rows): Row 1 2 Col 1 0.167765 0.878515 2 0.834745 0.834827 3 0.522844 0.496098 4 0.504533 0.163541 5 0.183085 0.628051 6 0.181541 0.374185 7 0.507580 0.832903 8 0.834187 0.167704 9 0.842937 0.503352 10 0.166800 0.122784 TEST04 CVT computes a Centroidal Voronoi Tessellation. Repeat test 1, but with Halton sampling. Step SEED L2-Change Energy 0 123456789 0.343097E-01 1 29252052 0.890558 0.343249E-01 2 29262052 0.294653 0.228400E-01 3 29272052 0.196791 0.206273E-01 4 29282052 0.149052 0.196306E-01 5 29292052 0.117047 0.190654E-01 6 29302052 0.975388E-01 0.187154E-01 7 29312052 0.873418E-01 0.184721E-01 8 29322052 0.812065E-01 0.182544E-01 9 29332052 0.754804E-01 0.180372E-01 10 29342052 0.697301E-01 0.178618E-01 11 29352052 0.624813E-01 0.177162E-01 12 29362052 0.543370E-01 0.175969E-01 13 29372052 0.487347E-01 0.175152E-01 14 29382052 0.455322E-01 0.174390E-01 15 29392052 0.420863E-01 0.173858E-01 16 29402052 0.378334E-01 0.173413E-01 17 29412052 0.367790E-01 0.173118E-01 18 29422052 0.385578E-01 0.172580E-01 19 29432052 0.364872E-01 0.172285E-01 20 29442052 0.343201E-01 0.171788E-01 21 29452052 0.344728E-01 0.171618E-01 22 29462052 0.286870E-01 0.171296E-01 23 29472052 0.225448E-01 0.171019E-01 24 29482052 0.212899E-01 0.171036E-01 25 29492052 0.178806E-01 0.170790E-01 26 29502052 0.156606E-01 0.170798E-01 27 29512052 0.138195E-01 0.170728E-01 28 29522052 0.101995E-01 0.170707E-01 29 29532052 0.954090E-02 0.170745E-01 30 29542052 0.102057E-01 0.170601E-01 31 29552052 0.887253E-02 0.170691E-01 32 29562052 0.844319E-02 0.170647E-01 33 29572052 0.422509E-02 0.170650E-01 34 29582052 0.710075E-02 0.170591E-01 35 29592052 0.710168E-02 0.170693E-01 36 29602052 0.570035E-02 0.170602E-01 37 29612052 0.763221E-02 0.170606E-01 38 29622052 0.658065E-02 0.170677E-01 39 29632052 0.664932E-02 0.170602E-01 40 29642052 0.729668E-02 0.170655E-01 Dimension DIM_NUM = 2 Number of points N = 10 Initial SEED = 123456789 Current SEED = 29652052 INIT = "uniform". Max iterations IT_MAX = 40 IT_FIXED (fixed samples) 1 Iterations IT_NUM = 40 Difference IT_DIFF = 0.729668E-02 CVT ENERGY = 0.170655E-01 SAMPLE = "halton". Samples SAMPLE_NUM = 10000 Sampling BATCH size = 1000 EPSILON (unit roundoff) = 0.222045E-15 Generators (rows): Row 1 2 Col 1 0.168095 0.878286 2 0.834748 0.835029 3 0.522741 0.497604 4 0.504804 0.163962 5 0.182069 0.627496 6 0.181252 0.372989 7 0.505663 0.832959 8 0.834007 0.168284 9 0.842571 0.503020 10 0.166927 0.122462 TEST05 CVT computes a Centroidal Voronoi Tessellation. Repeat test 1, but with grid sampling. Step SEED L2-Change Energy 0 123456789 0.343146E-01 1 29252052 0.894182 0.343146E-01 2 29262052 0.289683 0.227960E-01 3 29272052 0.192593 0.206112E-01 4 29282052 0.151542 0.196286E-01 5 29292052 0.113430 0.190577E-01 6 29302052 0.953549E-01 0.187215E-01 7 29312052 0.883581E-01 0.184768E-01 8 29322052 0.802806E-01 0.182516E-01 9 29332052 0.738208E-01 0.180469E-01 10 29342052 0.692118E-01 0.178695E-01 11 29352052 0.612144E-01 0.177196E-01 12 29362052 0.541458E-01 0.176043E-01 13 29372052 0.493226E-01 0.175162E-01 14 29382052 0.451789E-01 0.174462E-01 15 29392052 0.408350E-01 0.173893E-01 16 29402052 0.400852E-01 0.173437E-01 17 29412052 0.370572E-01 0.173011E-01 18 29422052 0.393376E-01 0.172607E-01 19 29432052 0.335672E-01 0.172201E-01 20 29442052 0.346096E-01 0.171853E-01 21 29452052 0.329736E-01 0.171526E-01 22 29462052 0.299827E-01 0.171248E-01 23 29472052 0.255714E-01 0.171040E-01 24 29482052 0.192152E-01 0.170893E-01 25 29492052 0.166202E-01 0.170801E-01 26 29502052 0.163850E-01 0.170733E-01 27 29512052 0.104304E-01 0.170680E-01 28 29522052 0.886042E-02 0.170659E-01 29 29532052 0.779801E-02 0.170643E-01 30 29542052 0.971012E-02 0.170628E-01 31 29552052 0.972930E-02 0.170605E-01 32 29562052 0.788473E-02 0.170584E-01 33 29572052 0.601319E-02 0.170569E-01 34 29582052 0.479507E-02 0.170562E-01 35 29592052 0.293378E-02 0.170556E-01 36 29602052 0.126444E-02 0.170554E-01 37 29612052 0.193406E-02 0.170553E-01 38 29622052 0.182779E-02 0.170552E-01 39 29632052 0.148689E-02 0.170551E-01 40 29642052 0.188772E-02 0.170550E-01 Dimension DIM_NUM = 2 Number of points N = 10 Initial SEED = 123456789 Current SEED = 29652052 INIT = "uniform". Max iterations IT_MAX = 40 IT_FIXED (fixed samples) 1 Iterations IT_NUM = 40 Difference IT_DIFF = 0.188772E-02 CVT ENERGY = 0.170550E-01 SAMPLE = "grid". Samples SAMPLE_NUM = 10000 Sampling BATCH size = 1000 EPSILON (unit roundoff) = 0.222045E-15 Generators (rows): Row 1 2 Col 1 0.168101 0.877639 2 0.835000 0.835000 3 0.521731 0.496450 4 0.505149 0.163885 5 0.182271 0.626967 6 0.181448 0.372936 7 0.506241 0.833254 8 0.834054 0.168216 9 0.841670 0.503786 10 0.167071 0.122429 TEST06 CVT computes a Centroidal Voronoi Tessellation. Repeat test 1, with FORTRAN90 RANDOM_NUMBER sampling. Step SEED L2-Change Energy 0 123456789 0.344954E-01 1 29262052 0.884165 0.338552E-01 2 29282052 0.302118 0.226992E-01 3 29302052 0.222015 0.207643E-01 4 29322052 0.144794 0.196236E-01 5 29342052 0.137524 0.190493E-01 6 29362052 0.873487E-01 0.188341E-01 7 29382052 0.927746E-01 0.185623E-01 8 29402052 0.101195 0.182930E-01 9 29422052 0.105100 0.180857E-01 10 29442052 0.805787E-01 0.178387E-01 11 29462052 0.656002E-01 0.177410E-01 12 29482052 0.779175E-01 0.180097E-01 13 29502052 0.770351E-01 0.173875E-01 14 29522052 0.785944E-01 0.172862E-01 15 29542052 0.814296E-01 0.174516E-01 16 29562052 0.475482E-01 0.174093E-01 17 29582052 0.595356E-01 0.170314E-01 18 29602052 0.610694E-01 0.172289E-01 19 29622052 0.558575E-01 0.171896E-01 20 29642052 0.378768E-01 0.172038E-01 21 29662052 0.409548E-01 0.172010E-01 22 29682052 0.414104E-01 0.169993E-01 23 29702052 0.416602E-01 0.170319E-01 24 29722052 0.467810E-01 0.171707E-01 25 29742052 0.513598E-01 0.170172E-01 26 29762052 0.457678E-01 0.172050E-01 27 29782052 0.322662E-01 0.169840E-01 28 29802052 0.410102E-01 0.170331E-01 29 29822052 0.415421E-01 0.170189E-01 30 29842052 0.474144E-01 0.168242E-01 31 29862052 0.582895E-01 0.169290E-01 32 29882052 0.493208E-01 0.171863E-01 33 29902052 0.406251E-01 0.170954E-01 34 29922052 0.478630E-01 0.170457E-01 35 29942052 0.317595E-01 0.170808E-01 36 29962052 0.374083E-01 0.171537E-01 37 29982052 0.566266E-01 0.170216E-01 38 30002052 0.436309E-01 0.170259E-01 39 30022052 0.528204E-01 0.171164E-01 40 30042052 0.277773E-01 0.171476E-01 Dimension DIM_NUM = 2 Number of points N = 10 Initial SEED = 123456789 Current SEED = 30062052 INIT = "uniform". Max iterations IT_MAX = 40 IT_FIXED (fixed samples) 1 Iterations IT_NUM = 40 Difference IT_DIFF = 0.277773E-01 CVT ENERGY = 0.171476E-01 SAMPLE = "RANDOM". Samples SAMPLE_NUM = 10000 Sampling BATCH size = 1000 EPSILON (unit roundoff) = 0.222045E-15 Generators (rows): Row 1 2 Col 1 0.171613 0.877802 2 0.838659 0.839419 3 0.521948 0.485116 4 0.509032 0.161088 5 0.180875 0.631655 6 0.176225 0.375828 7 0.508620 0.823755 8 0.838593 0.174374 9 0.837919 0.514760 10 0.165807 0.118388 TEST07 CVT computes a Centroidal Voronoi Tessellation. Repeat test 1, but with a different seed. Step SEED L2-Change Energy 0 987654321 0.375864E-01 1 1703739660 0.783121 0.386543E-01 2 1956211282 0.347315 0.232336E-01 3 703763026 0.250606 0.204697E-01 4 878628297 0.213945 0.191143E-01 5 1552444705 0.166292 0.185077E-01 6 390357255 0.122423 0.178363E-01 7 689833854 0.873161E-01 0.173483E-01 8 1751610365 0.578908E-01 0.173027E-01 9 852519522 0.745691E-01 0.173766E-01 10 309030697 0.569482E-01 0.170224E-01 11 1316592958 0.572571E-01 0.170560E-01 12 1105896658 0.462554E-01 0.171196E-01 13 1045195875 0.515713E-01 0.173048E-01 14 1897548003 0.577570E-01 0.169453E-01 15 1870465968 0.666145E-01 0.171019E-01 16 789370021 0.538122E-01 0.171234E-01 17 444524362 0.471216E-01 0.172111E-01 18 689805975 0.433941E-01 0.171347E-01 19 1558080092 0.481297E-01 0.171542E-01 20 1381530061 0.496368E-01 0.171866E-01 21 157511124 0.545409E-01 0.171291E-01 22 16143349 0.413219E-01 0.171619E-01 23 1155157008 0.533501E-01 0.171342E-01 24 25788210 0.457358E-01 0.170686E-01 25 1456770040 0.480610E-01 0.170587E-01 26 922599524 0.451128E-01 0.171621E-01 27 910060503 0.298797E-01 0.170532E-01 28 1036849277 0.437739E-01 0.169254E-01 29 1558644233 0.533064E-01 0.170381E-01 30 192063102 0.455004E-01 0.173339E-01 31 61275996 0.513876E-01 0.168389E-01 32 1575049278 0.443013E-01 0.170992E-01 33 339755903 0.388937E-01 0.171641E-01 34 1360647670 0.351732E-01 0.169650E-01 35 1802859403 0.456284E-01 0.169189E-01 36 2066534888 0.496651E-01 0.169849E-01 37 101044121 0.500693E-01 0.169682E-01 38 340723523 0.554669E-01 0.173174E-01 39 1321862545 0.485230E-01 0.171405E-01 40 1047111519 0.370525E-01 0.169869E-01 Dimension DIM_NUM = 2 Number of points N = 10 Initial SEED = 987654321 Current SEED = 403266883 INIT = "uniform". Max iterations IT_MAX = 40 IT_FIXED (fixed samples) 1 Iterations IT_NUM = 40 Difference IT_DIFF = 0.370525E-01 CVT ENERGY = 0.169869E-01 SAMPLE = "uniform". Samples SAMPLE_NUM = 10000 Sampling BATCH size = 1000 EPSILON (unit roundoff) = 0.222045E-15 Generators (rows): Row 1 2 Col 1 0.826932 0.838484 2 0.497395 0.525672 3 0.120213 0.170507 4 0.502184 0.840736 5 0.372854 0.184051 6 0.878276 0.172642 7 0.627165 0.181649 8 0.835662 0.513581 9 0.160801 0.831680 10 0.167605 0.507742 TEST08 CVT computes a Centroidal Voronoi Tessellation. Repeat test 1 with a different batch size. Step SEED L2-Change Energy 0 123456789 0.341992E-01 1 396524255 0.881218 0.339084E-01 2 1804077568 0.311239 0.232023E-01 3 58279807 0.191678 0.204726E-01 4 890554932 0.151504 0.195617E-01 5 214536894 0.106324 0.189207E-01 6 1592773222 0.117191 0.189294E-01 7 1874550474 0.957184E-01 0.184396E-01 8 2041615702 0.977054E-01 0.183635E-01 9 2056573688 0.910568E-01 0.182822E-01 10 1522461057 0.839276E-01 0.177141E-01 11 1338654111 0.821802E-01 0.178577E-01 12 758290347 0.918929E-01 0.176410E-01 13 1895558656 0.753929E-01 0.176485E-01 14 553115540 0.706886E-01 0.171265E-01 15 1150573110 0.665998E-01 0.174222E-01 16 1130930562 0.696999E-01 0.173089E-01 17 1398936087 0.706177E-01 0.172119E-01 18 1351174075 0.581978E-01 0.170374E-01 19 1271270653 0.480012E-01 0.173369E-01 20 383674773 0.556872E-01 0.172128E-01 21 533903734 0.486140E-01 0.173451E-01 22 412527942 0.574532E-01 0.169264E-01 23 498761874 0.554653E-01 0.170907E-01 24 1767322623 0.594193E-01 0.172429E-01 25 8802113 0.421582E-01 0.171463E-01 26 337656644 0.357121E-01 0.169909E-01 27 982491412 0.473630E-01 0.170611E-01 28 1661195312 0.410137E-01 0.171046E-01 29 1727223671 0.378915E-01 0.170059E-01 30 1934076831 0.424840E-01 0.171542E-01 31 286968437 0.555760E-01 0.170516E-01 32 1651892980 0.339753E-01 0.173517E-01 33 1093660108 0.453935E-01 0.172445E-01 34 1312296485 0.495113E-01 0.170867E-01 35 1505118054 0.449818E-01 0.172817E-01 36 1225827697 0.469830E-01 0.173189E-01 37 2103031061 0.460460E-01 0.170275E-01 38 237994259 0.451147E-01 0.170597E-01 39 308344071 0.470322E-01 0.171993E-01 40 878651545 0.402539E-01 0.169572E-01 Dimension DIM_NUM = 2 Number of points N = 10 Initial SEED = 123456789 Current SEED = 336554041 INIT = "uniform". Max iterations IT_MAX = 40 IT_FIXED (fixed samples) 1 Iterations IT_NUM = 40 Difference IT_DIFF = 0.402539E-01 CVT ENERGY = 0.169572E-01 SAMPLE = "uniform". Samples SAMPLE_NUM = 10000 Sampling BATCH size = 5 EPSILON (unit roundoff) = 0.222045E-15 Generators (rows): Row 1 2 Col 1 0.173177 0.883037 2 0.835817 0.827856 3 0.517580 0.498432 4 0.509540 0.165024 5 0.182392 0.622352 6 0.179416 0.374999 7 0.513340 0.833788 8 0.834205 0.165368 9 0.834444 0.494964 10 0.173130 0.122778 TEST09 CVT computes a Centroidal Voronoi Tessellation. Repeat test 1 with a fixed sample set. (IT_FIXED = IT_MAX) Step SEED L2-Change Energy 0 123456789 0.341992E-01 1 396524255 0.881218 0.339084E-01 2 396524255 0.288490 0.228074E-01 3 396524255 0.181606 0.206921E-01 4 396524255 0.143931 0.198189E-01 5 396524255 0.116477 0.192979E-01 6 396524255 0.984893E-01 0.189683E-01 7 396524255 0.907669E-01 0.187162E-01 8 396524255 0.908941E-01 0.184940E-01 9 396524255 0.922005E-01 0.182592E-01 10 396524255 0.859268E-01 0.180474E-01 11 396524255 0.796400E-01 0.178549E-01 12 396524255 0.800619E-01 0.176870E-01 13 396524255 0.651065E-01 0.175387E-01 14 396524255 0.587560E-01 0.174414E-01 15 396524255 0.546665E-01 0.173577E-01 16 396524255 0.487698E-01 0.172830E-01 17 396524255 0.433497E-01 0.172254E-01 18 396524255 0.379348E-01 0.171841E-01 19 396524255 0.321769E-01 0.171514E-01 20 396524255 0.276877E-01 0.171297E-01 21 396524255 0.246237E-01 0.171122E-01 22 396524255 0.238782E-01 0.170986E-01 23 396524255 0.220947E-01 0.170852E-01 24 396524255 0.228963E-01 0.170739E-01 25 396524255 0.191269E-01 0.170615E-01 26 396524255 0.196181E-01 0.170519E-01 27 396524255 0.192246E-01 0.170418E-01 28 396524255 0.177851E-01 0.170333E-01 29 396524255 0.122718E-01 0.170279E-01 30 396524255 0.131588E-01 0.170243E-01 31 396524255 0.118481E-01 0.170205E-01 32 396524255 0.947071E-02 0.170174E-01 33 396524255 0.932686E-02 0.170154E-01 34 396524255 0.109052E-01 0.170129E-01 35 396524255 0.108748E-01 0.170100E-01 36 396524255 0.751463E-02 0.170078E-01 37 396524255 0.794483E-02 0.170064E-01 38 396524255 0.789481E-02 0.170045E-01 39 396524255 0.912418E-02 0.170024E-01 40 396524255 0.101925E-01 0.169997E-01 Dimension DIM_NUM = 2 Number of points N = 10 Initial SEED = 123456789 Current SEED = 1804077568 INIT = "uniform". Max iterations IT_MAX = 40 IT_FIXED (fixed samples) 40 Iterations IT_NUM = 40 Difference IT_DIFF = 0.101925E-01 CVT ENERGY = 0.169997E-01 SAMPLE = "uniform". Samples SAMPLE_NUM = 10000 Sampling BATCH size = 1000 EPSILON (unit roundoff) = 0.222045E-15 Generators (rows): Row 1 2 Col 1 0.172379 0.845471 2 0.835383 0.828787 3 0.637518 0.482235 4 0.499658 0.155444 5 0.110493 0.510304 6 0.368913 0.506616 7 0.506491 0.828206 8 0.828248 0.154796 9 0.882482 0.479897 10 0.172396 0.172063 TEST10 CVT computes a Centroidal Voronoi Tessellation. Compute 100 points in 3D. Step SEED L2-Change Energy 0 123456789 0.190278E-01 1 914652320 5.71800 0.190795E-01 2 1182014520 2.37674 0.141520E-01 3 734877640 1.72475 0.131748E-01 4 1808590614 1.46361 0.128552E-01 5 1376754667 1.36806 0.125019E-01 6 39225810 1.39710 0.124375E-01 7 1835665933 1.35677 0.122988E-01 8 683368702 1.36786 0.121452E-01 9 1192180323 1.21410 0.120655E-01 10 1451564322 1.19831 0.120431E-01 11 1525473369 1.21355 0.120139E-01 12 1258784351 1.23666 0.119742E-01 13 1939780355 1.21704 0.119680E-01 14 1792150340 1.20606 0.119348E-01 15 1077745899 1.20416 0.118143E-01 16 1381981537 1.23797 0.118297E-01 17 2135871724 1.09343 0.118518E-01 18 733709902 1.16174 0.118788E-01 19 1511987296 1.19554 0.118066E-01 20 1849799130 1.08146 0.117540E-01 21 1011187751 1.21392 0.117771E-01 22 430762708 1.19615 0.118007E-01 23 2076182163 1.15623 0.118430E-01 24 662061697 1.20557 0.117898E-01 25 253385434 1.17947 0.117468E-01 26 836175832 1.18853 0.117202E-01 27 529591688 1.17680 0.117259E-01 28 154633765 1.12259 0.117726E-01 29 778723351 1.05924 0.117629E-01 30 1341651177 1.12813 0.117690E-01 31 11115418 1.16819 0.117621E-01 32 181250519 1.20900 0.117427E-01 33 628655915 1.14151 0.117160E-01 34 1574892130 1.10617 0.117580E-01 35 1853258335 1.13261 0.117641E-01 36 1106992749 1.15861 0.116933E-01 37 2113411425 1.14066 0.117347E-01 38 1552081075 1.21592 0.117278E-01 39 921807728 1.14794 0.117371E-01 40 1707197995 1.19541 0.117497E-01 Dimension DIM_NUM = 3 Number of points N = 100 Initial SEED = 123456789 Current SEED = 1382742100 INIT = "uniform". Max iterations IT_MAX = 40 IT_FIXED (fixed samples) 1 Iterations IT_NUM = 40 Difference IT_DIFF = 1.19541 CVT ENERGY = 0.117497E-01 SAMPLE = "uniform". Samples SAMPLE_NUM = 10000 Sampling BATCH size = 1000 EPSILON (unit roundoff) = 0.222045E-15 First 10 generators (rows): Row 1 2 3 Col 1 0.329372 0.902757 0.880838 2 0.442548 0.416484 0.945376E-01 3 0.329721 0.110790 0.768365E-01 4 0.705992 0.107839 0.631626 5 0.471560 0.722953 0.880901 6 0.127223 0.919186 0.370276 7 0.861713E-01 0.101029 0.900459 8 0.901934 0.100079 0.119446 9 0.328626 0.920465 0.116447 10 0.440031 0.908374 0.314557 TEST11 CVT computes a Centroidal Voronoi Tessellation. In this test, we initialize the generators to grid points; this is an unstable CVT solution. Initial generators (rows): Row 1 2 Col 1 0.125000 0.125000 2 0.125000 0.375000 3 0.125000 0.625000 4 0.125000 0.875000 5 0.375000 0.125000 6 0.375000 0.375000 7 0.375000 0.625000 8 0.375000 0.875000 9 0.625000 0.125000 10 0.625000 0.375000 11 0.625000 0.625000 12 0.625000 0.875000 13 0.875000 0.125000 14 0.875000 0.375000 15 0.875000 0.625000 16 0.875000 0.875000 Step SEED L2-Change Energy 0 123456789 0.104645E-01 1 1328042058 0.161104 0.109193E-01 2 1539989115 0.199621 0.105351E-01 3 200382020 0.209617 0.102604E-01 4 1958948863 0.259336 0.101123E-01 5 691823909 0.213672 0.105139E-01 6 655816003 0.235274 0.109313E-01 7 450926852 0.261111 0.108770E-01 8 1273398721 0.250138 0.105427E-01 9 1462772409 0.241644 0.103570E-01 10 1371189933 0.223922 0.110969E-01 11 845524934 0.228212 0.105100E-01 12 2013025619 0.244735 0.107247E-01 13 1021685265 0.239236 0.106548E-01 14 204970068 0.246327 0.106566E-01 15 778034624 0.216993 0.103616E-01 16 1938972586 0.242286 0.107049E-01 17 1519833303 0.273886 0.103971E-01 18 2077270047 0.230644 0.105690E-01 19 1923497586 0.160058 0.104248E-01 20 167698528 0.221340 0.103729E-01 21 1655479938 0.196210 0.103319E-01 22 475051709 0.171397 0.100915E-01 23 83696007 0.193487 0.101762E-01 24 1805594913 0.190997 0.102699E-01 25 508154254 0.213475 0.102875E-01 26 1861070898 0.225397 0.106690E-01 27 1661633678 0.200201 0.104139E-01 28 914834335 0.206418 0.103726E-01 29 715643788 0.185243 0.106074E-01 30 810642687 0.190593 0.104067E-01 31 471503978 0.251515 0.105335E-01 32 1019331795 0.266192 0.108879E-01 33 2118751796 0.215274 0.109702E-01 34 1168375939 0.221011 0.108212E-01 35 26482542 0.207715 0.102305E-01 36 2056065863 0.203282 0.108559E-01 37 414064984 0.253232 0.104830E-01 38 1342855162 0.173817 0.101891E-01 39 444956424 0.198617 0.100251E-01 40 2115448041 0.228793 0.108514E-01 Dimension DIM_NUM = 2 Number of points N = 16 Initial SEED = 123456789 Current SEED = 118254244 INIT = "user initialization". Max iterations IT_MAX = 40 IT_FIXED (fixed samples) 1 Iterations IT_NUM = 40 Difference IT_DIFF = 0.228793 CVT ENERGY = 0.108514E-01 SAMPLE = "uniform". Samples SAMPLE_NUM = 1000 Sampling BATCH size = 1000 EPSILON (unit roundoff) = 0.222045E-15 Final generators (rows): Row 1 2 Col 1 0.121439 0.124202 2 0.133557 0.392903 3 0.118580 0.656708 4 0.129644 0.880618 5 0.362683 0.119017 6 0.385495 0.322734 7 0.361119 0.594092 8 0.377020 0.858389 9 0.607368 0.127296 10 0.643546 0.378122 11 0.600855 0.632838 12 0.647725 0.875136 13 0.884853 0.117756 14 0.877764 0.343945 15 0.872582 0.597904 16 0.882169 0.869201 TEST12 The "random" initialization option calls the system random number generator. There is some question about whether this works correctly. The test is as follows: CVT call #1: DIM_NUM = 2 N = 10 INIT = -1 IT_MAX = 0 SEED = 100000 Print output values of SEED and R #1. CVT call #2: (jump SEED) DIM_NUM = 2 N = 10 INIT = -1 IT_MAX = 0 SEED = 200000. Print output values of SEED and R #2. CVT call #3: (restore SEED) DIM_NUM = 2 N = 10 INIT = -1 IT_MAX = 0 SEED_INIT = 100000 Print output values of SEED and R #3. We expect that: * the values of R #1 and R #2 differ; AND * the values of R #1 and R #3 agree. Step SEED L2-Change Energy 0 100000 0.285261E-01 Dimension DIM_NUM = 2 Number of points N = 10 Initial SEED = 100000 Current SEED = 1647317536 INIT = "random". Max iterations IT_MAX = 0 IT_FIXED (fixed samples) 1 Iterations IT_NUM = 0 Difference IT_DIFF = 0.00000 CVT ENERGY = 0.285261E-01 SAMPLE = "uniform". Samples SAMPLE_NUM = 10000 Sampling BATCH size = 1000 EPSILON (unit roundoff) = 0.222045E-15 R #1: Row 1 2 Col 1 0.163557 0.128330 2 0.728153 0.318377 3 0.890122 0.622915E-01 4 0.133334 0.652257 5 0.413010 0.582457 6 0.804597 0.946490 7 0.683923E-01 0.967262 8 0.213065 0.428098 9 0.575803 0.123041 10 0.126983 0.689060 Step SEED L2-Change Energy 0 200000 0.320037E-01 Dimension DIM_NUM = 2 Number of points N = 10 Initial SEED = 200000 Current SEED = 568835027 INIT = "random". Max iterations IT_MAX = 0 IT_FIXED (fixed samples) 1 Iterations IT_NUM = 0 Difference IT_DIFF = 0.00000 CVT ENERGY = 0.320037E-01 SAMPLE = "uniform". Samples SAMPLE_NUM = 10000 Sampling BATCH size = 1000 EPSILON (unit roundoff) = 0.222045E-15 R #2: Row 1 2 Col 1 0.212044 0.927361 2 0.223788 0.106077 3 0.207957 0.438716 4 0.443317E-01 0.215025 5 0.626695 0.877505 6 0.432569 0.317168 7 0.563816 0.136534 8 0.856702 0.704179 9 0.598387E-01 0.831223 10 0.242783 0.662314 Step SEED L2-Change Energy 0 100000 0.285261E-01 Dimension DIM_NUM = 2 Number of points N = 10 Initial SEED = 100000 Current SEED = 1647317536 INIT = "random". Max iterations IT_MAX = 0 IT_FIXED (fixed samples) 1 Iterations IT_NUM = 0 Difference IT_DIFF = 0.00000 CVT ENERGY = 0.285261E-01 SAMPLE = "uniform". Samples SAMPLE_NUM = 10000 Sampling BATCH size = 1000 EPSILON (unit roundoff) = 0.222045E-15 R #3: Row 1 2 Col 1 0.163557 0.128330 2 0.728153 0.318377 3 0.890122 0.622915E-01 4 0.133334 0.652257 5 0.413010 0.582457 6 0.804597 0.946490 7 0.683923E-01 0.967262 8 0.213065 0.428098 9 0.575803 0.123041 10 0.126983 0.689060 TEST13 CVT computes a Centroidal Voronoi Tessellation. In this example, we call the "USER" routine, which allows the user to define the geometry and density implicitly, by returning sample points. Step SEED L2-Change Energy 0 123456789 0.109796E-01 1 123456789 4.52097 0.108136E-01 2 123456789 1.81246 0.684193E-02 3 123456789 1.28087 0.623030E-02 4 123456789 1.03029 0.604953E-02 5 123456789 0.950720 0.590688E-02 6 123456789 0.917005 0.587044E-02 7 123456789 0.878544 0.584668E-02 8 123456789 0.817168 0.573118E-02 9 123456789 0.897585 0.572022E-02 10 123456789 0.846693 0.566964E-02 11 123456789 0.835267 0.557785E-02 12 123456789 0.794614 0.554170E-02 13 123456789 0.742569 0.552092E-02 14 123456789 0.769460 0.558744E-02 15 123456789 0.767798 0.545017E-02 16 123456789 0.808586 0.549936E-02 17 123456789 0.820742 0.543456E-02 18 123456789 0.805101 0.544275E-02 19 123456789 0.776684 0.542899E-02 20 123456789 0.779336 0.549018E-02 21 123456789 0.698205 0.548805E-02 22 123456789 0.815561 0.542540E-02 23 123456789 0.697345 0.541393E-02 24 123456789 0.718998 0.537188E-02 25 123456789 0.736981 0.539708E-02 26 123456789 0.801012 0.540335E-02 27 123456789 0.773023 0.537162E-02 28 123456789 0.698042 0.534773E-02 29 123456789 0.700234 0.534826E-02 30 123456789 0.742788 0.538991E-02 31 123456789 0.814323 0.537255E-02 32 123456789 0.758652 0.537922E-02 33 123456789 0.807691 0.532668E-02 34 123456789 0.728024 0.527605E-02 35 123456789 0.735531 0.532839E-02 36 123456789 0.730882 0.539136E-02 37 123456789 0.793012 0.534060E-02 38 123456789 0.762100 0.531227E-02 39 123456789 0.733865 0.534647E-02 40 123456789 0.690335 0.529368E-02 Dimension DIM_NUM = 2 Number of points N = 100 Initial SEED = 123456789 Current SEED = 123456789 INIT = "user". Max iterations IT_MAX = 40 IT_FIXED (fixed samples) 1 Iterations IT_NUM = 40 Difference IT_DIFF = 0.690335 CVT ENERGY = 0.529368E-02 SAMPLE = "user". Samples SAMPLE_NUM = 10000 Sampling BATCH size = 1000 EPSILON (unit roundoff) = 0.222045E-15 TEST14 Generate a CVT in the interval [0,1] using 10 points. Step SEED L2-Change Energy 0 123456789 0.157708E-02 1 529581864 0.163938 0.158402E-02 2 178038614 0.411995E-01 0.113365E-02 3 1965902505 0.274862E-01 0.110224E-02 4 1443447695 0.227305E-01 0.107841E-02 5 631623582 0.224778E-01 0.106463E-02 6 2068398041 0.211534E-01 0.105417E-02 7 606376808 0.209514E-01 0.104211E-02 8 1535406263 0.197632E-01 0.103653E-02 9 1705765739 0.195414E-01 0.102219E-02 10 684646033 0.175148E-01 0.100682E-02 11 75217582 0.183635E-01 0.100507E-02 12 2112829574 0.164685E-01 0.996846E-03 13 284023060 0.166623E-01 0.987710E-03 14 1401063163 0.167631E-01 0.978045E-03 15 1808424944 0.150309E-01 0.974279E-03 16 374964289 0.168720E-01 0.960649E-03 17 691171749 0.170846E-01 0.958608E-03 18 2012255916 0.158204E-01 0.950441E-03 19 588882215 0.147841E-01 0.944058E-03 20 1627928220 0.152235E-01 0.941969E-03 21 1109588173 0.138421E-01 0.931500E-03 22 1469366331 0.148544E-01 0.934350E-03 23 669813744 0.149486E-01 0.925906E-03 24 1768424989 0.122674E-01 0.918672E-03 25 1887537795 0.142293E-01 0.920140E-03 26 92804193 0.141567E-01 0.918602E-03 27 1379446580 0.113628E-01 0.907351E-03 28 1579108318 0.125077E-01 0.907878E-03 29 862792735 0.122400E-01 0.901990E-03 30 218821533 0.107816E-01 0.903236E-03 31 1848001137 0.133499E-01 0.899408E-03 32 413943529 0.133959E-01 0.894013E-03 33 352039253 0.117419E-01 0.892683E-03 34 1837212886 0.103276E-01 0.886067E-03 35 379495084 0.103393E-01 0.884781E-03 36 1226962997 0.973531E-02 0.883175E-03 37 958619815 0.111526E-01 0.878220E-03 38 2100271157 0.111416E-01 0.876237E-03 39 1835842199 0.105023E-01 0.878514E-03 40 1366517430 0.109901E-01 0.873609E-03 Dimension DIM_NUM = 1 Number of points N = 10 Initial SEED = 123456789 Current SEED = 1334318103 INIT = "uniform". Max iterations IT_MAX = 40 IT_FIXED (fixed samples) 1 Iterations IT_NUM = 40 Difference IT_DIFF = 0.109901E-01 CVT ENERGY = 0.873609E-03 SAMPLE = "uniform". Samples SAMPLE_NUM = 100000 Sampling BATCH size = 10000 EPSILON (unit roundoff) = 0.222045E-15 Generators (rows): Row 1 Col 1 0.300814 2 0.941731 3 0.826243 4 0.601668 5 0.495943 6 0.125042 7 0.395508 8 0.211475 9 0.409015E-01 10 0.711881 CVT_PRB Normal end of execution. September 8 2010 5:29:18.361 PM