February 15 2008   3:11:25.141 PM
 
EISPACK_PRB2
  Do some symmetric eigenproblem tests.
 
TEST01
  RS computes the eigenvalues and eigenvectors
  of a real symmetric matrix.
 
  Matrix order =        4
 
  Initialize random number generator using SEED =        12345
  RANDOM_INITIALIZE returns suggested seed =        12345
 
  SYMM_TEST will give us a symmetric matrix
  with known eigenstructure.
 
  The matrix A:
 
  Col          1             2             3             4      
  Row
 
    1    1.92644      0.316021     -0.267767      0.178532E-01
    2   0.316021       1.17753      0.156898     -0.169315    
    3  -0.267767      0.156898      0.495902       1.05229    
    4   0.178532E-01 -0.169315       1.05229      0.138188E-01
 
 
  The eigenvector matrix Q:
 
  Col          1             2             3             4      
  Row
 
    1  -0.921626      0.785585E-01 -0.102911      0.365846    
    2  -0.306868     -0.123536     -0.394218     -0.857416    
    3   0.204536      0.622736     -0.734700      0.174870    
    4   0.120840     -0.768614     -0.542417      0.316882    
 
 
  LAMBDA_MIN =  -0.867793    
  LAMBDA_MAX =    2.08874    
 
  The eigenvalue vector LAMBDA:
 
         1     2.0887432    
         2   -0.86779261    
         3     1.3194709    
         4     1.0732642    
 
  The column norms of A*Q:
 
         1     2.0887432    
         2    0.86779261    
         3     1.3194709    
         4     1.0732642    
 
  Now call EISPACK routine RS
  and see if it can recover Q and LAMBDA.
 
  LAMBDA_MIN =  -0.867793    
  LAMBDA_MAX =    2.08874    
 
  The computed eigenvalues Lambda:
 
         1   -0.86779261    
         2     1.0732642    
         3     1.3194709    
         4     2.0887432    
 
  The eigenvector matrix:
 
  Col          1             2             3             4      
  Row
 
    1   0.785585E-01  0.365846      0.102911     -0.921626    
    2  -0.123536     -0.857416      0.394218     -0.306868    
    3   0.622736      0.174870      0.734700      0.204536    
    4  -0.768614      0.316882      0.542417      0.120840    
 
 
  The residual (A-Lambda*I)*X:
 
  Col          1             2             3             4      
  Row
 
    1   0.416334E-16  0.222045E-15 -0.555112E-16 -0.444089E-15
    2  -0.277556E-16 -0.222045E-15  0.333067E-15  0.222045E-15
    3   0.111022E-15  0.166533E-15  0.555112E-15  0.166533E-15
    4  -0.333067E-15 -0.277556E-15  0.666134E-15 -0.555112E-16
 
 
  Setup time =   0.820000E-04
  Solve time =   0.470000E-04
 
TEST01
  RS computes the eigenvalues and eigenvectors
  of a real symmetric matrix.
 
  Matrix order =       16
 
  Initialize random number generator using SEED =    573802814
  RANDOM_INITIALIZE returns suggested seed =    573802814
 
  SYMM_TEST will give us a symmetric matrix
  with known eigenstructure.
 
  LAMBDA_MIN =  -0.423239    
  LAMBDA_MAX =    3.58584    
 
  Now call EISPACK routine RS
  and see if it can recover Q and LAMBDA.
 
  LAMBDA_MIN =  -0.423239    
  LAMBDA_MAX =    3.58584    
 
  Setup time =   0.432000E-03
  Solve time =   0.433000E-03
 
TEST01
  RS computes the eigenvalues and eigenvectors
  of a real symmetric matrix.
 
  Matrix order =       64
 
  Initialize random number generator using SEED =   1447861724
  RANDOM_INITIALIZE returns suggested seed =   1447861724
 
  SYMM_TEST will give us a symmetric matrix
  with known eigenstructure.
 
  LAMBDA_MIN =   -1.68551    
  LAMBDA_MAX =    3.11144    
 
  Now call EISPACK routine RS
  and see if it can recover Q and LAMBDA.
 
  LAMBDA_MIN =   -1.68551    
  LAMBDA_MAX =    3.11144    
 
  Setup time =   0.215230E-01
  Solve time =   0.158160E-01
 
TEST01
  RS computes the eigenvalues and eigenvectors
  of a real symmetric matrix.
 
  Matrix order =      256
 
  Initialize random number generator using SEED =   1983473816
  RANDOM_INITIALIZE returns suggested seed =   1983473816
 
  SYMM_TEST will give us a symmetric matrix
  with known eigenstructure.
 
  LAMBDA_MIN =   -1.16179    
  LAMBDA_MAX =    4.37245    
 
  Now call EISPACK routine RS
  and see if it can recover Q and LAMBDA.
 
  LAMBDA_MIN =   -1.16179    
  LAMBDA_MAX =    4.37245    
 
  Setup time =    1.80524    
  Solve time =   0.848499    
 
EISPACK_PRB2
  Normal end of execution.
 
February 15 2008   3:11:27.841 PM