17 September 2021 11:54:21.266 PM TETRAHEDRON_MONTE_CARLO_NONUNIT_TEST FORTRAN90 version Test the TETRAHEDRON_MONTE_CARLO_NONUNIT library. TEST01 Sample using TETRAHEDRON_UNIT_SAMPLE_01 Integrate TETRAHEDRON_UNIT_INTEGRAND_03 Integration region is the unit tetrahedron. Use an increasing number of points P_NUM. Note that the sample routine is a bad sampler. P_NUM X^2 X*Y X*Z Y^2 Y*Z Z^2 1 0.120763E-02 0.528746E-02 0.458634E-02 0.231505E-01 0.200807E-01 0.174180E-01 2 0.200559E-01 0.481267E-02 0.125281E-01 0.241942E-01 0.427564E-02 0.789573E-02 4 0.302268E-01 0.231009E-01 0.718381E-02 0.307148E-01 0.102074E-01 0.752060E-02 8 0.117702E-01 0.106626E-01 0.950056E-02 0.199444E-01 0.116801E-01 0.123572E-01 16 0.132896E-01 0.902844E-02 0.813869E-02 0.101117E-01 0.730063E-02 0.107560E-01 32 0.147432E-01 0.104011E-01 0.109008E-01 0.136417E-01 0.856791E-02 0.152403E-01 64 0.134063E-01 0.950749E-02 0.811380E-02 0.129571E-01 0.892010E-02 0.111060E-01 128 0.142174E-01 0.976202E-02 0.948439E-02 0.151230E-01 0.946042E-02 0.118131E-01 256 0.144301E-01 0.908613E-02 0.985396E-02 0.123813E-01 0.905926E-02 0.151476E-01 512 0.142517E-01 0.950238E-02 0.934722E-02 0.135979E-01 0.861832E-02 0.125197E-01 1024 0.139837E-01 0.943655E-02 0.915851E-02 0.142043E-01 0.938047E-02 0.139340E-01 2048 0.134590E-01 0.924911E-02 0.937593E-02 0.135310E-01 0.944727E-02 0.142939E-01 4096 0.134540E-01 0.931752E-02 0.937959E-02 0.136518E-01 0.928029E-02 0.139004E-01 8192 0.135961E-01 0.949325E-02 0.926899E-02 0.137752E-01 0.926235E-02 0.134578E-01 16384 0.138082E-01 0.925538E-02 0.937326E-02 0.136220E-01 0.932962E-02 0.136365E-01 32768 0.137786E-01 0.935584E-02 0.932879E-02 0.137130E-01 0.927055E-02 0.134594E-01 65536 0.137548E-01 0.937239E-02 0.940389E-02 0.136580E-01 0.935064E-02 0.137592E-01 TEST02 Sample using TETRAHEDRON_UNIT_SAMPLE_02 Integrate TETRAHEDRON_UNIT_INTEGRAND_03 Integration region is the unit tetrahedron. Use an increasing number of points P_NUM. Note that the sample routine is a good sampler. P_NUM X^2 X*Y X*Z Y^2 Y*Z Z^2 1 0.268006E-02 0.923216E-03 0.138385E-01 0.318026E-03 0.476701E-02 0.714547E-01 2 0.279388E-01 0.203075E-01 0.280425E-02 0.153809E-01 0.189399E-02 0.315035E-03 4 0.252599E-01 0.641717E-02 0.102405E-01 0.194421E-01 0.490141E-02 0.330667E-01 8 0.411817E-01 0.133023E-01 0.815115E-02 0.188628E-01 0.530985E-02 0.406394E-02 16 0.162033E-01 0.546692E-02 0.758590E-02 0.858323E-02 0.786344E-02 0.169424E-01 32 0.127310E-01 0.855491E-02 0.570481E-02 0.236266E-01 0.118318E-01 0.192896E-01 64 0.135433E-01 0.755872E-02 0.745470E-02 0.226530E-01 0.879154E-02 0.175110E-01 128 0.201186E-01 0.936979E-02 0.925194E-02 0.142647E-01 0.854525E-02 0.160820E-01 256 0.172159E-01 0.846661E-02 0.788278E-02 0.167369E-01 0.782860E-02 0.164737E-01 512 0.169019E-01 0.837163E-02 0.844684E-02 0.170946E-01 0.803322E-02 0.158228E-01 1024 0.176071E-01 0.835279E-02 0.851073E-02 0.166157E-01 0.810487E-02 0.162110E-01 2048 0.164696E-01 0.854378E-02 0.815693E-02 0.169153E-01 0.826185E-02 0.163875E-01 4096 0.164776E-01 0.856374E-02 0.823639E-02 0.171876E-01 0.854875E-02 0.168340E-01 8192 0.170230E-01 0.836067E-02 0.839444E-02 0.162770E-01 0.835697E-02 0.165180E-01 16384 0.166148E-01 0.840943E-02 0.831745E-02 0.165876E-01 0.837075E-02 0.168483E-01 32768 0.165732E-01 0.833576E-02 0.829823E-02 0.166856E-01 0.834954E-02 0.165646E-01 65536 0.166837E-01 0.832516E-02 0.834317E-02 0.167211E-01 0.828953E-02 0.166479E-01 TEST03 Sample using TETRAHEDRON_UNIT_SAMPLE_03 Integrate TETRAHEDRON_UNIT_INTEGRAND_03 Integration region is the unit tetrahedron. Use an increasing number of points P_NUM. Note that the sample routine is a good sampler. P_NUM X^2 X*Y X*Z Y^2 Y*Z Z^2 1 0.263700E-01 0.881743E-03 0.665648E-02 0.294831E-04 0.222575E-03 0.168027E-02 2 0.135749E-01 0.171669E-01 0.133528E-01 0.222438E-01 0.192891E-01 0.239396E-01 4 0.481305E-01 0.142144E-01 0.503491E-02 0.245073E-01 0.383041E-02 0.151703E-02 8 0.140793E-01 0.496968E-02 0.153538E-01 0.118848E-01 0.989156E-02 0.217565E-01 16 0.158204E-01 0.849140E-02 0.705630E-02 0.205308E-01 0.911986E-02 0.238999E-01 32 0.158471E-01 0.719333E-02 0.680947E-02 0.161667E-01 0.897796E-02 0.135959E-01 64 0.224717E-01 0.846362E-02 0.810511E-02 0.153362E-01 0.866248E-02 0.171957E-01 128 0.123761E-01 0.783250E-02 0.820489E-02 0.154607E-01 0.856506E-02 0.181744E-01 256 0.160873E-01 0.866279E-02 0.786049E-02 0.187582E-01 0.874858E-02 0.193828E-01 512 0.156889E-01 0.835191E-02 0.856512E-02 0.171408E-01 0.857428E-02 0.170312E-01 1024 0.173688E-01 0.809484E-02 0.820148E-02 0.152979E-01 0.800980E-02 0.160049E-01 2048 0.170081E-01 0.870171E-02 0.849897E-02 0.165675E-01 0.862529E-02 0.167969E-01 4096 0.165166E-01 0.821636E-02 0.825902E-02 0.164795E-01 0.859219E-02 0.162134E-01 8192 0.164993E-01 0.837596E-02 0.830946E-02 0.169553E-01 0.834807E-02 0.165908E-01 16384 0.167675E-01 0.837675E-02 0.842485E-02 0.164294E-01 0.828018E-02 0.167963E-01 32768 0.170165E-01 0.840034E-02 0.839122E-02 0.165714E-01 0.827328E-02 0.164725E-01 65536 0.166177E-01 0.835167E-02 0.836320E-02 0.166437E-01 0.834775E-02 0.167750E-01 TEST04 Sample using TETRAHEDRON_UNIT_SAMPLE_04 Integrate TETRAHEDRON_UNIT_INTEGRAND_03 Integration region is the unit tetrahedron. Use an increasing number of points P_NUM. Note that the sample routine is a good sampler. P_NUM X^2 X*Y X*Z Y^2 Y*Z Z^2 1 0.710712E-01 0.208659E-02 0.873222E-02 0.612602E-04 0.256370E-03 0.107289E-02 2 0.171202E-01 0.619286E-02 0.160657E-01 0.123334E-01 0.804945E-02 0.155724E-01 4 0.408874E-02 0.514367E-03 0.618915E-02 0.844536E-03 0.239784E-02 0.163394E-01 8 0.129251E-01 0.505735E-02 0.705661E-02 0.146145E-01 0.582353E-02 0.192570E-01 16 0.191599E-01 0.974345E-02 0.100214E-01 0.169617E-01 0.989801E-02 0.194349E-01 32 0.116544E-01 0.824662E-02 0.782170E-02 0.159309E-01 0.645544E-02 0.166135E-01 64 0.157147E-01 0.840319E-02 0.853668E-02 0.157488E-01 0.100792E-01 0.209082E-01 128 0.163329E-01 0.727833E-02 0.896483E-02 0.151086E-01 0.901500E-02 0.153849E-01 256 0.171697E-01 0.850398E-02 0.748752E-02 0.180581E-01 0.824735E-02 0.154249E-01 512 0.155455E-01 0.787818E-02 0.837353E-02 0.178996E-01 0.867139E-02 0.184081E-01 1024 0.169094E-01 0.809834E-02 0.814815E-02 0.162248E-01 0.817180E-02 0.169490E-01 2048 0.173290E-01 0.872283E-02 0.819597E-02 0.162348E-01 0.812960E-02 0.159704E-01 4096 0.166896E-01 0.833840E-02 0.836203E-02 0.166834E-01 0.832732E-02 0.164112E-01 8192 0.167690E-01 0.828496E-02 0.847933E-02 0.163504E-01 0.829380E-02 0.169832E-01 16384 0.166988E-01 0.825950E-02 0.826190E-02 0.170481E-01 0.846335E-02 0.164934E-01 32768 0.165493E-01 0.830354E-02 0.838069E-02 0.166469E-01 0.839485E-02 0.167151E-01 65536 0.164739E-01 0.829776E-02 0.825136E-02 0.167255E-01 0.830295E-02 0.165317E-01 TEST05 Sample using TETRAHEDRON_UNIT_SAMPLE_01 Integrate TETRAHEDRON_UNIT_INTEGRAND_USER Integration region is over a general tetrahedron. Use an increasing number of points P_NUM. Note that the sample routine is a bad sampler. Tetrahedron vertices: Row 1 2 3 Col 1 1.00000 2.00000 3.00000 2 4.00000 1.00000 2.00000 3 2.00000 4.00000 4.00000 4 3.00000 2.00000 5.00000 P_NUM 1 2.666667 7.677802 13.787802 103.740525 200.440978 1312.253750 2 2.666667 6.938161 11.816584 103.243196 150.068891 743.623269 4 2.666667 6.704359 9.832288 62.465590 104.450599 525.091316 8 2.666667 6.886201 12.933546 99.679981 164.233500 921.710407 16 2.666667 6.665851 13.444781 131.040428 198.773150 1137.920542 32 2.666667 6.560239 14.186029 115.167655 180.885655 1029.824747 64 2.666667 6.715444 13.164204 121.007074 186.443549 1056.415307 128 2.666667 6.708003 13.303313 114.847142 179.694899 1016.562371 256 2.666667 6.568157 14.341910 119.389100 188.171742 1089.711169 512 2.666667 6.648467 13.568521 118.347471 183.471437 1039.428166 1024 2.666667 6.670851 13.827819 116.517535 184.197580 1061.060697 2048 2.666667 6.664625 13.996550 118.496219 187.356269 1084.394072 4096 2.666667 6.668053 13.889247 118.206828 186.159822 1071.057169 8192 2.666667 6.679247 13.751302 117.712599 185.204417 1062.855240 16384 2.666667 6.657845 13.835744 118.013330 185.366053 1063.253727 32768 2.666667 6.666450 13.776989 117.856604 185.092830 1060.927415 65536 2.666667 6.658648 13.854599 117.775059 185.324895 1064.423344 TEST06 Sample using TETRAHEDRON_UNIT_SAMPLE_02 Integrate TETRAHEDRON_UNIT_INTEGRAND_USER Integration region is over a general tetrahedron. Use an increasing number of points P_NUM. Note that the sample routine is a good sampler. Tetrahedron vertices: Row 1 2 3 Col 1 1.00000 2.00000 3.00000 2 4.00000 1.00000 2.00000 3 2.00000 4.00000 4.00000 4 3.00000 2.00000 5.00000 P_NUM 1 2.666667 5.694113 28.442319 165.668296 291.701055 2034.199225 2 2.666667 6.514680 8.937770 118.133757 150.884520 718.697240 4 2.666667 6.199145 17.032660 109.214372 169.404021 979.344983 8 2.666667 6.074354 10.864381 91.163120 125.794511 597.643776 16 2.666667 6.582133 14.817567 138.959210 212.065688 1265.815645 32 2.666667 7.086171 13.712946 109.369437 181.802196 1108.657503 64 2.666667 6.993824 13.869338 117.882584 188.286150 1129.830457 128 2.666667 6.447030 14.202862 117.603156 180.628479 1039.809284 256 2.666667 6.668839 14.075610 122.582250 189.245134 1112.020605 512 2.666667 6.662429 13.948270 121.549236 186.055319 1075.615127 1024 2.666667 6.633790 14.086561 121.408880 186.875408 1089.281918 2048 2.666667 6.697528 14.013978 120.808052 187.435992 1097.251829 4096 2.666667 6.688150 14.089732 119.581451 186.116284 1089.330905 8192 2.666667 6.643817 14.137582 121.371246 187.517483 1095.103064 16384 2.666667 6.663181 14.168441 121.177565 187.748181 1098.663591 32768 2.666667 6.673128 14.108860 121.380329 187.791142 1097.591535 65536 2.666667 6.665841 14.128984 121.332467 187.571805 1095.700072 TEST07 Sample using TETRAHEDRON_UNIT_SAMPLE_03 Integrate TETRAHEDRON_UNIT_INTEGRAND_USER Integration region is over a general tetrahedron. Use an increasing number of points P_NUM. Note that the sample routine is a good sampler. Tetrahedron vertices: Row 1 2 3 Col 1 1.00000 2.00000 3.00000 2 4.00000 1.00000 2.00000 3 2.00000 4.00000 4.00000 4 3.00000 2.00000 5.00000 P_NUM 1 2.666667 5.646280 12.760581 179.010725 204.010937 944.926345 2 2.666667 6.590726 14.935970 75.882759 145.818766 877.503837 4 2.666667 6.113695 9.601741 90.959173 117.015621 540.867180 8 2.666667 6.227600 16.808507 126.994141 194.490847 1103.340771 16 2.666667 6.831215 15.060979 113.757266 189.221208 1181.612546 32 2.666667 6.758450 13.476394 124.801124 190.872900 1109.538339 64 2.666667 6.458926 14.327367 116.526385 181.916916 1076.530277 128 2.666667 6.736866 14.667483 127.749147 199.750429 1188.253632 256 2.666667 6.728306 14.418535 115.851672 184.599302 1098.269254 512 2.666667 6.703776 14.200561 121.123832 188.357604 1104.283444 1024 2.666667 6.622420 14.046198 123.947422 188.925505 1095.625916 2048 2.666667 6.643654 14.139573 119.072649 185.098911 1079.839320 4096 2.666667 6.681554 14.063617 121.846048 188.593370 1103.270378 8192 2.666667 6.680693 14.093908 120.889151 187.248768 1094.263908 16384 2.666667 6.653646 14.175600 121.389161 187.728982 1096.719417 32768 2.666667 6.650172 14.102911 121.016268 186.843567 1089.393131 65536 2.666667 6.662744 14.166102 121.437893 187.896932 1098.789017 TEST08 Sample using TETRAHEDRON_UNIT_SAMPLE_04 Integrate TETRAHEDRON_UNIT_INTEGRAND_USER Integration region is over a general tetrahedron. Use an increasing number of points P_NUM. Note that the sample routine is a good sampler. Tetrahedron vertices: Row 1 2 3 Col 1 1.00000 2.00000 3.00000 2 4.00000 1.00000 2.00000 3 2.00000 4.00000 4.00000 4 3.00000 2.00000 5.00000 P_NUM 1 2.666667 4.354428 12.227043 119.932177 117.919463 412.310379 2 2.666667 6.329140 15.148304 113.207805 174.273903 944.218753 4 2.666667 6.830304 16.787164 222.123933 319.213985 2031.767418 8 2.666667 6.771477 14.758498 132.574443 202.633107 1188.023345 16 2.666667 6.535086 14.810590 115.318790 183.953720 1102.360868 32 2.666667 6.859090 14.023154 129.133011 198.553146 1160.635898 64 2.666667 6.573301 15.111860 119.660683 190.144808 1138.226186 128 2.666667 6.644415 14.173193 126.706319 193.130570 1122.166906 256 2.666667 6.746496 13.683588 120.799047 186.569274 1092.276434 512 2.666667 6.725338 14.446927 120.380519 189.234633 1120.687078 1024 2.666667 6.653303 14.218105 122.353554 189.050661 1107.796121 2048 2.666667 6.651354 13.968910 120.940358 186.565773 1086.549983 4096 2.666667 6.671532 14.084542 121.305812 187.560627 1095.226604 8192 2.666667 6.648014 14.256495 122.065406 188.938007 1107.557730 16384 2.666667 6.682070 14.065573 120.757633 186.970964 1092.555809 32768 2.666667 6.667073 14.163238 121.352006 187.983853 1099.805288 65536 2.666667 6.679438 14.100596 121.788998 188.274269 1100.802684 TETRAHEDRON_MONTE_CARLO_NONUNIT_TEST Normal end of execution. 17 September 2021 11:54:21.382 PM