21 August    2023   8:51:47.566 AM      

asa053_test():
  FORTRAN77 version
  Test asa053().

TEST01
  Generate a single Wishart deviate.

  The number of variables is    3
  The number of degrees of freedom is    1
 
  The upper Cholesky factor:
 
  Col:       1             2             3       
  Row
  ---
    1    3.00000       2.00000       1.00000    
    2    0.00000       4.00000       2.00000    
    3    0.00000       0.00000       5.00000    
 
  The sample matrix:
 
  Col:       1             2             3       
  Row
  ---
    1    54.5541       43.2451      -3.81605    
    2    43.2451       53.5239     -0.155780    
    3   -3.81605     -0.155780       4.72559    

TEST02
  Average many Wishart samples.
  Compare to D' * D * np / n.

  The number of variables is    3
  The number of degrees of freedom is    2
 
  The upper Cholesky factor:
 
  Col:       1             2             3       
  Row
  ---
    1    3.00000       2.00000       1.00000    
    2    0.00000       4.00000       2.00000    
    3    0.00000       0.00000       5.00000    
 
  The averaged matrix:
 
  Col:       1             2             3       
  Row
  ---
    1    13.4591       8.97139       4.48088    
    2    8.97139       30.0573       14.9865    
    3    4.48088       14.9865       45.0242    
 
  Expected result:
 
  Col         1             2             3       
  Row
 
    1:   13.5000       9.00000       4.50000    
    2:   9.00000       30.0000       15.0000    
    3:   4.50000       15.0000       45.0000    

asa053_test():
  Normal end of execution.

21 August    2023   8:51:47.612 AM