15 July 2023 03:03:46 PM

square_symq_rule_test():
  C version
  Test square_symq_rule().

test01():
  Symmetric quadrature rule for a square.
  precision = 5

  Number of nodes N = 7

     J  W       X       Y

     0       0.0794152        0.508879       0.0170357
     1       0.0794152        0.491121        0.982964
     2         0.13795        0.889436        0.214215
     3         0.13795        0.110564        0.785785
     4        0.139778        0.114811        0.208516
     5        0.139778        0.885189        0.791484
     6        0.285714             0.5             0.5
  Weight sum  1

test02():
  Test a rule of precision 5.

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

  Degree  Maximum error

   0                         0
   1                         0
   2     1.110223024625157e-16
   3     5.551115123125783e-17
   4     5.551115123125783e-17
   5     8.326672684688674e-17
   6     0.0004624560744359296
   7      0.001248225890155455

test02():
  Test the precision of a quadrature rule for the unit square.
  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                         0                         0
   2     5.551115123125783e-17     1.665334536937735e-16
   3     5.551115123125783e-17     2.220446049250313e-16
   4     5.551115123125783e-17     2.775557561562891e-16
   5     1.110223024625157e-16     4.996003610813204e-16
   6     1.110223024625157e-16      3.33066907387547e-16
   7     1.110223024625157e-16     4.163336342344337e-16
   8     4.440892098500626e-16     4.440892098500626e-16
   9     4.440892098500626e-16     4.440892098500626e-16
  10     1.110223024625157e-16     4.996003610813204e-16
  11     2.220446049250313e-16      7.28583859910259e-16
  12     4.440892098500626e-16     5.828670879282072e-16
  13     4.440892098500626e-16     5.100087019371813e-16
  14      3.33066907387547e-16     1.040834085586084e-15
  15     1.110223024625157e-16     9.020562075079397e-16
  16     6.661338147750939e-16     1.165734175856414e-15
  17     2.220446049250313e-16     1.873501354054952e-15
  18     1.665334536937735e-16     8.586881206085196e-16
  19     1.110223024625157e-16     1.498801083243961e-15
  20     5.551115123125783e-16     2.914335439641036e-15

square_symq_rule_test():
  Normal end of execution.

15 July 2023 03:03:46 PM