14 July 2023 03:58:35 PM

quadrilateral_witherden_rule_test():
  C version
  Test quadrilateral_witherden_rule().

quadrilateral_witherden_rule_test01():
  Quadrature rule for the unit quadrilateral.
  Precision p = 5

  Number of nodes N = 8

     I      W          X           Y

     0    0.204082    0.841565         0.5
     1    0.204082         0.5    0.841565
     2    0.204082    0.158435         0.5
     3    0.204082         0.5    0.158435
     4   0.0459184    0.940959    0.940959
     5   0.0459184    0.940959   0.0590414
     6   0.0459184   0.0590414    0.940959
     7   0.0459184   0.0590414   0.0590414

  Weight Sum = 1

quadrilateral_witherden_rule_test02():
  Test the precision of a quadrature rule for the unit quadrilateral.

  Number of nodes N = 8

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

  Degree  Maximum error

   0     1.110223024625157e-16
   1     1.110223024625157e-16
   2     5.551115123125783e-17
   3     2.775557561562891e-17
   4     2.775557561562891e-17
   5     2.775557561562891e-17
   6     0.0003086419753086378
   7      0.000817901234567861

quadrilateral_witherden_rule_test03():
  Test the precision of quadrature rules for the unit quadrilateral.
  Check rules of precision p = 0 through' 21
  for error in approximating integrals of monomials.

              maximum                   maximum
   p          absolute                  relative
              error                     error

   0                         0                         0
   1                         0                         0
   2     5.551115123125783e-17     1.665334536937735e-16
   3     5.551115123125783e-17     1.665334536937735e-16
   4     1.110223024625157e-16     2.220446049250313e-16
   5     1.110223024625157e-16     2.220446049250313e-16
   6     5.551115123125783e-17     2.081668171172169e-16
   7     5.551115123125783e-17     2.081668171172169e-16
   8     2.220446049250313e-16     8.673617379884035e-16
   9     2.220446049250313e-16     1.040834085586084e-15
  10     4.440892098500626e-16     1.124100812432971e-15
  11     4.440892098500626e-16     1.249000902703301e-15
  12     6.661338147750939e-16     1.998401444325282e-15
  13     6.661338147750939e-16     1.998401444325282e-15
  14     2.775557561562891e-16     1.332267629550188e-15
  15     2.775557561562891e-16     1.332267629550188e-15
  16     6.661338147750939e-16     1.360023205165817e-15
  17     6.661338147750939e-16     1.360023205165817e-15
  18     6.661338147750939e-16     5.495603971894524e-15
  19     6.661338147750939e-16     5.533767888366015e-15
  20     1.110223024625157e-15     2.810252031082427e-15
  21     1.110223024625157e-15     2.810252031082427e-15

quadrilateral_witherden_rule_test()
  Normal end of execution.

14 July 2023 03:58:35 PM