27-Aug-2022 16:29:07

r8sr_test():
  MATLAB version
  Test r8sr().

R8SR_DIF2_TEST
  R8SR_DIF2 sets up an R8SR second difference matrix;

  Matrix order N = 5

  The R8SR second difference matrix:

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

R8SR_INDICATOR_TEST
  R8SR_INDICATOR sets up an R8SR indicator matrix;

  Matrix order N = 5

  The R8SR indicator matrix:

  Col:       1             2             3             4             5       
  Row
  ---
   1            11            12             0             0            15  
   2             0            22             0             0            25  
   3            31             0            33             0             0  
   4            41             0             0            44             0  
   5             0            52            53             0            55  

R8SR_MTV_TEST
  R8SR_MTV multiplies a vector by an R8SR matrix;

  Matrix order N = 5

  The R8SR matrix, in R8GE form:

  Col:       1             2             3             4             5       
  Row
  ---
   1      0.218418      0.956318             0             0      0.829509  
   2             0      0.561695             0             0      0.415307  
   3      0.257578             0     0.0661187             0             0  
   4      0.043829             0             0      0.109957             0  
   5             0     0.0617272      0.449539             0      0.633966  

  The vector x:

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

  The product A' * x:

     1:     0.218418
     2:      0.89459
     3:    -0.449539
     4:            0
     5:     0.195544

R8SR_MV_TEST
  R8SR_MV multiplies an R8SR matrix by a vector;

  Matrix order N = 5

  The R8SR matrix, in R8GE form:

  Col:       1             2             3             4             5       
  Row
  ---
   1      0.218418      0.956318             0             0      0.829509  
   2             0      0.561695             0             0      0.415307  
   3      0.257578             0     0.0661187             0             0  
   4      0.043829             0             0      0.109957             0  
   5             0     0.0617272      0.449539             0      0.633966  

  The vector x:

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

  The product A * x:

     1:    -0.611091
     2:    -0.415307
     3:     0.257578
     4:     0.043829
     5:    -0.633966

R8SR_PRINT_TEST
  R8SR_PRINT prints a R8SR matrix.

  Matrix order N = 5

  The R8SR matrix:

  Col:       1             2             3             4             5       
  Row
  ---
   1      0.218418      0.956318             0             0      0.829509  
   2             0      0.561695             0             0      0.415307  
   3      0.257578             0     0.0661187             0             0  
   4      0.043829             0             0      0.109957             0  
   5             0     0.0617272      0.449539             0      0.633966  

R8SR_PRINT_TEST
  R8SR_PRINT prints a R8SR matrix.

  Matrix order N = 5

  Rows 1:N, column 5

  Col:       5       
  Row
  ---
   1      0.829509  
   2      0.415307  
   3             0  
   4             0  
   5      0.633966  

R8SR_RANDOM_TEST
  R8SR_RANDOM randomizes a R8SR matrix

  Matrix order N = 5

  The matrix in R8GE form:

  Col:       1             2             3             4             5       
  Row
  ---
   1      0.218418      0.956318             0             0      0.829509  
   2             0      0.561695             0             0      0.415307  
   3      0.257578             0     0.0661187             0             0  
   4      0.043829             0             0      0.109957             0  
   5             0     0.0617272      0.449539             0      0.633966  

R8SR_TO_R8GE_TEST
  R8SR_TO_R8GE converts an R8SR matrix to R8GE format.

  Matrix order N = 5

  The R8SR matrix:

  Col:       1             2             3             4             5       
  Row
  ---
   1      0.218418      0.956318             0             0      0.829509  
   2             0      0.561695             0             0      0.415307  
   3      0.257578             0     0.0661187             0             0  
   4      0.043829             0             0      0.109957             0  
   5             0     0.0617272      0.449539             0      0.633966  

  The R8GE matrix:

  Col:       1             2             3             4             5       
  Row
  ---
   1      0.218418      0.956318             0             0      0.829509  
   2             0      0.561695             0             0      0.415307  
   3      0.257578             0     0.0661187             0             0  
   4      0.043829             0             0      0.109957             0  
   5             0     0.0617272      0.449539             0      0.633966  

R8SR_ZEROS_TEST
  R8SR_ZEROS sets up a zero R8SR matrix;

  Matrix order N = 5

  The R8SR zero matrix:

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

r8sr_test():
  Normal end of execution.

27-Aug-2022 16:29:08