Tue Oct 19 17:01:47 2021

pyramid_integrals_test():
  Python version: 3.6.9
  Test pyramid_integrals().

r8_mop_test():
  r8_mop() evaluates (-1.0)^I4 as an R8.

    I4  r8_mop(I4)

   -81  -1.0
   -95  -1.0
     5  -1.0
   -20   1.0
   -85  -1.0
    26   1.0
   -13  -1.0
    22   1.0
   -64   1.0
   -76   1.0

r8_mop_test():
  Normal end of execution.

pyramid01_monomial_integral_test
  Python version: 3.6.9
  pyramid01_monomial_integral returns the integral of a monomial
  over the unit pyramid in 3D.
  Compare to a Monte Carlo estimate.

  Number of sample points used is 500000

   E1  E2  E3     MC-Estimate      Exact           Error

   0   0   0         1.33333         1.33333           0
   2   0   0        0.266579        0.266667     8.8e-05
   4   0   0        0.114297        0.114286     1.2e-05
   6   0   0        0.063527       0.0634921     3.5e-05
   0   2   0        0.266346        0.266667     0.00032
   2   2   0        0.063447       0.0634921     4.5e-05
   4   2   0       0.0296463       0.0296296     1.7e-05
   0   4   0        0.114028        0.114286     0.00026
   2   4   0       0.0295906       0.0296296     3.9e-05
   0   6   0       0.0632972       0.0634921     0.00019
   0   0   1        0.333682        0.333333     0.00035
   2   0   1        0.044477       0.0444444     3.3e-05
   4   0   1       0.0142891       0.0142857     3.4e-06
   0   2   1       0.0444765       0.0444444     3.2e-05
   2   2   1      0.00794202      0.00793651     5.5e-06
   0   4   1       0.0142914       0.0142857     5.7e-06
   0   0   2        0.133545        0.133333     0.00021
   2   0   2       0.0127213       0.0126984     2.3e-05
   4   0   2      0.00317931       0.0031746     4.7e-06
   0   2   2       0.0127106       0.0126984     1.2e-05
   2   2   2      0.00176739      0.00176367     3.7e-06
   0   4   2      0.00317745       0.0031746     2.9e-06
   0   0   3        0.066795       0.0666667     0.00013
   2   0   3      0.00477375       0.0047619     1.2e-05
   0   2   3      0.00476631       0.0047619     4.4e-06
   0   0   4       0.0381765       0.0380952     8.1e-05
   2   0   4      0.00212258       0.0021164     6.2e-06
   0   2   4      0.00211862       0.0021164     2.2e-06
   0   0   5       0.0238636       0.0238095     5.4e-05
   0   0   6       0.0159108        0.015873     3.8e-05

pyramid01_monomial_integral_test:
  Normal end of execution.

pyramid01_sample_test()
  Python version: 3.6.9
  pyramid01_sample samples points from the unit pyramid.

  Unit pyramid points

  Row: 
      0       
      1       
      2       

  Col
      0 :
   -0.674756  
   -0.648947  
   0.0742593  

      1 :
   -0.243269  
   0.0475062  
    0.289185  

      2 :
    0.833337  
    0.316699  
    0.115197  

      3 :
   -0.841984  
    0.875385  
   0.0491976  

      4 :
   -0.596778  
   -0.224387  
    0.268963  

      5 :
   -0.370739  
   -0.885588  
  0.00294414  

      6 :
   -0.383839  
   0.0572241  
    0.232698  

      7 :
    0.438783  
  0.00330122  
  0.00259256  

      8 :
   -0.411389  
   -0.757258  
    0.117798  

      9 :
   -0.128211  
   -0.785722  
    0.139722  

     10 :
   -0.484444  
   -0.675923  
     0.24624  

     11 :
   -0.220869  
   -0.206541  
    0.560238  

     12 :
   -0.777504  
   -0.687117  
   0.0767844  

     13 :
   -0.247341  
    0.154715  
     0.29743  

     14 :
    0.316319  
   -0.270306  
    0.581244  

     15 :
   0.0456574  
   -0.574791  
     0.35614  

     16 :
    0.304278  
   -0.537551  
    0.400158  

     17 :
    0.522088  
   -0.161692  
   0.0865056  

     18 :
    0.545416  
   -0.308364  
   0.0454193  

     19 :
    0.585218  
    0.869256  
   0.0931414  


pyramid01_sample_test
  Normal end of execution.

pyramid01_volume_test()
  Python version: 3.6.9
  pyramid01_volume() returns the volume of the unit pyramid.

  pyramid01_volume() = 1.33333

pyramid01_volume_test
  Normal end of execution.

pyramid_integrals_test():
  Normal end of execution.
Tue Oct 19 17:02:02 2021