Wed Oct 8 08:48:05 2025 pyramid_integrals_test(): python version: 3.10.12 numpy version: 1.26.4 Test pyramid_integrals(). r8_mop_test(): r8_mop() evaluates (-1.0)^I4 as an R8. I4 r8_mop(I4) 39 -1.0 95 -1.0 83 -1.0 -60 1.0 14 1.0 -23 -1.0 83 -1.0 -64 1.0 -21 -1.0 -93 -1.0 pyramid01_monomial_integral_test pyramid01_monomial_integral returns the integral of a monomial over the unit pyramid in 3D. Compare to a Monte Carlo estimate. 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.266632 0.266667 3.4e-05 4 0 0 0.114305 0.114286 2e-05 6 0 0 0.063499 0.0634921 6.9e-06 0 2 0 0.266332 0.266667 0.00033 2 2 0 0.063427 0.0634921 6.5e-05 4 2 0 0.0295635 0.0296296 6.6e-05 0 4 0 0.114201 0.114286 8.4e-05 2 4 0 0.029599 0.0296296 3.1e-05 0 6 0 0.063506 0.0634921 1.4e-05 0 0 1 0.333197 0.333333 0.00014 2 0 1 0.0444534 0.0444444 9e-06 4 0 1 0.0142972 0.0142857 1.2e-05 0 2 1 0.0443527 0.0444444 9.2e-05 2 2 1 0.00793295 0.00793651 3.6e-06 0 4 1 0.0142615 0.0142857 2.4e-05 0 0 2 0.133223 0.133333 0.00011 2 0 2 0.0127094 0.0126984 1.1e-05 4 0 2 0.00317728 0.0031746 2.7e-06 0 2 2 0.0126681 0.0126984 3e-05 2 2 2 0.00176323 0.00176367 4.4e-07 0 4 2 0.00316754 0.0031746 7.1e-06 0 0 3 0.0665556 0.0666667 0.00011 2 0 3 0.00476965 0.0047619 7.7e-06 0 2 3 0.00475006 0.0047619 1.2e-05 0 0 4 0.0379843 0.0380952 0.00011 2 0 4 0.00212073 0.0021164 4.3e-06 0 2 4 0.00211056 0.0021164 5.8e-06 0 0 5 0.0237041 0.0238095 0.00011 0 0 6 0.0157759 0.015873 9.7e-05 pyramid01_sample_test() pyramid01_sample samples points from the unit pyramid. Unit pyramid points Row: 0 1 2 Col 0 : -0.818512 0.678079 0.119521 1 : -0.416863 0.119279 0.540248 2 : 0.0619124 -0.484026 0.0593566 3 : -0.277272 0.0554696 0.530571 4 : -0.341784 -0.921529 0.0428079 5 : -0.0740712 -0.203287 0.377582 6 : 0.233132 0.679394 0.00856659 7 : -0.321909 -0.0982764 0.600932 8 : 0.70316 0.189554 0.220206 9 : 0.0256629 0.356223 0.250923 10 : -0.442291 -0.843888 0.0752897 11 : -0.402119 0.427897 0.0430591 12 : -0.495506 0.195636 0.324873 13 : 0.208197 -0.459343 0.0341642 14 : -0.120759 0.372568 0.502976 15 : 0.48271 0.358155 0.209893 16 : -0.345747 -0.178034 0.0741858 17 : -0.327271 0.15836 0.378707 18 : -0.560817 0.620792 0.105151 19 : -0.0796307 -0.696134 0.113065 pyramid01_volume_test() pyramid01_volume() returns the volume of the unit pyramid. pyramid01_volume() = 1.33333 pyramid_integrals_test(): Normal end of execution. Wed Oct 8 08:48:18 2025