17 July 2023   1:13:56.455 PM

hexahedron_jaskowiec_rule_test():
  Fortran90 version
  Test hexahedron_jaskowiec_rule().

hexahedron_jaskowiec_rule_test01():
  Quadrature rule for the hexahedron,
  given in barycentric coordinates.
  Precision p =  5

  Number of nodes N =  14

     I      W          X           Y           Z

     1    0.110803    0.102089    0.500000    0.500000
     2    0.110803    0.500000    0.102089    0.500000
     3    0.110803    0.500000    0.500000    0.102089
     4    0.110803    0.897911    0.500000    0.500000
     5    0.110803    0.500000    0.897911    0.500000
     6    0.110803    0.500000    0.500000    0.897911
     7    0.041898    0.879393    0.879393    0.879393
     8    0.041898    0.120607    0.120607    0.120607
     9    0.041898    0.120607    0.879393    0.879393
    10    0.041898    0.879393    0.120607    0.879393
    11    0.041898    0.879393    0.879393    0.120607
    12    0.041898    0.120607    0.120607    0.879393
    13    0.041898    0.879393    0.120607    0.120607
    14    0.041898    0.120607    0.879393    0.120607

  Weight Sum =    1.00000    

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    0.2220446049250313E-15
   1    0.3330669073875470E-15
   2    0.1110223024625157E-15
   3    0.1665334536937735E-15
   4    0.1110223024625157E-15
   5    0.1110223024625157E-15
   6    0.4208754208754398E-03
   7    0.1235269360269312E-02

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.000000000000000         0.000000000000000    
   1     0.000000000000000         0.000000000000000    
   2    0.2220446049250313E-15    0.3330669073875470E-15
   3    0.2220446049250313E-15    0.3330669073875470E-15
   4    0.3330669073875470E-15    0.9992007221626409E-15
   5    0.3330669073875470E-15    0.9992007221626409E-15
   6    0.3330669073875470E-15    0.7771561172376096E-15
   7    0.3330669073875470E-15    0.7771561172376096E-15
   8    0.3330669073875470E-15    0.9992007221626409E-15
   9    0.3330669073875470E-15    0.9992007221626409E-15
  10    0.8326672684688674E-15    0.2498001805406602E-14
  11    0.8326672684688674E-15    0.2498001805406602E-14
  12    0.1554312234475219E-14    0.5204170427930421E-14
  13    0.1554312234475219E-14    0.5245803791353865E-14
  14    0.9992007221626409E-15    0.5027228633380787E-14
  15    0.9992007221626409E-15    0.5245803791353865E-14
  16    0.2498001805406602E-14    0.1360023205165817E-13
  17    0.2498001805406602E-14    0.1548761119352093E-13
  18    0.1110223024625157E-14    0.9617306950815419E-14
  19    0.1110223024625157E-14    0.9835882108788496E-14
  20    0.8770761894538737E-14    0.3863315917174148E-13
  21    0.8770761894538737E-14    0.4186321428401030E-13

hexahedron_jaskowiec_rule_test()
  Normal end of execution.

17 July 2023   1:13:56.523 PM