16 July 2023 07:29:36 AM

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():
  Get a quadrature rule for the symmetric square.
  Test its accuracy.
  Precision = 5

  Number of quadrature points = 7

  Degree  Maximum error

   0                         0
   1                         0
   2               1.11022e-16
   3               5.55112e-17
   4               5.55112e-17
   5               8.32667e-17
   6               0.000462456
   7                0.00124823

test03():
  Test the precision of quadrature rules for the unit quadrilateral.
  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.55112e-17               1.66533e-16
   3               5.55112e-17               2.22045e-16
   4               5.55112e-17               2.77556e-16
   5               1.11022e-16                 4.996e-16
   6               1.11022e-16               3.33067e-16
   7               1.11022e-16               4.16334e-16
   8               4.44089e-16               4.44089e-16
   9               4.44089e-16               4.44089e-16
  10               1.11022e-16                 4.996e-16
  11               2.22045e-16               7.28584e-16
  12               4.44089e-16               5.82867e-16
  13               4.44089e-16               5.10009e-16
  14               3.33067e-16               1.04083e-15
  15               1.11022e-16               9.02056e-16
  16               6.66134e-16               1.16573e-15
  17               2.22045e-16                1.8735e-15
  18               1.66533e-16               8.58688e-16
  19               1.11022e-16                1.4988e-15
  20               5.55112e-16               2.91434e-15

square_symq_rule_test():
  Normal end of execution.

16 July 2023 07:29:37 AM