8 June 2020   8:43:31.210 AM
 
QR_SOLVE_TEST
  FORTRAN90 version
  Test the QR_SOLVE library.
  QR_SOLVE needs the R8LIB library.
  This test also needs the TEST_LS library.
 
NORMAL_SOLVE_TEST
  NORMAL_SOLVE is a function with a simple interface which
  solves a linear system A*x = b in the least squares sense.
  Compare a tabulated solution X1 to the NORMAL_SOLVE result X2.
 
  NORMAL_SOLVE cannot be applied when N < M,
  or if the matrix does not have full column rank.
 
  Number of problems =    6
 
  Index     M     N     ||B||         ||X1 - X2||   ||X1||       ||X2||        ||R1||        ||R2||
 
      1     5     3     6.205        0.2861E-07     5.461         5.461         1.075         1.075    
      2     6     3     537.6        0.5416E-04     194.8         194.8         5.222         5.222    
      3     5     3     40.37      ------------     10.84      ------------    0.5477E-06  ------------
      4     3     5     232.1      ------------     7.416      ------------     0.000      ------------
      5    10    10     1.000      ------------    0.1534E+08  ------------    0.9673E-10  ------------
      6    10    10     1.000        0.5261E-15     1.000         1.000        0.4958E-15    0.1788E-15
 
QR_SOLVE_TEST
  QR_SOLVE is a function with a simple interface which
  solves a linear system A*x = b in the least squares sense.
  Compare a tabulated solution X1 to the QR_SOLVE result X2.
 
  Number of problems =    6
 
  Index     M     N     ||B||         ||X1 - X2||   ||X1||       ||X2||        ||R1||        ||R2||
 
      1     5     3     6.205        0.2861E-07     5.461         5.461         1.075         1.075    
      2     6     3     537.6        0.5416E-04     194.8         194.8         5.222         5.222    
      3     5     3     40.37         26.68         10.84         28.80        0.5477E-06    0.4936E-13
      4     3     5     232.1         10.37         7.416         12.75         0.000        0.1097E-12
      5    10    10     1.000         2182.        0.1534E+08    0.1534E+08    0.9673E-10    0.5318E-09
      6    10    10     1.000        0.5752E-15     1.000         1.000        0.4958E-15    0.3419E-15
 
SVD_SOLVE_TEST
  SVD_SOLVE is a function with a simple interface which
  solves a linear system A*x = b in the least squares sense.
  Compare a tabulated solution X1 to the SVD_SOLVE result X2.
 
  Number of problems =    6
 
  Index     M     N     ||B||         ||X1 - X2||   ||X1||       ||X2||        ||R1||        ||R2||
 
      1     5     3     6.205        0.2861E-07     5.461         5.461         1.075         1.075    
      2     6     3     537.6        0.5416E-04     194.8         194.8         5.222         5.222    
      3     5     3     40.37        0.5092E-07     10.84         10.84        0.5477E-06    0.4082E-13
      4     3     5     232.1        0.1462E-12     7.416         7.416         0.000        0.1015E-12
      5    10    10     1.000         1731.        0.1534E+08    0.1534E+08    0.9673E-10    0.5091E-09
      6    10    10     1.000        0.1165E-14     1.000         1.000        0.4958E-15    0.1077E-14
 
DQRLS_TEST
  DQRLS solves a linear system A*x = b 
  in the least squares sense.
 
  Coefficient matrix A:
 
  Col          1             2             3      
  Row
 
    1:      1.            1.            1.      
    2:      1.            2.            4.      
    3:      1.            3.            9.      
    4:      1.            4.           16.      
    5:      1.            5.           25.      
 
  Right hand side b:
 
         1:    1.0000000    
         2:    2.3000000    
         3:    4.6000000    
         4:    3.1000000    
         5:    1.2000000    
 
  Error code =   0
  Estimated matrix rank =   3
 
  Least squares solution x:
 
         1:   -3.0200000    
         2:    4.4914286    
         3:  -0.72857143    
 
  Residuals A*x-b
 
         1:   0.25714286    
         2:  -0.74857143    
         3:   0.70285714    
         4:  -0.18857143    
         5:  -0.22857143E-01
 
QR_SOLVE_TEST
  Normal end of execution.
 
 8 June 2020   8:43:31.211 AM