10 April 2023 10:12:10.327 AM r8sm_test(): FORTRAN90 version: Test rsm(). R8SM_INDICATOR_TEST R8SM_INDICATOR sets up an R8SM indicator matrix; M = 5 N = 4 The R8SM indicator matrix: Col: 1 2 3 4 Row --- 1 11.0000 12.0000 13.0000 14.0000 2 21.0000 22.0000 23.0000 24.0000 3 31.0000 32.0000 33.0000 34.0000 4 41.0000 42.0000 43.0000 44.0000 5 51.0000 52.0000 53.0000 54.0000 R8SM_ML_TEST R8SM_ML computes A*x or A'*X where A is a Sherman Morrison matrix. Matrix rows M = 7 Matrix columns N = 7 The Sherman Morrison matrix: Col: 1 2 3 4 5 Row --- 1 -0.523259 -0.519701 0.169118 0.505236 0.487381E-01 2 0.900681 -0.340474E-02 -0.452836E-01 0.979281E-01 0.209060 3 0.190252 0.912587E-01 0.356081 -0.281754 0.248217E-01 4 0.482296 -0.568034E-02 0.283504 0.224374 0.478787 5 -0.284400 -0.144488 -0.498077 0.595864 0.130874 6 -0.604187 -0.167760 -0.554103 -0.189652 -0.245859 7 -0.200361E-01 0.518989 0.623970 0.226007 0.622209 Col: 6 7 Row --- 1 0.134385 0.602509 2 0.154183 0.360426 3 0.361971 0.200264 4 0.347365 0.712898 5 -0.251987 0.103843 6 0.155440 0.286048 7 0.193524 0.819502 A*x and PLU*x 1 6.23320 6.23320 2 5.64311 5.64311 3 4.01178 4.01178 4 11.6874 11.6874 5 0.185197 0.185197 6 -1.65495 -1.65495 7 13.8026 13.8026 The Sherman Morrison matrix: Col: 1 2 3 4 5 Row --- 1 0.372211 0.446500E-01 0.806618 0.931129 0.533214 2 -0.234486 -0.116704 0.580868 0.879274 0.123585 3 0.256648 0.260123 0.212966 0.507635 0.397534E-01 4 -0.460832 -0.492993 0.751566 -0.487418E-01 0.681334E-01 5 0.554236 -0.477249E-01 0.716446 0.448273 0.646527 6 -0.803782 -0.483181 0.429272 0.202573 -0.129083 7 -0.627499 0.738779E-01 0.840658 0.833144 -0.217538 Col: 6 7 Row --- 1 0.577074 0.266591 2 0.242808 0.325956 3 0.473955 0.296694 4 0.397568 -0.274456 5 -0.348438E-01 0.651523 6 0.222917 0.143123 7 0.392929 0.434966 A'*x and (PLU)'*x 1 -7.61416 -7.61416 2 -4.00093 -4.00093 3 17.6560 17.6560 4 13.3064 13.3064 5 2.10755 2.10755 6 7.98862 7.98862 7 7.87188 7.87188 R8SM_MTV_TEST R8SM_MTV computes A'*x=b, where A is an R8SM matrix; M = 5 N = 4 The R8SM matrix: Col: 1 2 3 4 Row --- 1 11.0000 12.0000 13.0000 14.0000 2 21.0000 22.0000 23.0000 24.0000 3 31.0000 32.0000 33.0000 34.0000 4 41.0000 42.0000 43.0000 44.0000 5 51.0000 52.0000 53.0000 54.0000 The vector x: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 The product b=A'*x: 1 565.000 2 580.000 3 595.000 4 610.000 R8SM_MV_TEST R8SM_MV computes A*x=b, where A is an R8SM matrix; M = 5 N = 4 The R8SM matrix: Col: 1 2 3 4 Row --- 1 11.0000 12.0000 13.0000 14.0000 2 21.0000 22.0000 23.0000 24.0000 3 31.0000 32.0000 33.0000 34.0000 4 41.0000 42.0000 43.0000 44.0000 5 51.0000 52.0000 53.0000 54.0000 The vector x: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 The product b=A*x: 1 130.000 2 230.000 3 330.000 4 430.000 5 530.000 R8SM_PRINT_TEST R8SM_PRINT prints an R8SM matrix; M = 5 N = 4 The R8SM matrix: Col: 1 2 3 4 Row --- 1 11.0000 12.0000 13.0000 14.0000 2 21.0000 22.0000 23.0000 24.0000 3 31.0000 32.0000 33.0000 34.0000 4 41.0000 42.0000 43.0000 44.0000 5 51.0000 52.0000 53.0000 54.0000 R8SM_PRINT_SOME_TEST R8SM_PRINT_SOME prints some of an R8SM matrix; M = 9 N = 9 Rows 2-5, Cols 3-7: Col: 3 4 5 6 7 Row --- 2 23.0000 24.0000 25.0000 26.0000 27.0000 3 33.0000 34.0000 35.0000 36.0000 37.0000 4 43.0000 44.0000 45.0000 46.0000 47.0000 5 53.0000 54.0000 55.0000 56.0000 57.0000 R8SM_RANDOM_TEST R8SM_RANDOM sets up a random R8SM matrix; M = 5 N = 4 The R8SM matrix: Col: 1 2 3 4 Row --- 1 -0.565494 -0.315297E-01 -0.240356 -0.705102 2 0.189058 0.162004 0.153873 0.205581 3 0.717179 0.959643E-01 0.358019 0.249452 4 0.554841 0.429751E-01 0.752032 0.883629E-01 5 0.177785 0.604379 0.705757 -0.200583 R8SM_SL_TEST R8SM_SL solves B*x=b, where B is an R8SM matrix. Matrix rows M = 5 Matrix columns N = 5 The Sherman-Morrison matrix A: Col: 1 2 3 4 5 Row --- 1 -0.413081 -0.446389 -0.724118 -0.412241 0.269327E-01 2 0.877654 0.193737 0.351650 0.845924 0.737188 3 0.586159 -0.875398E-01 0.984781E-01 0.191185 -0.197573 4 -0.779730E-02 -0.418356 0.459899E-01 -0.278876 -0.742942 5 0.230434 0.483928 0.567228 -0.107606 0.166845E-01 The right hand side vector B: 1 -4.99251 2 9.38972 3 0.483389 4 -5.53675 5 2.55297 Solution to A * X = B: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 R8SM_SLT_TEST R8SM_SLT solves B'*x=b, where B is an R8SM matrix. Matrix rows M = 5 Matrix columns N = 5 The Sherman-Morrison matrix A: Col: 1 2 3 4 5 Row --- 1 -0.413081 -0.446389 -0.724118 -0.412241 0.269327E-01 2 0.877654 0.193737 0.351650 0.845924 0.737188 3 0.586159 -0.875398E-01 0.984781E-01 0.191185 -0.197573 4 -0.779730E-02 -0.418356 0.459899E-01 -0.278876 -0.742942 5 0.230434 0.483928 0.567228 -0.107606 0.166845E-01 The right hand side vector B: 1 4.22168 2 0.424679 3 3.29472 4 0.199627 5 -1.97976 Solution to A' * X = B: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 R8SM_TO_R8GE_TEST R8SM_TO_R8GE converts an R8SM matrix to R8GE format. M = 5 N = 4 The R8SM matrix: Col: 1 2 3 4 Row --- 1 11.0000 12.0000 13.0000 14.0000 2 21.0000 22.0000 23.0000 24.0000 3 31.0000 32.0000 33.0000 34.0000 4 41.0000 42.0000 43.0000 44.0000 5 51.0000 52.0000 53.0000 54.0000 The R8GE matrix: Col: 1 2 3 4 Row --- 1 11.0000 12.0000 13.0000 14.0000 2 21.0000 22.0000 23.0000 24.0000 3 31.0000 32.0000 33.0000 34.0000 4 41.0000 42.0000 43.0000 44.0000 5 51.0000 52.0000 53.0000 54.0000 R8SM_ZEROS_TEST R8SM_ZEROS sets up a zero R8SM matrix; M = 5 N = 4 The R8SM zero matrix: Col: 1 2 3 4 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 r8sm_test(): Normal end of execution. 10 April 2023 10:12:10.328 AM