16 July 2016 07:15:15 AM

R85_PRB
  C++ version
  Test the R85 library.

R85_DIF2_TEST
  R85_DIF2 sets up an R85 second difference matrix.

  Matrix order N = 5

  The R85 indicator matrix:

  Col:        1             2             3             4             5       
  Row
  ---
     1           2          -1           0                          
     2          -1           2          -1           0              
     3           0          -1           2          -1           0  
     4                       0          -1           2          -1  
     5                                   0          -1           2  


R85_INDICATOR_TEST
  R85_INDICATOR sets up an R85 indicator matrix.

  Matrix order N = 5

  The R85 indicator matrix:

  Col:        1             2             3             4             5       
  Row
  ---
     1          11          12          13                          
     2          21          22          23          24              
     3          31          32          33          34          35  
     4                      42          43          44          45  
     5                                  53          54          55  


R85_MTV_TEST
  R85_MTV computes b=A'*x where A is an R85 matrix.

  Matrix order N = 5

  The R85 matrix:

  Col:        1             2             3             4             5       
  Row
  ---
     1          11          12          13                          
     2          21          22          23          24              
     3          31          32          33          34          35  
     4                      42          43          44          45  
     5                                  53          54          55  


  The vector x:

     1               1
     2               2
     3               3
     4               4
     5               5

  The product b=A'*x:

     1              11
     2             112
     3             395
     4             508
     5             275

R85_MV_TEST
  R85_MV computes b=A*x where A is an R85 matrix.

  Matrix order N = 5

  The R85 matrix:

  Col:        1             2             3             4             5       
  Row
  ---
     1          11          12          13                          
     2          21          22          23          24              
     3          31          32          33          34          35  
     4                      42          43          44          45  
     5                                  53          54          55  


  The vector x:

     1               1
     2               2
     3               3
     4               4
     5               5

  The product b=A*x:

     1             146
     2             320
     3             595
     4             596
     5             560

R85_NP_FS_TEST
  R85_NP_FS factors and solves a pentadiagonal
  linear system, with no pivoting.

  Matrix order N = 10

  The pentadiagonal matrix:

  Col:        1             2             3             4             5       
  Row
  ---
     1    0.350752    0.043829    0.218418                          
     2    0.351629   0.0945448    0.633966    0.956318              
     3    0.188955    0.822887   0.0136169   0.0617272    0.829509  
     4                0.761492    0.267132    0.859097    0.449539  
     5                            0.396988    0.692066    0.840847  
     6                                        0.185314    0.561662  
     7                                                    0.574366  


  Col:        6             7             8             9            10       
  Row
  ---
     4    0.561695                                                  
     5    0.401306    0.415307                                      
     6    0.123104    0.754673   0.0661187                          
     7    0.861216  0.00751236    0.797287    0.257578              
     8    0.367027    0.453794    0.260303  0.00183837    0.109957  
     9                0.617205    0.911977    0.912484    0.897504  
    10                            0.361529    0.597917    0.113664  


  Right hand side:

     1         1.62087
     2         5.74755
     3         4.58067
     4         10.1064
     5         15.8815
     6         13.9567
     7         15.8424
     8         19.8832
     9         16.0093
    10         10.0938

  Solution to A*x=b:

     1               1
     2               2
     3               3
     4               4
     5               5
     6               6
     7               7
     8               8
     9               9
    10              10

R85_PRINT_TEST
  R85_PRINT prints an R85 matrix.

  Matrix order N = 5

  The R85 matrix:

  Col:        1             2             3             4             5       
  Row
  ---
     1          11          12          13                          
     2          21          22          23          24              
     3          31          32          33          34          35  
     4                      42          43          44          45  
     5                                  53          54          55  


R85_PRINT_SOME_TEST
  R85_PRINT_SOME prints some of an R85 matrix.

  Matrix order N = 8

  Rows 1-5, Cols 2-4:

  Col:        2             3             4       
  Row
  ---
     1          12          13              
     2          22          23          24  
     3          32          33          34  
     4          42          43          44  
     5                      53          54  


R85_RANDOM_TEST
  R85_RANDOM randomizes an R85 matrix.

  Matrix order N = 5

  The R85 matrix:

  Col:        1             2             3             4             5       
  Row
  ---
     1    0.109957    0.561695    0.218418                          
     2    0.401306    0.043829    0.415307    0.956318              
     3    0.897504    0.754673    0.633966   0.0661187    0.829509  
     4                0.350752    0.797287   0.0617272    0.257578  
     5                           0.0945448  0.00183837    0.449539  


R85_TO_R8GE_TEST
  R85_TO_R8GE converts an R85 matrix to R8GE format.

  Matrix order N = 5

  The R85 matrix:

  Col:        1             2             3             4             5       
  Row
  ---
     1          11          12          13                          
     2          21          22          23          24              
     3          31          32          33          34          35  
     4                      42          43          44          45  
     5                                  53          54          55  


  The R8GE matrix:

  Col:          1             2             3             4             5       
  Row
  ---
    1            11            21            31             0             0  
    2            12            22            32            42             0  
    3            13            23            33            43            53  
    4             0            24            34            44            54  
    5             0             0            35            45            55  

R85_ZEROS_TEST
  R85_ZEROS zeros an R85 indicator matrix.

  Matrix order N = 5

  The R85 zero matrix:

  Col:        1             2             3             4             5       
  Row
  ---
     1           0           0           0                          
     2           0           0           0           0              
     3           0           0           0           0           0  
     4                       0           0           0           0  
     5                                   0           0           0  


R85_PRB
  Normal end of execution.

16 July 2016 07:15:15 AM