17 July 2023 11:26:57 AM

hexahedron_jaskowiec_rule_test():
  C version
  Test hexahedron_jaskowiec_rule().

hexahedron_jaskowiec_rule_test01():
  Quadrature rule for the unit hexahedron.
  Precision p = 5

  Number of nodes N = 14

     I      W          X           Y           Z

     0    0.110803    0.102089         0.5         0.5
     1    0.110803         0.5    0.102089         0.5
     2    0.110803         0.5         0.5    0.102089
     3    0.110803    0.897911         0.5         0.5
     4    0.110803         0.5    0.897911         0.5
     5    0.110803         0.5         0.5    0.897911
     6   0.0418975    0.879393    0.879393    0.879393
     7   0.0418975    0.120607    0.120607    0.120607
     8   0.0418975    0.120607    0.879393    0.879393
     9   0.0418975    0.879393    0.120607    0.879393
    10   0.0418975    0.879393    0.879393    0.120607
    11   0.0418975    0.120607    0.120607    0.879393
    12   0.0418975    0.879393    0.120607    0.120607
    13   0.0418975    0.120607    0.879393    0.120607

  Weight Sum = 1

hexahedron_jaskowiec_rule_test02():
  Test the precision of a quadrature rule for the unit hexahedron.

  Number of nodes N = 14

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

  Degree  Maximum error

   0     2.220446049250313e-16
   1      3.33066907387547e-16
   2     1.110223024625157e-16
   3     1.665334536937735e-16
   4     1.110223024625157e-16
   5     1.110223024625157e-16
   6     0.0004208754208754398
   7      0.001235269360269312

hexahedron_jaskowiec_rule_test03():
  Test the precision of quadrature rules for the unit hexahedron.
  Check rules of precision p = 0 through 21
  for error in approximating integrals of monomials.

              maximum                   maximum
   p          absolute                  relative
              error                     error

   0                         0                         0
   1                         0                         0
   2     2.220446049250313e-16      3.33066907387547e-16
   3     2.220446049250313e-16      3.33066907387547e-16
   4      3.33066907387547e-16     9.992007221626409e-16
   5      3.33066907387547e-16     9.992007221626409e-16
   6      3.33066907387547e-16     7.771561172376096e-16
   7      3.33066907387547e-16     7.771561172376096e-16
   8      3.33066907387547e-16     9.992007221626409e-16
   9      3.33066907387547e-16     9.992007221626409e-16
  10     8.326672684688674e-16     2.498001805406602e-15
  11     8.326672684688674e-16     2.498001805406602e-15
  12     1.554312234475219e-15     5.620504062164855e-15
  13     1.554312234475219e-15     5.634381849972669e-15
  14     9.992007221626409e-16     5.027228633380787e-15
  15     9.992007221626409e-16     5.027228633380787e-15
  16     2.498001805406602e-15     1.348920974919565e-14
  17     2.498001805406602e-15     1.515454428613339e-14
  18     1.110223024625157e-15     9.908740494779522e-15
  19     1.110223024625157e-15     1.015854067532018e-14
  20     8.770761894538737e-15     3.844190590851504e-14
  21     8.770761894538737e-15     4.174438572590589e-14

hexahedron_jaskowiec_rule_test()
  Normal end of execution.

17 July 2023 11:26:57 AM