10 July 2023 08:24:11 PM

triangle_symq_rule_test():
  C++ version
  Test triangle_symq_rule().

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

  Number of nodes N = 7

     I      W          A           B           C

     0    0.125939    0.101287    0.101287    0.797427
     1    0.125939    0.797427    0.101287    0.101287
     2    0.125939    0.101287    0.797427    0.101287
     3    0.132394    0.470142   0.0597159    0.470142
     4    0.132394    0.470142    0.470142   0.0597159
     5    0.132394   0.0597159    0.470142    0.470142
     6       0.225    0.333333    0.333333    0.333333

  Weight Sum = 1

triangle_symq_rule_test02():
  Test the precision of a quadrature rule for the unit triangle.

  Number of nodes N = 7

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

  Degree  Maximum error

   0               1.11022e-16
   1               5.55112e-17
   2               5.55112e-17
   3               3.46945e-17
   4               2.77556e-17
   5               2.77556e-17
   6               8.18846e-05
   7               0.000254052

triangle_symq_rule_test03():
  Test the precision of a quadrature rule for the unit triangle.
  Check rules of precision p = 0 through 50
  for error in approximating integrals of monomials.

              maximum                   maximum
   p          absolute                  relative
              error                     error

   0                         0                         0
   1               5.55112e-17               3.33067e-16
   2               2.77556e-17               6.66134e-16
   3               1.11022e-16               1.04083e-15
   4               1.11022e-16               1.45717e-15
   5               1.11022e-16               1.63931e-15
   6               1.11022e-16                1.7486e-15
   7               1.11022e-16               1.96718e-15
   8               5.55112e-17               2.81025e-15
   9               1.11022e-16               3.22008e-15
  10               1.11022e-16               2.63331e-15
  11               1.11022e-16               3.16283e-15
  12               8.32667e-16               4.42008e-15
  13               5.55112e-17               4.29181e-15
  14               8.32667e-17                6.7394e-15
  15               1.66533e-16               5.30825e-15
  16               2.22045e-16                6.2372e-15
  17               2.22045e-16                4.8659e-15
  18               3.33067e-16               6.30355e-15
  19               2.22045e-16               8.19657e-15
  20               1.94289e-16               8.81586e-15
  21               8.32667e-17               1.00806e-14
  22               4.44089e-16               1.03687e-14
  23               2.22045e-16               9.88792e-15
  24               2.22045e-16               1.14166e-14
  25               1.66533e-16               1.52222e-14
  26               1.66533e-16               2.98355e-14
  27               1.66533e-16               1.60448e-14
  28               6.66134e-16               8.48723e-15
  29               7.77156e-16               2.13761e-14
  30               4.44089e-16               2.71033e-14
  31               1.11022e-16               3.02258e-14
  32               6.66134e-16               1.11437e-14
  33               5.55112e-17               2.60621e-14
  34               1.83187e-15               4.01632e-14
  35               6.66134e-16               1.94284e-14
  36               1.94289e-16               2.94207e-14
  37               1.66533e-15               9.43215e-15
  38               1.66533e-15               2.45247e-14
  39               3.33067e-16               6.27666e-14
  40               3.33067e-16               2.89384e-14
  41               5.55112e-16               2.28299e-14
  42               2.22045e-16               2.15388e-14
  43               1.66533e-15               1.74735e-14
  44               9.99201e-16               8.77521e-14
  45               1.94289e-16               3.87317e-14
  46               8.88178e-16               4.17036e-14
  47               2.22045e-16               7.03812e-14
  48               4.44089e-16               1.44635e-13
  49               8.88178e-16                1.1045e-13
  50               3.33067e-16               3.86728e-14

triangle_symq_rule_test()
  Normal end of execution.

10 July 2023 08:24:12 PM