14 July 2023   1:40:39.481 PM

prism_witherden_rule_test():
  Fortran90 version
  Test prism_witherden_rule().

prism_witherden_rule_test01():
  Quadrature rule for the unit prism,
  Precision p =  5

  Number of nodes N =  16

     I      W          X           Y           Z

     1    0.177864        0.333333        0.333333        0.500000    
     2    0.561775E-01    0.487300        0.254001E-01    0.500000    
     3    0.561775E-01    0.254001E-01    0.487300        0.500000    
     4    0.561775E-01    0.487300        0.487300        0.500000    
     5    0.464154E-01    0.445598        0.108804        0.644985E-01
     6    0.464154E-01    0.108804        0.445598        0.644985E-01
     7    0.464154E-01    0.445598        0.445598        0.644985E-01
     8    0.464154E-01    0.445598        0.108804        0.935501    
     9    0.464154E-01    0.108804        0.445598        0.935501    
    10    0.464154E-01    0.445598        0.445598        0.935501    
    11    0.625186E-01    0.100859        0.798282        0.214787    
    12    0.625186E-01    0.798282        0.100859        0.214787    
    13    0.625186E-01    0.100859        0.100859        0.214787    
    14    0.625186E-01    0.100859        0.798282        0.785213    
    15    0.625186E-01    0.798282        0.100859        0.785213    
    16    0.625186E-01    0.100859        0.100859        0.785213    

  Weight Sum    1.00000    

prism_witherden_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.5551115123125783E-16
   1    0.2775557561562891E-16
   2    0.1387778780781446E-16
   3    0.2775557561562891E-16
   4    0.5551115123125783E-16
   5    0.4163336342344337E-16
   6    0.8681433560950649E-04
   7    0.2425567824553866E-03

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

              maximum                   maximum
   p          absolute                  relative
              error                     error

   0     0.000000000000000         0.000000000000000    
   1    0.5551115123125783E-16    0.3330669073875470E-15
   2    0.2775557561562891E-16    0.1665334536937735E-15
   3    0.1110223024625157E-15    0.2220446049250313E-15
   4    0.5551115123125783E-16    0.3330669073875470E-15
   5    0.5551115123125783E-16    0.6245004513516506E-15
   6    0.6106226635438361E-15    0.2137179322403426E-14
   7    0.1387778780781446E-15    0.1578598363138894E-14
   8    0.1665334536937735E-15    0.1040834085586084E-14
   9    0.1665334536937735E-15    0.1748601263784622E-14
  10    0.3330669073875470E-15    0.2629732369363591E-14

prism_witherden_rule_test():
  Normal end of execution.

14 July 2023   1:40:39.482 PM