14 April 2014 09:21:52 AM
PYRAMID_INTEGRALS_PRB
C++ version
Test the PYRAMID_INTEGRALS library.
TEST01
Compare exact and estimated integrals
over the unit pyramid in 3D.
Number of sample points used is 500000
E1 E2 E3 MC-Estimate Exact Error
0 0 0 1.33333 1.33333 0
2 0 0 0.265432 0.266667 0.00123512
4 0 0 0.113616 0.114286 0.000669883
6 0 0 0.0630826 0.0634921 0.000409441
0 2 0 0.266404 0.266667 0.000263013
2 2 0 0.0631172 0.0634921 0.000374887
4 2 0 0.0294246 0.0296296 0.000204996
0 4 0 0.114151 0.114286 0.000134259
2 4 0 0.0294431 0.0296296 0.00018654
0 6 0 0.0633977 0.0634921 9.43871e-05
0 0 1 0.333908 0.333333 0.000575074
2 0 1 0.0443021 0.0444444 0.000142295
4 0 1 0.0142281 0.0142857 5.75918e-05
0 2 1 0.0445419 0.0444444 9.74194e-05
2 2 1 0.00791511 0.00793651 2.14007e-05
0 4 1 0.0143255 0.0142857 3.97702e-05
0 0 2 0.133741 0.133333 0.000407178
2 0 2 0.0126683 0.0126984 3.01348e-05
4 0 2 0.00316495 0.0031746 9.65166e-06
0 2 2 0.0127463 0.0126984 4.78405e-05
2 2 2 0.00176164 0.00176367 2.03053e-06
0 4 2 0.00318874 0.0031746 1.41372e-05
0 0 3 0.066929 0.0666667 0.000262307
2 0 3 0.00475339 0.0047619 8.51487e-06
0 2 3 0.0047852 0.0047619 2.32905e-05
0 0 4 0.0382648 0.0380952 0.000169557
2 0 4 0.00211345 0.0021164 2.94978e-06
0 2 4 0.00212875 0.0021164 1.23501e-05
0 0 5 0.0239229 0.0238095 0.000113393
0 0 6 0.0159526 0.015873 7.96185e-05
TEST02
Print sample points in the unit pyramid in 3D.
Unit pyramid points
Row: 0 1 2
Col
0: -0.529145 0.857507 0.0604056
1: 0.0498951 -0.0684939 0.595634
2: -0.170945 -0.27504 0.647423
3: 0.205248 -0.671476 0.233952
4: -0.183031 0.472301 0.0727328
5: -0.702641 0.560667 0.294766
6: -0.77088 -0.924746 0.0493645
7: 0.133511 -0.147631 0.804149
8: -0.23222 0.399618 0.515596
9: -0.191112 0.415901 0.355966
10: 0.36547 0.117333 0.0485835
11: -0.0778534 0.694145 0.157545
12: -0.457207 0.384369 0.265047
13: -0.450616 0.106489 0.284023
14: 0.139976 -0.165375 0.402856
15: 0.286934 -0.511457 0.331066
16: 0.257074 0.256803 0.604505
17: 0.387721 -0.757334 0.0802734
18: 0.473565 -0.369359 0.0347348
19: 0.390594 0.387917 0.258939
PYRAMID_INTEGRALS_PRB
Normal end of execution.
14 April 2014 09:21:52 AM