Wed Oct 8 09:16:14 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.333135 0.333333 0.000197884 1 1 0 0.0832902 0.0833333 4.31051e-05 2 1 0 0.0333224 0.0333333 1.0911e-05 3 1 0 0.0166691 0.0166667 2.39139e-06 4 1 0 0.00953118 0.00952381 7.37427e-06 5 1 0 0.00596123 0.00595238 8.85313e-06 0 2 0 0.166419 0.166667 0.000247463 1 2 0 0.0332945 0.0333333 3.88687e-05 2 2 0 0.0110971 0.0111111 1.40525e-05 3 2 0 0.00475614 0.0047619 5.762e-06 4 2 0 0.00237866 0.00238095 2.29733e-06 0 3 0 0.0997729 0.1 0.000227094 1 3 0 0.0166394 0.0166667 2.72209e-05 2 3 0 0.00475363 0.0047619 8.2764e-06 3 3 0 0.00178231 0.00178571 3.40559e-06 0 4 0 0.0664744 0.0666667 0.000192262 1 4 0 0.00950521 0.00952381 1.86037e-05 2 4 0 0.00237639 0.00238095 4.56204e-06 0 5 0 0.0474593 0.047619 0.000159716 1 5 0 0.00593946 0.00595238 1.29193e-05 0 6 0 0.0355811 0.0357143 0.000133178 0 1 1 -0.000221732 0 0.000221732 1 1 1 -7.74111e-05 0 7.74111e-05 2 1 1 -2.99751e-05 0 2.99751e-05 3 1 1 -1.46266e-05 0 1.46266e-05 4 1 1 -8.17663e-06 0 8.17663e-06 0 2 1 -0.000121788 0 0.000121788 1 2 1 -4.97439e-05 0 4.97439e-05 2 2 1 -1.66382e-05 0 1.66382e-05 3 2 1 -6.66935e-06 0 6.66935e-06 0 3 1 -4.16893e-05 0 4.16893e-05 1 3 1 -3.00025e-05 0 3.00025e-05 2 3 1 -9.52162e-06 0 9.52162e-06 0 4 1 1.00878e-05 0 1.00878e-05 1 4 1 -1.75661e-05 0 1.75661e-05 0 5 1 4.10982e-05 0 4.10982e-05 0 1 2 0.110929 0.111111 0.000181743 1 1 2 0.0277201 0.0277778 5.76895e-05 2 1 2 0.0110795 0.0111111 3.16176e-05 3 1 2 0.005538 0.00555556 1.75531e-05 0 2 2 0.055491 0.0555556 6.45959e-05 1 2 2 0.0111038 0.0111111 7.32662e-06 2 2 2 0.00369843 0.0037037 5.2691e-06 0 3 2 0.033301 0.0333333 3.22966e-05 1 3 2 0.00555684 0.00555556 1.2798e-06 0 4 2 0.0222046 0.0222222 1.76111e-05 0 1 3 -2.33822e-05 0 2.33822e-05 1 1 3 -2.67156e-05 0 2.67156e-05 2 1 3 -1.68988e-05 0 1.68988e-05 0 2 3 -4.05889e-06 0 4.05889e-06 1 2 3 -1.72582e-05 0 1.72582e-05 0 3 3 1.88359e-05 0 1.88359e-05 0 1 4 0.0665338 0.0666667 0.000132894 1 1 4 0.016607 0.0166667 5.96365e-05 0 2 4 0.0332851 0.0333333 4.82568e-05 0 1 5 4.6107e-05 0 4.6107e-05 wedge01_sample_test() wedge01_sample() samples the unit wedge. sample points in the unit wedge. Row: 0 1 2 Col 0 : 0.303575 0.395838 -0.878963 1 : 0.745241 0.169339 -0.989862 2 : 0.339618 0.0946823 -0.706132 3 : 0.181354 0.666026 -0.54914 4 : 0.236024 0.527413 0.172447 5 : 0.377796 0.034943 0.894296 6 : 0.0495765 0.705261 0.704512 7 : 0.378577 0.521163 0.931453 8 : 0.323712 0.586094 0.0995401 9 : 0.175444 0.321332 0.160284 wedge01_volume_test(): wedge01_volume() returns the volume of the unit wedge. wedge01_volume() = 1 wedge_integrals_test(): Normal end of execution. Wed Oct 8 09:16:46 2025