10 July 2023 06:06:57 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        0.0000000000000001
   1        0.0000000000000001
   2        0.0000000000000001
   3        0.0000000000000000
   4        0.0000000000000000
   5        0.0000000000000000
   6        0.0000818846056941
   7        0.0002540522381792

triangle_symq_rule_test03():
  Test the precision of quadrature rules 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.551115123125783e-17      3.33066907387547e-16
   2     2.775557561562891e-17     6.661338147750939e-16
   3     1.110223024625157e-16     1.040834085586084e-15
   4     1.110223024625157e-16     1.457167719820518e-15
   5     1.110223024625157e-16     1.639313684798083e-15
   6     1.110223024625157e-16     1.748601263784622e-15
   7     1.110223024625157e-16     1.967176421757699e-15
   8     5.551115123125783e-17     2.810252031082428e-15
   9     1.110223024625157e-16     3.220080452281949e-15
  10     1.110223024625157e-16     2.633310236532793e-15
  11     1.110223024625157e-16     3.162834577574714e-15
  12     8.326672684688674e-16     4.420075416788904e-15
  13     5.551115123125783e-17     4.291814299783869e-15
  14     8.326672684688674e-17     6.739400704169896e-15
  15     1.665334536937735e-16      5.30825383648903e-15
  16     2.220446049250313e-16      6.23719825787461e-15
  17     2.220446049250313e-16     4.865899350114944e-15
  18      3.33066907387547e-16     6.303551430830723e-15
  19     2.220446049250313e-16     8.196568423990414e-15
  20     1.942890293094024e-16     8.815864704914134e-15
  21     8.326672684688674e-17      1.00806407492271e-14
  22     4.440892098500626e-16     1.036865905634787e-14
  23     2.220446049250313e-16       9.8879238130678e-15
  24     2.220446049250313e-16     1.141664887627236e-14
  25     1.665334536937735e-16     1.522219850169648e-14
  26     1.665334536937735e-16     2.983550906332511e-14
  27     1.665334536937735e-16     1.604478268996123e-14
  28     6.661338147750939e-16     8.487230426818688e-15
  29     7.771561172376096e-16     2.137613003272421e-14
  30     4.440892098500626e-16     2.710331958866163e-14
  31     1.110223024625157e-16     3.022582184541989e-14
  32     6.661338147750939e-16     1.114369728019744e-14
  33     5.551115123125783e-17     2.606205182220656e-14
  34     1.831867990631508e-15     4.016318527755303e-14
  35     6.661338147750939e-16     1.942842859248978e-14
  36     1.942890293094024e-16     2.942069331213215e-14
  37     1.665334536937735e-15     9.432152324809207e-15
  38     1.665334536937735e-15     2.452465314162211e-14
  39      3.33066907387547e-16     6.276663216953082e-14
  40      3.33066907387547e-16      2.89384401858106e-14
  41     5.551115123125783e-16     2.282986218691066e-14
  42     2.220446049250313e-16     2.153876035859703e-14
  43     1.665334536937735e-15      1.74735395400305e-14
  44     9.992007221626409e-16     8.775207123445926e-14
  45     1.942890293094024e-16     3.873168265603431e-14
  46     8.881784197001252e-16     4.170361972422043e-14
  47     2.220446049250313e-16     7.038120086733102e-14
  48     4.440892098500626e-16     1.446352803149975e-13
  49     8.881784197001252e-16     1.104503858165296e-13
  50      3.33066907387547e-16     3.867284097125445e-14

triangle_symq_rule_test()
  Normal end of execution.

10 July 2023 06:06:57 PM