14 July 2023 10:01:51 AM

prism_witherden_rule_test():
  C version
  Test prism_witherden_rule().

prism_witherden_rule_test01():
  Quadrature rule for the unit prism.
  Precision p = 5

  Number of nodes N = 16

     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 = 1

prism_witherden_rule_test02():
  Test the precision of a quadrature rule for the unit prism.

  Number of nodes N = 16

  Stated precision of rule    = 5
  Number of quadrature points = 16

  Degree  Maximum error

   0     5.551115123125783e-17
   1     2.775557561562891e-17
   2     1.387778780781446e-17
   3     2.775557561562891e-17
   4     5.551115123125783e-17
   5     4.163336342344337e-17
   6     8.681433560950649e-05
   7     0.0002425567824553884

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     5.551115123125783e-17     6.245004513516506e-16
   6     6.106226635438361e-16     2.137179322403426e-15
   7     1.387778780781446e-16     1.554312234475219e-15
   8     1.665334536937735e-16     9.992007221626409e-16
   9     1.665334536937735e-16     1.942890293094024e-15
  10      3.33066907387547e-16     2.629732369363591e-15

prism_witherden_rule_test()
  Normal end of execution.

14 July 2023 10:01:51 AM