13 July 2023   9:39:05.930 PM

prism_jaskowiec_rule_test():
  Fortran90 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

     1    0.207143        0.333333        0.333333        0.500000    
     2    0.380756E-01    0.517646E-01    0.517646E-01    0.500000    
     3    0.380756E-01    0.517646E-01    0.896471        0.500000    
     4    0.380756E-01    0.896471        0.517646E-01    0.500000    
     5    0.763743E-01    0.166397        0.166397        0.903582    
     6    0.763743E-01    0.166397        0.667206        0.903582    
     7    0.763743E-01    0.667206        0.166397        0.903582    
     8    0.763743E-01    0.166397        0.166397        0.964183E-01
     9    0.763743E-01    0.166397        0.667206        0.964183E-01
    10    0.763743E-01    0.667206        0.166397        0.964183E-01
    11    0.367308E-01    0.497665        0.497665        0.301369    
    12    0.367308E-01    0.497665        0.467002E-02    0.301369    
    13    0.367308E-01    0.467002E-02    0.497665        0.301369    
    14    0.367308E-01    0.497665        0.497665        0.698631    
    15    0.367308E-01    0.497665        0.467002E-02    0.698631    
    16    0.367308E-01    0.467002E-02    0.497665        0.698631    

  Weight Sum    1.00000    

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.000000000000000    
   1    0.2775557561562891E-16
   2    0.5551115123125783E-16
   3    0.8326672684688674E-16
   4    0.4163336342344337E-16
   5    0.4163336342344337E-16
   6    0.1192533706373122E-03
   7    0.4173867972306067E-03

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.000000000000000         0.000000000000000    
   1     0.000000000000000         0.000000000000000    
   2    0.1110223024625157E-15    0.2220446049250313E-15
   3    0.5551115123125783E-16    0.2775557561562891E-15
   4    0.5551115123125783E-16    0.4683753385137379E-15
   5    0.8326672684688674E-16    0.6661338147750939E-15
   6    0.3330669073875470E-15    0.1554312234475219E-14
   7    0.1387778780781446E-15    0.1578598363138894E-14
   8    0.3885780586188048E-15    0.1405126015541214E-14
   9    0.4996003610813204E-15    0.1366094737331736E-14
  10    0.3330669073875470E-15    0.2592164764086968E-14
  11    0.2775557561562891E-15    0.3278627369596165E-14
  12    0.3608224830031759E-15    0.7310124727766265E-14
  13    0.3330669073875470E-15    0.3770421475035589E-14
  14    0.7771561172376096E-15    0.9742207041085749E-14
  15    0.8881784197001252E-15    0.7313594174718219E-14
  16    0.1221245327087672E-14    0.2557329348285009E-13
  17    0.5551115123125783E-15    0.1987472686426628E-13
  18    0.4996003610813204E-14    0.2979387569990166E-13
  19    0.6383782391594650E-15    0.1859596461192825E-13
  20    0.4440892098500626E-15    0.4414437565491980E-13

prism_jaskowiec_rule_test():
  Normal end of execution.

13 July 2023   9:39:05.977 PM