Thu Jul 13 13:46:22 2023 pyramid_witherden_rule_test(): Python version: 3.8.10 Test pyramid_witherden_rule(). pyramid_witherden_rule_test01(): Quadrature rule for the unit pyramid, Precision p = 5 I W X Y Z 0 0.0684354 0 0 0.729858 1 0.169397 0 0 0.300401 2 0.0587045 0 0 6.49177e-09 3 0.0764413 0.706526 0 0.125 4 0.0764413 0 0.706526 0.125 5 0.0764413 -0.706526 0 0.125 6 0.0764413 0 -0.706526 0.125 7 0.0396709 0.705117 0.705117 0.0611119 8 0.0396709 0.705117 -0.705117 0.0611119 9 0.0396709 -0.705117 0.705117 0.0611119 10 0.0396709 -0.705117 -0.705117 0.0611119 11 0.0597536 0.432883 0.432883 0.423601 12 0.0597536 0.432883 -0.432883 0.423601 13 0.0597536 -0.432883 0.432883 0.423601 14 0.0597536 -0.432883 -0.432883 0.423601 Weight Sum 1.0000000000000002 pyramid_witherden_rule_test02(): Test the precision of a quadrature rule for the unit pyramid. Stated precision of rule = 5 Number of quadrature points = 15 Degree Maximum error 0 2.220446049250313e-16 1 1.110223024625157e-16 2 5.551115123125783e-17 3 1.52655665885959e-16 4 1.266348137463069e-16 5 3.122502256758253e-17 6 0.0100361111063078 7 0.0007024040490240255 pyramid_witherden_rule_test03(): Test the precision of quadrature rules for the unit pyramid. 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 1.665334536937735e-16 4.996003610813202e-16 2 6.661338147750939e-16 4.996003610813204e-16 3 1.249000902703301e-16 1.873501354054955e-15 4 2.220446049250313e-16 9.835882108788589e-15 5 2.220446049250313e-16 9.972491582521764e-15 6 6.661338147750939e-16 7.445216318458507e-14 7 4.345482307321902e-16 7.389206434227901e-14 8 1.110223024625157e-15 6.166454066123704e-13 9 1.311450947838466e-15 1.067657924960871e-12 10 1.299741564375623e-15 4.703851596540448e-12 pyramid_witherden_rule_test(): Normal end of execution. Thu Jul 13 13:46:22 2023