11 July      2016  10:26:57.164 PM      
 
LINPLUS_PRB
  FORTRAN77 version:
  Test the LINPLUS library.
 
TEST01
  C83_CR_FA factors a complex tridiagonal matrix;
  C83_CR_SLS solves 1 or more factored systems.
 
  Matrix order N =       10
 
  Solution:
 
  Col:          1          
  Row
  ---
    1:  1.00      10.0    
    2:  2.00      20.0    
    3:  3.00      30.0    
    4:  4.00      40.0    
    5:  5.00      50.0    
    6:  6.00      60.0    
    7:  7.00      70.0    
    8:  8.00      80.0    
    9:  9.00      90.0    
   10:  10.0      100.    
 
TEST0196
  For a general matrix,
  R8VEC_TO_R8GE converts a real vector to an R8GE matrix.
  R8GE_TO_R8VEC converts an R8GE matrix to a real vector.
 
  Matrix rows M =           4
  Matrix columns N =        6
 
  The R8GE indicator matrix:
 
  Col:        1             2             3             4             5       
  Row
  ---
    1    11.0000       12.0000       13.0000       14.0000       15.0000    
    2    21.0000       22.0000       23.0000       24.0000       25.0000    
    3    31.0000       32.0000       33.0000       34.0000       35.0000    
    4    41.0000       42.0000       43.0000       44.0000       45.0000    
 
  Col:        6       
  Row
  ---
    1    16.0000    
    2    26.0000    
    3    36.0000    
    4    46.0000    
       1       1       1   11.0000    
       2       1       2   21.0000    
       3       1       3   31.0000    
       4       1       4   41.0000    
       1       2       5   12.0000    
       2       2       6   22.0000    
       3       2       7   32.0000    
       4       2       8   42.0000    
       1       3       9   13.0000    
       2       3      10   23.0000    
       3       3      11   33.0000    
       4       3      12   43.0000    
       1       4      13   14.0000    
       2       4      14   24.0000    
       3       4      15   34.0000    
       4       4      16   44.0000    
       1       5      17   15.0000    
       2       5      18   25.0000    
       3       5      19   35.0000    
       4       5      20   45.0000    
       1       6      21   16.0000    
       2       6      22   26.0000    
       3       6      23   36.0000    
       4       6      24   46.0000    
 
  The recovered R8GE indicator matrix:
 
  Col:        1             2             3             4             5       
  Row
  ---
    1    11.0000       12.0000       13.0000       14.0000       15.0000    
    2    21.0000       22.0000       23.0000       24.0000       25.0000    
    3    31.0000       32.0000       33.0000       34.0000       35.0000    
    4    41.0000       42.0000       43.0000       44.0000       45.0000    
 
  Col:        6       
  Row
  ---
    1    16.0000    
    2    26.0000    
    3    36.0000    
    4    46.0000    
 
R8BTO_SL_TEST
  R8BTO_SL solves a block Toeplitz system.
 
  Block order M =         2
  Block number L =        3
  Matrix order N =        6
 
  The block Toeplitz matrix:
 
  Col:       1             2             3             4             5       
  Row
  ---
    1       9.            1.            3.            4.            5.      
    2       2.            8.            6.            6.            7.      
    3       7.            8.            9.            1.            3.      
    4       8.            8.            2.            8.            6.      
    5       9.            0.            7.            8.            9.      
    6       9.            9.            8.            8.            2.      
 
  Col:       6       
  Row
  ---
    1       6.      
    2       7.      
    3       4.      
    4       6.      
    5       1.      
    6       8.      
 
  The right hand side B:
 
         1     97.000000    
         2     137.00000    
         3     93.000000    
         4     128.00000    
         5     113.00000    
         6     141.00000    
 
  The computed solution X:
 
         1     1.0000000    
         2     2.0000000    
         3     3.0000000    
         4     4.0000000    
         5     5.0000000    
         6     6.0000000    
 
R8GE_DET_TEST
  R8GE_DET computes the determinant
  of an R8GE matrix.
 
  Matrix order N =        4
 
  R8GE_DET computes the determinant =    112.000    
  Exact determinant =                    112.000    
 
