Tue May 20 22:25:09 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) -15 -1.0 6 1.0 -53 -1.0 -3 -1.0 15 -1.0 -71 -1.0 -65 -1.0 9 -1.0 -68 1.0 37 -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.266265 0.266667 0.0004 4 0 0 0.114052 0.114286 0.00023 6 0 0 0.0633533 0.0634921 0.00014 0 2 0 0.266158 0.266667 0.00051 2 2 0 0.0634443 0.0634921 4.8e-05 4 2 0 0.0296398 0.0296296 1e-05 0 4 0 0.114016 0.114286 0.00027 2 4 0 0.0296458 0.0296296 1.6e-05 0 6 0 0.0633335 0.0634921 0.00016 0 0 1 0.333259 0.333333 7.4e-05 2 0 1 0.0443523 0.0444444 9.2e-05 4 0 1 0.0142575 0.0142857 2.8e-05 0 2 1 0.0443462 0.0444444 9.8e-05 2 2 1 0.00791707 0.00793651 1.9e-05 0 4 1 0.0142401 0.0142857 4.6e-05 0 0 2 0.133355 0.133333 2.2e-05 2 0 2 0.0126625 0.0126984 3.6e-05 4 0 2 0.00316591 0.0031746 8.7e-06 0 2 2 0.0126783 0.0126984 2e-05 2 2 2 0.00175837 0.00176367 5.3e-06 0 4 2 0.00316532 0.0031746 9.3e-06 0 0 3 0.0667296 0.0666667 6.3e-05 2 0 3 0.00474695 0.0047619 1.5e-05 0 2 3 0.00475931 0.0047619 2.6e-06 0 0 4 0.0381668 0.0380952 7.2e-05 2 0 4 0.00210945 0.0021164 7e-06 0 2 4 0.00211711 0.0021164 7.1e-07 0 0 5 0.0238777 0.0238095 6.8e-05 0 0 6 0.0159341 0.015873 6.1e-05 pyramid01_sample_test() pyramid01_sample samples points from the unit pyramid. Unit pyramid points Row: 0 1 2 Col 0 : -0.177578 0.046901 0.237956 1 : -0.160188 0.0469698 0.824409 2 : 0.142021 0.574213 0.211501 3 : -0.348298 -0.0378194 0.61313 4 : -0.127889 0.075787 0.268402 5 : 0.430292 -0.101843 0.442664 6 : 0.0689281 0.146816 0.851143 7 : -0.359803 0.567387 0.325468 8 : 0.228758 -0.124012 0.163705 9 : 0.309093 -0.36456 0.157244 10 : 0.7704 -0.258528 0.177129 11 : -0.121594 -0.535706 0.44439 12 : -0.176869 -0.313598 0.0400958 13 : -0.226502 0.458833 0.462016 14 : -0.139229 0.244872 0.0617008 15 : -0.30229 -0.389349 0.210707 16 : -0.0223588 -0.0626591 0.575056 17 : 0.786511 0.881328 0.0977132 18 : 0.460309 0.196665 0.287708 19 : 0.464635 0.416494 0.418359 pyramid01_volume_test() pyramid01_volume() returns the volume of the unit pyramid. pyramid01_volume() = 1.33333 pyramid_integrals_test(): Normal end of execution. Tue May 20 22:25:21 2025