22 July      2014   9:26:21.987 PM      
 
ST_IO_PRB:
  FORTRAN77 version
  Test the ST_IO library.
 
TEST01
  ST_WRITE writes an ST file.
 
  Sparse Triplet (ST) header information:
 
  Minimum row index I_MIN =        0
  Maximum row index I_MAX =        4
  Minimum col index J_MIN =        0
  Maximum col index J_MAX =        4
  Number of rows        M =        5
  Number of columns     N =        5
  Number of nonzeros  NST =       11
 
  Sparse Triplet (ST) data:
     #     I     J       A
  ----  ----  ----  --------------
 
     1     4     0     51.000000    
     2     0     1     12.000000    
     3     0     0     11.000000    
     4     2     2     33.000000    
     5     0     4     15.000000    
     6     4     2     53.000000    
     7     4     4     55.000000    
     8     1     1     22.000000    
     9     2     4     35.000000    
    10     3     3     44.000000    
    11     1     0     21.000000    
 
  Wrote the matrix data to "a5by5.st".
 
TEST02
  ST_HEADER_READ reads the header from an ST file.
  ST_DATA_READ reads the data from an ST file.
 
  Read the data from "kershaw.st".
 
  Sparse Triplet (ST) header information:
 
  Minimum row index I_MIN =        0
  Maximum row index I_MAX =        3
  Minimum col index J_MIN =        0
  Maximum col index J_MAX =        3
  Number of rows        M =        4
  Number of columns     N =        4
  Number of nonzeros  NST =       12
 
  Sparse Triplet (ST) data read from file
     #     I     J       A
  ----  ----  ----  --------------
 
     1     0     0     3.0000000    
     2     1     0    -2.0000000    
     3     3     0     2.0000000    
     4     0     1    -2.0000000    
     5     1     1     3.0000000    
     6     2     1    -2.0000000    
     7     1     2    -2.0000000    
     8     2     2     3.0000000    
     9     3     2    -2.0000000    
    10     0     3     2.0000000    
    11     2     3    -2.0000000    
    12     3     3     3.0000000    
 
TEST03
  ST_SORT_A sorts an ST matrix by columns.
 
  Sparse Triplet (ST) header information:
 
  Minimum row index I_MIN =        1
  Maximum row index I_MAX =        5
  Minimum col index J_MIN =        1
  Maximum col index J_MAX =        5
  Number of rows        M =        5
  Number of columns     N =        5
  Number of nonzeros  NST =       11
 
  Matrix data before sorting:
     #     I     J       A
  ----  ----  ----  --------------
 
     1     5     1     51.000000    
     2     1     2     12.000000    
     3     1     1     11.000000    
     4     3     3     33.000000    
     5     1     5     15.000000    
     6     5     3     53.000000    
     7     5     5     55.000000    
     8     2     2     22.000000    
     9     3     5     35.000000    
    10     4     4     44.000000    
    11     2     1     21.000000    
 
  Matrix data sorted by column:
     #     I     J       A
  ----  ----  ----  --------------
 
     1     1     1     11.000000    
     2     2     1     21.000000    
     3     5     1     51.000000    
     4     1     2     12.000000    
     5     2     2     22.000000    
     6     3     3     33.000000    
     7     5     3     53.000000    
     8     4     4     44.000000    
     9     1     5     15.000000    
    10     3     5     35.000000    
    11     5     5     55.000000    
 
ST_IO_PRB:
  Normal end of execution.
 
22 July      2014   9:26:22.012 PM