Wed Jul 12 20:30:50 2023

pyramid_jaskowiec_rule_test():
  Python version: 3.8.10
  Test pyramid_jaskowiec_rule().

pyramid_jaskowiec_rule_test01():
  Quadrature rule for the unit pyramid,
  Precision p =  5

     I      W          X           Y           Z

     0    0.067734    0.000000    0.000000    0.730709
     1    0.064709    0.000000    0.000000    0.006197
     2    0.177272    0.000000    0.000000    0.268446
     3    0.059108    0.753441    0.000000    0.125000
     4    0.059108    0.000000    0.753441    0.125000
     5    0.059108   -0.753441    0.000000    0.125000
     6    0.059108    0.000000   -0.753441    0.125000
     7    0.065375    0.417152    0.417152    0.421822
     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.048088    0.674023    0.674023    0.065796
    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

  Weight Sum    1.0

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

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

  Degree  Maximum error

   0     2.220446049250313e-16
   1     2.775557561562891e-16
   2     2.775557561562891e-17
   3     1.249000902703301e-16
   4     1.196959198423997e-16
   5     3.469446951953614e-17
   6      0.008772414128785593
   7     0.0006056836049176022

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                         0
   1     1.665334536937735e-16     4.996003610813202e-16
   2     2.220446049250313e-16     6.661338147750936e-16
   3     1.665334536937735e-16     1.665334536937738e-15
   4     2.220446049250313e-16     9.699272635055414e-15
   5     2.775557561562891e-16     9.426053687589064e-15
   6     2.315855840429037e-16     7.294945897352004e-14
   7     4.380176776841438e-16     7.389206434227901e-14
   8      7.07116656895046e-16     6.125398040357083e-13
   9     1.303211011327576e-15     1.065229447619772e-12
  10     1.304945734803553e-15     4.705225162108368e-12
  11     1.456734038951524e-15     5.078292353393403e-12
  12     5.464378949326942e-15     1.782945339099073e-11
  13     1.311342527621218e-14      3.55595313924017e-11
  14     1.311299159534318e-14     2.409713164897596e-10
  15     1.621370138843448e-14     8.912396433407355e-10
  16     1.621348454799998e-14     1.501765729988434e-09
  17     7.491262384723418e-14      4.89270635762246e-09
  18     1.361616846834099e-13      2.31787203929049e-08
  19     1.361616711308827e-13     7.986468636387497e-08
  20     1.544540523582515e-13     1.478863315684453e-07

pyramid_jaskowiec_rule_test():
  Normal end of execution.
Wed Jul 12 20:30:51 2023