12 July 2023   3:54:47.306 PM

pyramid_jaskowiec_rule_test():
  Fortran90 version
  Test pyramid_jaskowiec_rule().

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

  Number of nodes N =  15

     I      W          X           Y           Z

     1    0.067734    0.000000    0.000000    0.730709
     2    0.064709    0.000000    0.000000    0.006197
     3    0.177272    0.000000    0.000000    0.268446
     4    0.059108    0.753441    0.000000    0.125000
     5    0.059108    0.000000    0.753441    0.125000
     6    0.059108   -0.753441    0.000000    0.125000
     7    0.059108    0.000000   -0.753441    0.125000
     8    0.065375    0.417152    0.417152    0.421822
     9    0.065375   -0.417152    0.417152    0.421822
    10    0.065375    0.417152   -0.417152    0.421822
    11    0.065375   -0.417152   -0.417152    0.421822
    12    0.048088    0.674023    0.674023    0.065796
    13    0.048088   -0.674023    0.674023    0.065796
    14    0.048088    0.674023   -0.674023    0.065796
    15    0.048088   -0.674023   -0.674023    0.065796

  Weight Sum =    1.00000    

pyramid_jaskowiec_rule_test02():
  Test the precision of a quadrature rule for the unit pyramid,

  Number of nodes N =  15

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

  Degree  Maximum error

   0    0.2220446049250313E-15
   1    0.2775557561562891E-15
   2    0.2775557561562891E-16
   3    0.8326672684688674E-16
   4    0.1196959198423997E-15
   5    0.3469446951953614E-16
   6    0.8772414128785586E-02
   7    0.6056836049176022E-03

pyramid_jaskowiec_rule_test03():
  Test the precision of quadrature rules for the unit pyramid,
  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.1665334536937735E-15    0.4996003610813202E-15
   2    0.2220446049250313E-15    0.6661338147750936E-15
   3    0.1665334536937735E-15    0.1665334536937738E-14
   4    0.2220446049250313E-15    0.9699272635055414E-14
   5    0.2775557561562891E-15    0.9426053687589064E-14
   6    0.2324529457808922E-15    0.7322267792098640E-13
   7    0.4371503159461554E-15    0.7389206434227901E-13
   8    0.8881784197001252E-15    0.6125398040357083E-12
   9    0.1303211011327576E-14    0.1064840354565121E-11
  10    0.1303644692196571E-14    0.4705835635694111E-11
  11    0.1455866677213535E-14    0.5078631505385482E-11
  12    0.5471317843230850E-14    0.1782910999959874E-10
  13    0.1311385895708117E-13    0.3555919840680948E-10
  14    0.1311212423360519E-13    0.2409720516017026E-09
  15    0.1621391822886897E-13    0.8912399120581262E-09
  16    0.1621391822886897E-13    0.1501765970065100E-08
  17    0.7491294910788593E-13    0.4892707874773605E-08
  18    0.1361616304733013E-12    0.2317872039290490E-07
  19    0.1361616169207741E-12    0.7986468657124887E-07
  20    0.1544540726870422E-12    0.1478863312362768E-06

pyramid_jaskowiec_rule_test()
  Normal end of execution.

12 July 2023   3:54:47.355 PM