13 July 2023 07:07:13 PM

prism_jaskowiec_rule_test():
  C version
  Test prism_jaskowiec_rule().

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

  Number of nodes N = 16

     I      W          X           Y           Z

     0    0.207143    0.333333    0.333333         0.5
     1   0.0380756   0.0517646   0.0517646         0.5
     2   0.0380756   0.0517646    0.896471         0.5
     3   0.0380756    0.896471   0.0517646         0.5
     4   0.0763743    0.166397    0.166397    0.903582
     5   0.0763743    0.166397    0.667206    0.903582
     6   0.0763743    0.667206    0.166397    0.903582
     7   0.0763743    0.166397    0.166397   0.0964183
     8   0.0763743    0.166397    0.667206   0.0964183
     9   0.0763743    0.667206    0.166397   0.0964183
    10   0.0367308    0.497665    0.497665    0.301369
    11   0.0367308    0.497665  0.00467002    0.301369
    12   0.0367308  0.00467002    0.497665    0.301369
    13   0.0367308    0.497665    0.497665    0.698631
    14   0.0367308    0.497665  0.00467002    0.698631
    15   0.0367308  0.00467002    0.497665    0.698631

  Weight Sum = 1

prism_jaskowiec_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                         0
   1     2.775557561562891e-17
   2     5.551115123125783e-17
   3     8.326672684688674e-17
   4     2.775557561562891e-17
   5     4.163336342344337e-17
   6     0.0001192533706372984
   7     0.0004173867972305859

prism_jaskowiec_rule_test03():
  Test the precision of quadrature rules for the unit prism.
  Check rules of precision p = 0 through 20
  for error in approximating integrals of monomials.

              maximum                   maximum
   p          absolute                  relative
              error                     error

   0                         0                         0
   1                         0                         0
   2     1.110223024625157e-16     2.220446049250313e-16
   3     5.551115123125783e-17     2.775557561562891e-16
   4     5.551115123125783e-17     4.683753385137379e-16
   5     8.326672684688674e-17     6.661338147750939e-16
   6      3.33066907387547e-16     1.554312234475219e-15
   7     1.387778780781446e-16     1.554312234475219e-15
   8     3.885780586188048e-16     1.092875789865388e-15
   9     4.996003610813204e-16     1.405126015541214e-15
  10      3.33066907387547e-16     3.042976027406441e-15
  11     2.775557561562891e-16     3.278627369596165e-15
  12     3.608224830031759e-16     7.310124727766265e-15
  13      3.33066907387547e-16      4.08006961549745e-15
  14     7.771561172376096e-16     9.742207041085749e-15
  15     8.881784197001252e-16     7.077671781985373e-15
  16     1.221245327087672e-15     2.604687299179176e-14
  17     5.551115123125783e-16     1.987472686426628e-14
  18     4.996003610813204e-15     2.979387569990166e-14
  19      6.38378239159465e-16     1.845508609214092e-14
  20     4.440892098500626e-16     4.431351119382754e-14

prism_jaskowiec_rule_test()
  Normal end of execution.

13 July 2023 07:07:13 PM