Tue May 20 22:27:38 2025 r8sm_test(): python version: 3.10.12 numpy version: 1.26.4 Test r8sm(). R8SM_INDICATOR_TEST R8SM_INDICATOR sets up an R8SM indicator matrix M = 5 N = 4 The R8SM indicator matrix: Col: 0 1 2 3 Row 0 : 11 12 13 14 1 : 21 22 23 24 2 : 31 32 33 34 3 : 41 42 43 44 4 : 51 52 53 54 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: 0 1 2 3 4 Row 0 : -0.477252 0.328188 -0.508719 -0.552461 -0.0488291 1 : 0.608136 0.0344015 -0.217362 0.542245 0.14612 2 : 0.332458 0.72197 0.311996 0.0307884 0.414131 3 : 0.0513574 0.126306 0.432506 0.661439 0.912905 4 : -0.0771845 0.304584 -0.153432 0.840788 0.880882 5 : -0.329604 0.679348 -0.13912 -0.620351 0.0294452 6 : -0.0849737 0.52607 -0.364791 0.169264 -0.185535 Col: 5 6 Row 0 : 0.242799 -0.250765 1 : 0.444736 -0.198758 2 : 0.314718 0.340365 3 : -0.0264301 0.387526 4 : 0.64686 0.517384 5 : -0.0181551 -0.183843 6 : 0.29801 -0.543375 A*x and PLU*x 0: -4.09958 -4.09958 1: 4.20154 4.20154 2: 9.17705 9.17705 3: 11.3659 11.3659 4: 15.3421 15.3421 5: -3.11828 -3.11828 6: -2.39339 -2.39339 The Sherman Morrison matrix: Col: 0 1 2 3 4 Row 0 : -0.0250682 0.49273 0.647701 0.449501 0.0695067 1 : 0.162477 -0.168487 -0.0370388 -0.340765 0.361968 2 : 0.511761 0.538428 0.48114 0.183381 0.480138 3 : 0.434243 0.115365 0.0318934 0.412913 0.255405 4 : 0.359541 -0.0122373 0.0519866 0.529672 0.224986 5 : -0.0110074 0.422812 0.75213 0.4977 0.884939 6 : 0.676613 0.196783 0.562453 0.14529 0.72589 Col: 5 6 Row 0 : 0.406707 0.25502 1 : -0.0473622 0.531337 2 : 0.191185 0.015902 3 : -0.174167 0.456008 4 : 0.4134 0.065638 5 : 0.505431 0.293554 6 : -0.191574 0.00622485 A'*x and (PLU)'*x 0: 10.0401 10.0401 1: 6.08567 6.08567 2: 10.8545 10.8545 3: 8.62136 8.62136 4: 14.7713 14.7713 5: 3.94743 3.94743 6: 5.32252 5.32252 R8SM_MTV_TEST R8SM_MTV computes A'*x=b, where A is an R8SM matrix M = 5 N = 4 The R8SM matrix: Col: 0 1 2 3 Row 0 : 11 12 13 14 1 : 21 22 23 24 2 : 31 32 33 34 3 : 41 42 43 44 4 : 51 52 53 54 The vector x: 0: 1 1: 2 2: 3 3: 4 4: 5 The product b=A*x: 0: 565 1: 580 2: 595 3: 610 R8SM_MV_TEST R8SM_MV computes A*x=b, where A is an R8SM matrix M = 5 N = 4 The R8SM matrix: Col: 0 1 2 3 Row 0 : 11 12 13 14 1 : 21 22 23 24 2 : 31 32 33 34 3 : 41 42 43 44 4 : 51 52 53 54 The vector x: 0: 1 1: 2 2: 3 3: 4 The product b=A*x: 0: 130 1: 230 2: 330 3: 430 4: 530 R8SM_PRINT_TEST R8SM_PRINT prints an R8SM matrix M = 5 N = 4 The R8SM matrix: Col: 0 1 2 3 Row 0 : 11 12 13 14 1 : 21 22 23 24 2 : 31 32 33 34 3 : 41 42 43 44 4 : 51 52 53 54 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 : 34 35 36 37 38 3 : 44 45 46 47 48 4 : 54 55 56 57 58 5 : 64 65 66 67 68 R8SM_RANDOM_TEST R8SM_RANDOM sets up an R8SM random matrix M = 5 N = 4 The R8SM random matrix: Col: 0 1 2 3 Row 0 : 0.0310923 0.806756 0.364206 0.82415 1 : 0.419155 0.38898 0.248968 0.423122 2 : 0.724025 0.347063 0.249229 0.893791 3 : 0.124904 0.0231781 -0.0640952 0.854416 4 : 0.544734 0.859279 -0.270115 0.583247 R8SM_SL_TEST R8SM_SL implements the Sherman-Morrison method for solving a perturbed linear system. Matrix rows M = 5 Matrix columns N = 5 The Sherman-Morrison matrix A: Col: 0 1 2 3 4 Row 0 : 0.387814 0.237066 0.735596 0.806735 0.058022 1 : 0.586674 0.457573 0.192665 0.0338868 0.767461 2 : 0.349546 -0.0584846 0.756639 0.0550642 -0.168821 3 : 0.290842 -0.311814 0.159215 0.182219 0.0814848 4 : 0.854179 0.485078 0.700005 -0.198945 0.771415 The right hand side vector B: 0: 8.04407 1: 2.1549 2: 7.52772 3: 0.773852 4: 5.2695 Solution to A * X = B: 0: 1 1: 2 2: 3 3: 4 4: 5 The Sherman-Morrison matrix A: Col: 0 1 2 3 4 Row 0 : -0.3835 0.0324105 -0.0691325 0.283456 0.0624166 1 : 0.43991 0.2857 0.724146 0.84497 0.334315 2 : -0.0838368 0.203386 -0.160915 0.00753862 0.123489 3 : -0.0794646 0.0761025 0.174748 0.391299 0.365701 4 : 0.309504 0.0507996 -0.137372 0.277075 0.00637695 The right hand side vector B: 0: 0.91983 1: 8.2352 2: 0.487789 3: 3.99069 4: 1.13917 Solution to A * X = B: 0: 1 1: 2 2: 3 3: 4 4: 5 R8SM_TO_R8GE_TEST R8SM_TO_R8GE converts an R8SM matrix to R8GE format. M = 5 N = 4 The R8SM matrix: Col: 0 1 2 3 Row 0 : 11 12 13 14 1 : 21 22 23 24 2 : 31 32 33 34 3 : 41 42 43 44 4 : 51 52 53 54 The R8GE matrix: [[11. 12. 13. 14.] [21. 22. 23. 24.] [31. 32. 33. 34.] [41. 42. 43. 44.] [51. 52. 53. 54.]] R8SM_ZEROS_TEST R8SM_ZEROS sets up an R8SM zero matrix M = 5 N = 4 The R8SM zeros matrix: Col: 0 1 2 3 Row 0 : 0 0 0 0 1 : 0 0 0 0 2 : 0 0 0 0 3 : 0 0 0 0 4 : 0 0 0 0 r8sm_test(): Normal end of execution. Tue May 20 22:27:38 2025