Fri Jul 14 15:49:02 2023 prism_witherden_rule_test(): Python version: 3.8.10 Test prism_witherden_rule(). prism_witherden_rule_test01(): Quadrature rule for the unit prism, Precision p = 5 I W X Y Z 0 0.177864 0.333333 0.333333 0.5 1 0.0561775 0.4873 0.0254001 0.5 2 0.0561775 0.0254001 0.4873 0.5 3 0.0561775 0.4873 0.4873 0.5 4 0.0464154 0.445598 0.108804 0.0644985 5 0.0464154 0.108804 0.445598 0.0644985 6 0.0464154 0.445598 0.445598 0.0644985 7 0.0464154 0.445598 0.108804 0.935501 8 0.0464154 0.108804 0.445598 0.935501 9 0.0464154 0.445598 0.445598 0.935501 10 0.0625186 0.100859 0.798282 0.214787 11 0.0625186 0.798282 0.100859 0.214787 12 0.0625186 0.100859 0.100859 0.214787 13 0.0625186 0.100859 0.798282 0.785213 14 0.0625186 0.798282 0.100859 0.785213 15 0.0625186 0.100859 0.100859 0.785213 Weight Sum 0.9999999999999998 prism_witherden_rule_test02(): Test the precision of a quadrature rule for the unit prism. Stated precision of rule = 5 Number of quadrature points = 16 Degree Maximum error 0 1.110223024625157e-16 1 2.775557561562891e-17 2 2.775557561562891e-17 3 5.551115123125783e-17 4 5.551115123125783e-17 5 4.163336342344337e-17 6 8.681433560950649e-05 7 0.0002425567824553866 prism_witherden_rule_test03(): Test the precision of quadrature rules for the unit prism. Check rules of precision p = 0 through 10 for error in approximating integrals of monomials. maximum maximum p absolute relative error error 0 0 0 1 5.551115123125783e-17 3.33066907387547e-16 2 2.775557561562891e-17 1.665334536937735e-16 3 1.110223024625157e-16 2.220446049250313e-16 4 5.551115123125783e-17 3.33066907387547e-16 5 1.110223024625157e-16 6.245004513516506e-16 6 6.661338147750939e-16 2.137179322403426e-15 7 3.33066907387547e-16 1.311450947838466e-15 8 1.387778780781446e-16 9.107298248878237e-16 9 1.110223024625157e-16 1.748601263784622e-15 10 2.220446049250313e-16 2.629732369363591e-15 prism_witherden_rule_test(): Normal end of execution. Fri Jul 14 15:49:02 2023