TEST295
  For a matrix in general storage,
  R8GE_DILU returns the DILU factors.
 
  Matrix rows M =           9
  Matrix columns N =        9
 
  Matrix A:
 
  Col:        1             2             3             4             5       
  Row
  ---
    1    4.00000      -1.00000       0.00000      -1.00000       0.00000    
    2   -1.00000       4.00000      -1.00000       0.00000      -1.00000    
    3    0.00000      -1.00000       4.00000      -1.00000       0.00000    
    4   -1.00000       0.00000      -1.00000       4.00000      -1.00000    
    5    0.00000      -1.00000       0.00000      -1.00000       4.00000    
    6    0.00000       0.00000      -1.00000       0.00000      -1.00000    
    7    0.00000       0.00000       0.00000      -1.00000       0.00000    
    8    0.00000       0.00000       0.00000       0.00000      -1.00000    
    9    0.00000       0.00000       0.00000       0.00000       0.00000    
 
  Col:        6             7             8             9       
  Row
  ---
    1    0.00000       0.00000       0.00000       0.00000    
    2    0.00000       0.00000       0.00000       0.00000    
    3   -1.00000       0.00000       0.00000       0.00000    
    4    0.00000      -1.00000       0.00000       0.00000    
    5   -1.00000       0.00000      -1.00000       0.00000    
    6    4.00000      -1.00000       0.00000      -1.00000    
    7   -1.00000       4.00000      -1.00000       0.00000    
    8    0.00000      -1.00000       4.00000      -1.00000    
    9   -1.00000       0.00000      -1.00000       4.00000    
 
  DILU factor:
 
         1    0.25000000    
         2    0.26666667    
         3    0.26785714    
         4    0.28717949    
         5    0.29017857    
         6    0.29053178    
         7    0.29220211    
         8    0.29260134    
         9    0.29266578    
 
TEST31
  For a matrix in general storage,
  R8GE_FA computes the LU factors,
  R8GE_SL solves a factored system.
 
  Matrix order N =        5
 
  The matrix:
 
  Col:        1             2             3             4             5       
  Row
  ---
    1   0.218418      0.661187E-01  0.617272E-01  0.183837E-02  0.859097    
    2   0.956318      0.257578      0.449539      0.897504      0.840847    
    3   0.829509      0.109957      0.401306      0.350752      0.123104    
    4   0.561695      0.438290E-01  0.754673      0.945448E-01  0.751236E-02
    5   0.415307      0.633966      0.797287      0.136169E-01  0.260303    
 
  The compressed LU factors:
 
  Col         1             2             3             4             5       
  Row
 
    1   0.956318      0.257578      0.449539      0.897504      0.840847    
    2  -0.228395      0.522106      0.602062     -0.376149     -0.104858    
    3  -0.867399      0.217324      0.614552     -0.510026     -0.507943    
    4  -0.587352      0.205820     -0.231419     -0.391459     -0.511487    
    5  -0.434277     -0.139612E-01  0.803036E-01 -0.610161      0.939815    
 
  The pivot vector P:
 
         1             2
         2             5
         3             4
         4             4
         5             5
 
  Solution:
 
         1     1.0000000    
         2     2.0000000    
         3     3.0000000    
         4     4.0000000    
         5     5.0000000    
 
