18 August 2014 12:34:33 PM
WEDGE_MONTE_CARLO_PRB
C++ version
Test the WEDGE_MONTE_CARLO library.
TEST01
Use WEDGE01_SAMPLE for a Monte Carlo estimate of an
integral over the interior of the unit wedge in 3D.
N 1 X Y Z X^2 XY Z^2 X^3
1 1 0.867886 0.0254803 0.123391 0.753226 0.022114 0.0152253 0.653715
2 1 0.334286 0.310076 -0.440504 0.136338 0.0662521 0.30936 0.062017
4 1 0.277164 0.14829 -0.812702 0.139137 0.0375949 0.749292 0.0793129
8 1 0.315984 0.337347 -0.141956 0.161208 0.0681841 0.130774 0.103442
16 1 0.307463 0.315327 0.19375 0.140347 0.0812156 0.447464 0.0778203
32 1 0.308952 0.355892 -0.145337 0.135733 0.0916267 0.315629 0.0723115
64 1 0.32878 0.304804 0.130211 0.166761 0.0759253 0.314265 0.105102
128 1 0.338684 0.318797 -0.0144669 0.174893 0.0732268 0.361755 0.108481
256 1 0.326424 0.350059 -0.0452982 0.165104 0.0847417 0.312929 0.10076
512 1 0.310707 0.343038 0.0156344 0.151195 0.0779385 0.314584 0.0892454
1024 1 0.336508 0.328384 -0.020843 0.171212 0.0817678 0.314866 0.104494
2048 1 0.337542 0.334927 -0.0171698 0.168322 0.0868239 0.336835 0.100226
4096 1 0.332109 0.335048 -0.00508344 0.165275 0.0832678 0.326477 0.0988976
8192 1 0.332154 0.329244 0.00490948 0.165365 0.0824993 0.338957 0.0989495
16384 1 0.331892 0.336805 -0.00264762 0.166759 0.0828449 0.333197 0.100773
32768 1 0.332006 0.332594 0.000155373 0.165295 0.0828243 0.330301 0.0988825
65536 1 0.332382 0.334483 -0.00713768 0.165773 0.0835444 0.333493 0.0991654
Exact 1 0.333333 0.333333 0 0.166667 0.0833333 0.333333 0.1
WEDGE_MONTE_CARLO_PRB
Normal end of execution.
18 August 2014 12:34:33 PM