11 July 2023   1:13:42.485 PM

tetrahedron_jaskowiec_rule_test():
  Fortran90 version
  Test tetrahedron_jaskowiec_rule().

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

  Number of nodes N =  14

     I      W          A           B           C           D

     1    0.112688    0.310886    0.310886    0.310886    0.067342
     2    0.112688    0.067342    0.310886    0.310886    0.310886
     3    0.112688    0.310886    0.067342    0.310886    0.310886
     4    0.112688    0.310886    0.310886    0.067342    0.310886
     5    0.073493    0.092735    0.092735    0.092735    0.721794
     6    0.073493    0.721794    0.092735    0.092735    0.092735
     7    0.073493    0.092735    0.721794    0.092735    0.092735
     8    0.073493    0.092735    0.092735    0.721794    0.092735
     9    0.042546    0.045504    0.454496    0.454496    0.045504
    10    0.042546    0.454496    0.045504    0.454496    0.045504
    11    0.042546    0.045504    0.045504    0.454496    0.454496
    12    0.042546    0.454496    0.454496    0.045504    0.045504
    13    0.042546    0.045504    0.454496    0.045504    0.454496
    14    0.042546    0.454496    0.045504    0.045504    0.454496

  Weight Sum =    1.00000    

tetrahedron_jaskowiec_rule_test02():
  Test the precision of a quadrature rule for the unit tetrahedron,

  Number of nodes N =  14

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

  Degree  Maximum error

   0     0.000000000000000    
   1    0.6938893903907228E-17
   2    0.3469446951953614E-17
   3    0.1734723475976807E-17
   4    0.1734723475976807E-17
   5    0.8673617379884035E-18
   6    0.1361833211602648E-04
   7    0.3762506581633752E-04

tetrahedron_jaskowiec_rule_test02():
  Test the precision of a quadrature rule for the unit tetrahedron,
  Check rules of precision p =  0 through 20
  for error in approximating integrals of monomials.

              maximum                   maximum
   p          absolute                  relative
              error                     error

   0     0.000000000000000         0.000000000000000    
   1     0.000000000000000         0.000000000000000    
   2    0.6938893903907228E-17    0.6938893903907228E-17
   3    0.2775557561562891E-16    0.2775557561562891E-16
   4    0.5551115123125783E-16    0.5551115123125783E-16
   5    0.6938893903907228E-17    0.6938893903907228E-17
   6    0.8326672684688674E-16    0.8326672684688674E-16
   7    0.8326672684688674E-16    0.8326672684688674E-16
   8    0.1110223024625157E-15    0.1110223024625157E-15
   9    0.1040834085586084E-16    0.1040834085586084E-16
  10    0.1942890293094024E-15    0.1942890293094024E-15
  11    0.1942890293094024E-15    0.1942890293094024E-15
  12    0.2775557561562891E-16    0.2775557561562891E-16
  13    0.1387778780781446E-15    0.1387778780781446E-15
  14    0.3053113317719180E-15    0.3053113317719180E-15
  15    0.2498001805406602E-15    0.2498001805406602E-15
  16    0.4440892098500626E-15    0.4440892098500626E-15
  17    0.8326672684688674E-16    0.8326672684688674E-16
  18    0.2498001805406602E-15    0.2498001805406602E-15
  19    0.5551115123125783E-16    0.5551115123125783E-16
  20    0.1942890293094024E-15    0.1942890293094024E-15

tetrahedron_jaskowiec_rule_test04():
  Integrate 1/sqrt(r) over the reference tetrahedron.
  Exact integral value is   0.2400589101620030    
  Volume of tetrahedron is   0.1666666666666667    

   P    N     Q     |Q-Exact]

   0    1    0.2532785618838642        0.1321965172186121E-01
   1    1    0.2532785618838642        0.1321965172186121E-01
   2    4    0.2442781387638714        0.4219228601868463E-02
   3    8    0.2421415445769921        0.2082634414989182E-02
   4   14    0.2410193813886114        0.9604712266084448E-03
   5   14    0.2414426895710490        0.1383779409046021E-02
   6   24    0.2403540555991646        0.2951454371616369E-03
   7   35    0.2396527439280942        0.4061662339087113E-03
   8   46    0.2404603979555745        0.4014877935715477E-03
   9   59    0.2402098487285089        0.1509385665059704E-03
  10   81    0.2399908438432391        0.6806631876382641E-04
  11  110    0.2401675962446256        0.1086860826226554E-03
  12  168    0.2400320171196075        0.2689304239547607E-04
  13  172    0.2400773798220388        0.1846966003585249E-04
  14  204    0.2401064662715015        0.4755610949852551E-04
  15  264    0.2400749490947738        0.1603893277082880E-04
  16  304    0.2400527172492860        0.6192912716945775E-05
  17  364    0.2400730008057125        0.1409064370958002E-04
  18  436    0.2400708087778065        0.1189861580358498E-04
  19  487    0.2400533785020640        0.5531659938912137E-05
  20  552    0.2400588444706481        0.6569135482803468E-07

tetrahedron_jaskowiec_rule_test()
  Normal end of execution.

11 July 2023   1:13:42.537 PM