>> cd tetrahedron_monte_carlo >> tetrahedron_monte_carlo_test 16-Aug-2009 20:10:13 TETRAHEDRON_MONTE_CARLO_TEST MATLAB version TETRAHEDRON_MONTE_CARLO_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 ??? Input argument "p_num" is undefined. Error in ==> tetrahedron_unit_sample_01 at 39 for j = 1, p_num Error in ==> tetrahedron_monte_carlo_test01 at 47 [ result, seed ] = tetrahedron_monte_carlo ( t, p_num, f_num, ... Error in ==> tetrahedron_monte_carlo_test at 27 tetrahedron_monte_carlo_test01 ( ); >> tetrahedron_monte_carlo_test 16-Aug-2009 20:11:44 TETRAHEDRON_MONTE_CARLO_TEST MATLAB version TETRAHEDRON_MONTE_CARLO_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.001208 0.005287 0.004586 0.023150 0.020081 0.017418 2 0.020056 0.004813 0.012528 0.024194 0.004276 0.007896 4 0.030227 0.023101 0.007184 0.030715 0.010207 0.007521 8 0.011770 0.010663 0.009501 0.019944 0.011680 0.012357 16 0.013290 0.009028 0.008139 0.010112 0.007301 0.010756 32 0.014743 0.010401 0.010901 0.013642 0.008568 0.015240 64 0.013406 0.009507 0.008114 0.012957 0.008920 0.011106 128 0.014217 0.009762 0.009484 0.015123 0.009460 0.011813 256 0.014430 0.009086 0.009854 0.012381 0.009059 0.015148 512 0.014252 0.009502 0.009347 0.013598 0.008618 0.012520 1024 0.013984 0.009437 0.009159 0.014204 0.009380 0.013934 2048 0.013459 0.009249 0.009376 0.013531 0.009447 0.014294 4096 0.013454 0.009318 0.009380 0.013652 0.009280 0.013900 8192 0.013596 0.009493 0.009269 0.013775 0.009262 0.013458 16384 0.013808 0.009255 0.009373 0.013622 0.009330 0.013637 32768 0.013779 0.009356 0.009329 0.013713 0.009271 0.013459 65536 0.013755 0.009372 0.009404 0.013658 0.009351 0.013759 TETRAHEDRON_MONTE_CARLO_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.001208 0.005287 0.004586 0.023150 0.020081 0.017418 2 0.020056 0.004813 0.012528 0.024194 0.004276 0.007896 4 0.030227 0.023101 0.007184 0.030715 0.010207 0.007521 8 0.011770 0.010663 0.009501 0.019944 0.011680 0.012357 16 0.013290 0.009028 0.008139 0.010112 0.007301 0.010756 32 0.014743 0.010401 0.010901 0.013642 0.008568 0.015240 64 0.013406 0.009507 0.008114 0.012957 0.008920 0.011106 128 0.014217 0.009762 0.009484 0.015123 0.009460 0.011813 256 0.014430 0.009086 0.009854 0.012381 0.009059 0.015148 512 0.014252 0.009502 0.009347 0.013598 0.008618 0.012520 1024 0.013984 0.009437 0.009159 0.014204 0.009380 0.013934 2048 0.013459 0.009249 0.009376 0.013531 0.009447 0.014294 4096 0.013454 0.009318 0.009380 0.013652 0.009280 0.013900 8192 0.013596 0.009493 0.009269 0.013775 0.009262 0.013458 16384 0.013808 0.009255 0.009373 0.013622 0.009330 0.013637 32768 0.013779 0.009356 0.009329 0.013713 0.009271 0.013459 65536 0.013755 0.009372 0.009404 0.013658 0.009351 0.013759 TETRAHEDRON_MONTE_CARLO_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.001208 0.005287 0.004586 0.023150 0.020081 0.017418 2 0.020056 0.004813 0.012528 0.024194 0.004276 0.007896 4 0.030227 0.023101 0.007184 0.030715 0.010207 0.007521 8 0.011770 0.010663 0.009501 0.019944 0.011680 0.012357 16 0.013290 0.009028 0.008139 0.010112 0.007301 0.010756 32 0.014743 0.010401 0.010901 0.013642 0.008568 0.015240 64 0.013406 0.009507 0.008114 0.012957 0.008920 0.011106 128 0.014217 0.009762 0.009484 0.015123 0.009460 0.011813 256 0.014430 0.009086 0.009854 0.012381 0.009059 0.015148 512 0.014252 0.009502 0.009347 0.013598 0.008618 0.012520 1024 0.013984 0.009437 0.009159 0.014204 0.009380 0.013934 2048 0.013459 0.009249 0.009376 0.013531 0.009447 0.014294 4096 0.013454 0.009318 0.009380 0.013652 0.009280 0.013900 8192 0.013596 0.009493 0.009269 0.013775 0.009262 0.013458 16384 0.013808 0.009255 0.009373 0.013622 0.009330 0.013637 32768 0.013779 0.009356 0.009329 0.013713 0.009271 0.013459 65536 0.013755 0.009372 0.009404 0.013658 0.009351 0.013759 TETRAHEDRON_MONTE_CARLO_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.001208 0.005287 0.004586 0.023150 0.020081 0.017418 2 0.020056 0.004813 0.012528 0.024194 0.004276 0.007896 4 0.030227 0.023101 0.007184 0.030715 0.010207 0.007521 8 0.011770 0.010663 0.009501 0.019944 0.011680 0.012357 16 0.013290 0.009028 0.008139 0.010112 0.007301 0.010756 32 0.014743 0.010401 0.010901 0.013642 0.008568 0.015240 64 0.013406 0.009507 0.008114 0.012957 0.008920 0.011106 128 0.014217 0.009762 0.009484 0.015123 0.009460 0.011813 256 0.014430 0.009086 0.009854 0.012381 0.009059 0.015148 512 0.014252 0.009502 0.009347 0.013598 0.008618 0.012520 1024 0.013984 0.009437 0.009159 0.014204 0.009380 0.013934 2048 0.013459 0.009249 0.009376 0.013531 0.009447 0.014294 4096 0.013454 0.009318 0.009380 0.013652 0.009280 0.013900 8192 0.013596 0.009493 0.009269 0.013775 0.009262 0.013458 16384 0.013808 0.009255 0.009373 0.013622 0.009330 0.013637 32768 0.013779 0.009356 0.009329 0.013713 0.009271 0.013459 65536 0.013755 0.009372 0.009404 0.013658 0.009351 0.013759 TETRAHEDRON_MONTE_CARLO_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.000000 2.000000 3.000000 2 4.000000 1.000000 2.000000 3 2.000000 4.000000 4.000000 4 3.000000 2.000000 5.000000 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 TETRAHEDRON_MONTE_CARLO_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.000000 2.000000 3.000000 2 4.000000 1.000000 2.000000 3 2.000000 4.000000 4.000000 4 3.000000 2.000000 5.000000 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 TETRAHEDRON_MONTE_CARLO_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.000000 2.000000 3.000000 2 4.000000 1.000000 2.000000 3 2.000000 4.000000 4.000000 4 3.000000 2.000000 5.000000 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 TETRAHEDRON_MONTE_CARLO_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.000000 2.000000 3.000000 2 4.000000 1.000000 2.000000 3 2.000000 4.000000 4.000000 4 3.000000 2.000000 5.000000 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 TETRAHEDRON_MONTE_CARLO_TEST Normal end of execution. 16-Aug-2009 20:27:10 >>