TEST315
  For a matrix in general storage,
  R8GE_ILU returns the ILU factors.
 
  Matrix rows M =           9
  Matrix columns N =        9
 
  Matrix A:
 
  Col:        1             2             3             4             5       
  Row
  ---
    1    4.00000      -1.00000       0.00000      -1.00000       0.00000    
    2   -1.00000       4.00000      -1.00000       0.00000      -1.00000    
    3    0.00000      -1.00000       4.00000      -1.00000       0.00000    
    4   -1.00000       0.00000      -1.00000       4.00000      -1.00000    
    5    0.00000      -1.00000       0.00000      -1.00000       4.00000    
    6    0.00000       0.00000      -1.00000       0.00000      -1.00000    
    7    0.00000       0.00000       0.00000      -1.00000       0.00000    
    8    0.00000       0.00000       0.00000       0.00000      -1.00000    
    9    0.00000       0.00000       0.00000       0.00000       0.00000    
 
  Col:        6             7             8             9       
  Row
  ---
    1    0.00000       0.00000       0.00000       0.00000    
    2    0.00000       0.00000       0.00000       0.00000    
    3   -1.00000       0.00000       0.00000       0.00000    
    4    0.00000      -1.00000       0.00000       0.00000    
    5   -1.00000       0.00000      -1.00000       0.00000    
    6    4.00000      -1.00000       0.00000      -1.00000    
    7   -1.00000       4.00000      -1.00000       0.00000    
    8    0.00000      -1.00000       4.00000      -1.00000    
    9   -1.00000       0.00000      -1.00000       4.00000    
 
  Factor L:
 
  Col:        1             2             3             4             5       
  Row
  ---
    1    1.00000       0.00000       0.00000       0.00000       0.00000    
    2  -0.250000       1.00000       0.00000       0.00000       0.00000    
    3    0.00000     -0.266667       1.00000       0.00000       0.00000    
    4  -0.250000       0.00000     -0.267857       1.00000       0.00000    
    5    0.00000     -0.266667       0.00000     -0.287179       1.00000    
    6    0.00000       0.00000     -0.267857       0.00000     -0.290179    
    7    0.00000       0.00000       0.00000     -0.287179       0.00000    
    8    0.00000       0.00000       0.00000       0.00000     -0.290179    
    9    0.00000       0.00000       0.00000       0.00000       0.00000    
 
  Col:        6             7             8             9       
  Row
  ---
    1    0.00000       0.00000       0.00000       0.00000    
    2    0.00000       0.00000       0.00000       0.00000    
    3    0.00000       0.00000       0.00000       0.00000    
    4    0.00000       0.00000       0.00000       0.00000    
    5    0.00000       0.00000       0.00000       0.00000    
    6    1.00000       0.00000       0.00000       0.00000    
    7  -0.290532       1.00000       0.00000       0.00000    
    8    0.00000     -0.292202       1.00000       0.00000    
    9  -0.290532       0.00000     -0.292601       1.00000    
 
  Factor U:
 
  Col:        1             2             3             4             5       
  Row
  ---
    1    4.00000      -1.00000       0.00000      -1.00000       0.00000    
    2    0.00000       3.75000      -1.00000       0.00000      -1.00000    
    3    0.00000       0.00000       3.73333      -1.00000       0.00000    
    4    0.00000       0.00000       0.00000       3.48214      -1.00000    
    5    0.00000       0.00000       0.00000       0.00000       3.44615    
    6    0.00000       0.00000       0.00000       0.00000       0.00000    
    7    0.00000       0.00000       0.00000       0.00000       0.00000    
    8    0.00000       0.00000       0.00000       0.00000       0.00000    
    9    0.00000       0.00000       0.00000       0.00000       0.00000    
 
  Col:        6             7             8             9       
  Row
  ---
    1    0.00000       0.00000       0.00000       0.00000    
    2    0.00000       0.00000       0.00000       0.00000    
    3   -1.00000       0.00000       0.00000       0.00000    
    4    0.00000      -1.00000       0.00000       0.00000    
    5   -1.00000       0.00000      -1.00000       0.00000    
    6    3.44196      -1.00000       0.00000      -1.00000    
    7    0.00000       3.42229      -1.00000       0.00000    
    8    0.00000       0.00000       3.41762      -1.00000    
    9    0.00000       0.00000       0.00000       3.41687    
 
  Product L*U:
 
  Col:        1             2             3             4             5       
  Row
  ---
    1    4.00000      -1.00000       0.00000      -1.00000       0.00000    
    2   -1.00000       4.00000      -1.00000      0.250000      -1.00000    
    3    0.00000      -1.00000       4.00000      -1.00000      0.266667    
    4   -1.00000      0.250000      -1.00000       4.00000      -1.00000    
    5    0.00000      -1.00000      0.266667      -1.00000       4.00000    
    6    0.00000       0.00000      -1.00000      0.267857      -1.00000    
    7    0.00000       0.00000       0.00000      -1.00000      0.287179    
    8    0.00000       0.00000       0.00000       0.00000      -1.00000    
    9    0.00000       0.00000       0.00000       0.00000       0.00000    
 
  Col:        6             7             8             9       
  Row
  ---
    1    0.00000       0.00000       0.00000       0.00000    
    2    0.00000       0.00000       0.00000       0.00000    
    3   -1.00000       0.00000       0.00000       0.00000    
    4   0.267857      -1.00000       0.00000       0.00000    
    5   -1.00000      0.287179      -1.00000       0.00000    
    6    4.00000      -1.00000      0.290179      -1.00000    
    7   -1.00000       4.00000      -1.00000      0.290532    
    8   0.290179      -1.00000       4.00000      -1.00000    
    9   -1.00000      0.290532      -1.00000       4.00000    
 
TEST317
  For a matrix in general storage,
  R8GE_INDICATOR sets up the indicator matrix.
 
  Matrix rows M =           7
  Matrix columns N =        5
 
  The R8GE indicator matrix:
 
  Col:        1             2             3             4             5       
  Row
  ---
    1    11.0000       12.0000       13.0000       14.0000       15.0000    
    2    21.0000       22.0000       23.0000       24.0000       25.0000    
    3    31.0000       32.0000       33.0000       34.0000       35.0000    
    4    41.0000       42.0000       43.0000       44.0000       45.0000    
    5    51.0000       52.0000       53.0000       54.0000       55.0000    
    6    61.0000       62.0000       63.0000       64.0000       65.0000    
    7    71.0000       72.0000       73.0000       74.0000       75.0000    
 
TEST34
  For a matrix in general storage,
  R8GE_FS factors and solves a linear system.
 
  Matrix order N =       10
 
  Solution:
 
         1     1.0000000    
         2     2.0000000    
         3     3.0000000    
         4     4.0000000    
         5     5.0000000    
         6     6.0000000    
         7     7.0000000    
         8     8.0000000    
         9     9.0000000    
        10     10.000000    
 
TEST345
  For a matrix in general storage,
  R8GE_FSS factors and solves multiple linear system.
 
  Matrix order N =       10
 
  Solutions:
 
  Col         1             2             3       
  Row
 
    1    1.00000       1.00000       1.00000    
    2    1.00000       2.00000       2.00000    
    3    1.00000       3.00000       3.00000    
    4    1.00000       4.00000       1.00000    
    5    1.00000       5.00000       2.00000    
    6    1.00000       6.00000       3.00000    
    7    1.00000       7.00000       1.00000    
    8    1.00000       8.00000       2.00000    
    9    1.00000       9.00000       3.00000    
   10    1.00000       10.0000       1.00000    
 
LINPLUS_PRB
  Normal end of execution.
 
11 July      2016  10:26:57.167 PM