Tue May 20 22:50:55 2025 wedge_integrals_test(): python version: 3.10.12 numpy version: 1.26.4 Test wedge_integrals(). wedge01_monomial_integral_test: wedge01_monomial_integral computes the integral of a monomial over the interior of the unit wedge in 3D. Compare with a Monte Carlo estimate. Number of sample points used is 500000 E1 E2 E3 MC-Estimate Exact Error 0 1 0 0.333193 0.333333 0.000139958 1 1 0 0.0833124 0.0833333 2.09464e-05 2 1 0 0.0332991 0.0333333 3.42668e-05 3 1 0 0.0166282 0.0166667 3.85068e-05 4 1 0 0.00948784 0.00952381 3.59703e-05 5 1 0 0.00592125 0.00595238 3.11334e-05 0 2 0 0.166403 0.166667 0.000263401 1 2 0 0.0333112 0.0333333 2.21144e-05 2 2 0 0.0111036 0.0111111 7.48022e-06 3 2 0 0.0047554 0.0047619 6.50839e-06 4 2 0 0.00237513 0.00238095 5.82403e-06 0 3 0 0.0997083 0.1 0.000291697 1 3 0 0.016643 0.0166667 2.36781e-05 2 3 0 0.0047586 0.0047619 3.30574e-06 3 3 0 0.00178435 0.00178571 1.36232e-06 0 4 0 0.0663771 0.0666667 0.000289541 1 4 0 0.0095007 0.00952381 2.31102e-05 2 4 0 0.00237848 0.00238095 2.46787e-06 0 5 0 0.0473442 0.047619 0.000274813 1 5 0 0.00593102 0.00595238 2.13646e-05 0 6 0 0.0354597 0.0357143 0.000254581 0 1 1 -0.000504575 0 0.000504575 1 1 1 -0.00011204 0 0.00011204 2 1 1 -4.75269e-05 0 4.75269e-05 3 1 1 -2.67764e-05 0 2.67764e-05 4 1 1 -1.73689e-05 0 1.73689e-05 0 2 1 -0.000329273 0 0.000329273 1 2 1 -4.57298e-05 0 4.57298e-05 2 2 1 -1.11001e-05 0 1.11001e-05 3 2 1 -4.34198e-06 0 4.34198e-06 0 3 1 -0.000242951 0 0.000242951 1 3 1 -2.88945e-05 0 2.88945e-05 2 3 1 -4.62492e-06 0 4.62492e-06 0 4 1 -0.000183781 0 0.000183781 1 4 1 -2.1136e-05 0 2.1136e-05 0 5 1 -0.00013972 0 0.00013972 0 1 2 0.111011 0.111111 0.00010037 1 1 2 0.027741 0.0277778 3.67666e-05 2 1 2 0.011089 0.0111111 2.20966e-05 3 1 2 0.00553675 0.00555556 1.88041e-05 0 2 2 0.0554237 0.0555556 0.000131867 1 2 2 0.0110821 0.0111111 2.90003e-05 2 2 2 0.00369492 0.0037037 8.78217e-06 0 3 2 0.0332003 0.0333333 0.00013302 1 3 2 0.00553199 0.00555556 2.35659e-05 0 4 2 0.0220954 0.0222222 0.000126792 0 1 3 -0.000335875 0 0.000335875 1 1 3 -7.95863e-05 0 7.95863e-05 2 1 3 -3.23691e-05 0 3.23691e-05 0 2 3 -0.000218622 0 0.000218622 1 2 3 -3.15709e-05 0 3.15709e-05 0 3 3 -0.000159515 0 0.000159515 0 1 4 0.0665805 0.0666667 8.62027e-05 1 1 4 0.0166366 0.0166667 3.00301e-05 0 2 4 0.0332471 0.0333333 8.62304e-05 0 1 5 -0.000253302 0 0.000253302 wedge01_sample_test() wedge01_sample() samples the unit wedge. sample points in the unit wedge. Row: 0 1 2 Col 0 : 0.532444 0.138659 -0.572463 1 : 0.0235064 0.46044 -0.525377 2 : 0.560104 0.188683 -0.682717 3 : 0.219528 0.455945 -0.0953069 4 : 0.502218 0.140575 -0.176278 5 : 0.617489 0.145192 0.0557084 6 : 0.352097 0.423943 0.019731 7 : 0.596342 0.345853 0.667925 8 : 0.355955 0.356404 0.313608 9 : 0.523031 0.226913 -0.917112 wedge01_volume_test(): wedge01_volume() returns the volume of the unit wedge. wedge01_volume() = 1 wedge_integrals_test(): Normal end of execution. Tue May 20 22:51:27 2025