22 June 2015   9:07:53.608 PM
 
PYRAMID_INTEGRALS_PRB
  FORTRAN90 version
  Test the PYRAMID_INTEGRALS library.
 
TEST01
  Compare exact and estimated integrals 
  over the unit pyramid in 3D.

  Number of sample points used is 500000

   E1  E2  E3     MC-Estimate      Exact           Error

   0   0   0     1.33333         1.33333        0.00E+00
   2   0   0    0.265432        0.266667        0.12E-02
   4   0   0    0.113616        0.114286        0.67E-03
   6   0   0    0.630826E-01    0.634921E-01    0.41E-03
   0   2   0    0.266404        0.266667        0.26E-03
   2   2   0    0.631172E-01    0.634921E-01    0.37E-03
   4   2   0    0.294246E-01    0.296296E-01    0.20E-03
   0   4   0    0.114151        0.114286        0.13E-03
   2   4   0    0.294431E-01    0.296296E-01    0.19E-03
   0   6   0    0.633977E-01    0.634921E-01    0.94E-04
   0   0   1    0.333908        0.333333        0.58E-03
   2   0   1    0.443021E-01    0.444444E-01    0.14E-03
   4   0   1    0.142281E-01    0.142857E-01    0.58E-04
   0   2   1    0.445419E-01    0.444444E-01    0.97E-04
   2   2   1    0.791511E-02    0.793651E-02    0.21E-04
   0   4   1    0.143255E-01    0.142857E-01    0.40E-04
   0   0   2    0.133741        0.133333        0.41E-03
   2   0   2    0.126683E-01    0.126984E-01    0.30E-04
   4   0   2    0.316495E-02    0.317460E-02    0.97E-05
   0   2   2    0.127463E-01    0.126984E-01    0.48E-04
   2   2   2    0.176164E-02    0.176367E-02    0.20E-05
   0   4   2    0.318874E-02    0.317460E-02    0.14E-04
   0   0   3    0.669290E-01    0.666667E-01    0.26E-03
   2   0   3    0.475339E-02    0.476190E-02    0.85E-05
   0   2   3    0.478520E-02    0.476190E-02    0.23E-04
   0   0   4    0.382648E-01    0.380952E-01    0.17E-03
   2   0   4    0.211345E-02    0.211640E-02    0.29E-05
   0   2   4    0.212875E-02    0.211640E-02    0.12E-04
   0   0   5    0.239229E-01    0.238095E-01    0.11E-03
   0   0   6    0.159526E-01    0.158730E-01    0.80E-04
 
TEST02
  Print sample points in the unit pyramid in 3D.
 
  Unit pyramid points
 
  Row          1             2             3      
  Col
 
    1: -0.529145      0.857507      0.604056E-01
    2:  0.498951E-01 -0.684939E-01  0.595634    
    3: -0.170945     -0.275040      0.647423    
    4:  0.205248     -0.671476      0.233952    
    5: -0.183031      0.472301      0.727328E-01
    6: -0.702641      0.560667      0.294766    
    7: -0.770880     -0.924746      0.493645E-01
    8:  0.133511     -0.147631      0.804149    
    9: -0.232220      0.399618      0.515596    
   10: -0.191112      0.415901      0.355966    
   11:  0.365470      0.117333      0.485835E-01
   12: -0.778534E-01  0.694145      0.157545    
   13: -0.457207      0.384369      0.265047    
   14: -0.450616      0.106489      0.284023    
   15:  0.139976     -0.165375      0.402856    
   16:  0.286934     -0.511457      0.331066    
   17:  0.257074      0.256803      0.604505    
   18:  0.387721     -0.757334      0.802734E-01
   19:  0.473565     -0.369359      0.347348E-01
   20:  0.390594      0.387917      0.258939    
 
PYRAMID_INTEGRALS_PRB
  Normal end of execution.
 
22 June 2015   9:07:53.965 PM