17 July 2023 12:29:37 PM

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.22045e-16
   1               3.33067e-16
   2               1.11022e-16
   3               1.66533e-16
   4               1.11022e-16
   5               1.11022e-16
   6               0.000420875
   7                0.00123527

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.22045e-16               3.33067e-16
   3               2.22045e-16               3.33067e-16
   4               3.33067e-16               9.99201e-16
   5               3.33067e-16               9.99201e-16
   6               3.33067e-16               7.77156e-16
   7               3.33067e-16               7.77156e-16
   8               3.33067e-16               9.99201e-16
   9               3.33067e-16               9.99201e-16
  10               8.32667e-16                 2.498e-15
  11               8.32667e-16                 2.498e-15
  12               1.55431e-15                5.6205e-15
  13               1.55431e-15               5.63438e-15
  14               9.99201e-16               5.02723e-15
  15               9.99201e-16               5.02723e-15
  16                 2.498e-15               1.34892e-14
  17                 2.498e-15               1.51545e-14
  18               1.11022e-15               9.90874e-15
  19               1.11022e-15               1.01585e-14
  20               8.77076e-15               3.84419e-14
  21               8.77076e-15               4.17444e-14

hexahedron_jaskowiec_rule_test()
  Normal end of execution.

17 July 2023 12:29:38 PM