Sat Aug 27 10:57:29 2022 r8sm_test(): Python version: 3.6.9 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.354053 0.00813291 0.40513 0.364666 -0.439075 1 : 0.30222 0.379509 -0.343119 -0.45965 -0.378464 2 : -0.099354 -0.397294 0.487696 -0.17681 0.0588551 3 : 0.27067 -0.316114 -0.350701 -0.210731 -0.315484 4 : 0.149632 -0.264498 -0.0452115 -0.407094 0.0246587 5 : 0.0421225 -0.134432 -0.219335 -0.0758321 0.0133776 6 : 0.106366 0.541913 0.181841 0.553473 0.0523058 Col: 5 6 Row 0 : 0.665061 -0.143448 1 : 0.386484 -0.264483 2 : 0.658623 -0.585023 3 : 0.548863 -0.345215 4 : 0.247486 -0.217992 5 : 0.335677 0.657945 6 : 0.465356 0.642589 A*x and PLU*x 0: 3.83523 3.83523 1: -3.23152 -3.23152 2: 0.0127578 0.0127578 3: -2.95734 -2.95734 4: -2.06111 -2.06111 5: 5.49849 5.49849 6: 11.5014 11.5014 The Sherman Morrison matrix: Col: 0 1 2 3 4 Row 0 : 0.786825 0.375228 -0.0904913 -0.110343 -0.0426225 1 : -0.0457519 0.117927 -0.272812 -0.569809 0.241773 2 : -0.505884 -0.379006 0.0902134 0.0404657 0.090571 3 : 0.915423 0.675858 0.705935 0.626778 -0.00715722 4 : 0.386827 -0.286005 0.205325 -0.249856 0.145109 5 : 0.389476 -0.301176 -0.481372 -0.539264 0.106564 6 : -0.316722 -0.141304 0.0837026 0.169911 -0.173944 Col: 5 6 Row 0 : 0.651282 0.540535 1 : -0.00753499 0.382868 2 : 0.0160905 -0.113202 3 : 0.320778 0.616754 4 : 0.134085 -0.293463 5 : 0.684186 -0.104831 6 : 0.208505 0.0662695 A'*x and (PLU)'*x 0: 4.8933 4.8933 1: -2.04872 -2.04872 2: 1.18257 1.18257 3: -1.91694 -1.91694 4: 0.831334 0.831334 5: 8.20267 8.20267 6: 1.80126 1.80126 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.228919 -0.35722 0.358974 0.390231 1 : 0.982434 0.240895 0.878986 0.0913621 2 : 0.509169 0.292036 0.378554 -0.269201 3 : 0.897196 0.577992 0.418027 0.0671357 4 : 0.489954 0.584306 0.534768 0.707258 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.210721 0.0948478 0.071751 0.323773 0.371511 1 : 0.314321 0.215798 0.581387 0.197304 0.692286 2 : 0.609571 -0.107352 0.305478 0.22388 0.119362 3 : 0.190727 0.517195 0.550696 -0.0247004 -0.208612 4 : 0.828625 0.792905 0.687757 0.485453 0.467904 The right hand side vector B: 0: 7.57411 1: 6.23769 2: 7.79253 3: 3.71849 4: 3.61924 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.199196 -0.155122 -0.10288 0.11587 0.624883 1 : -0.234504 0.25939 0.473314 -0.05192 -0.281691 2 : -0.578583 0.513353 -0.334213 -0.558408 -0.363629 3 : -0.215253 0.609289 -0.0088167 0.432886 0.24979 4 : 0.657182 0.688585 0.192865 0.31855 0.0725715 The right hand side vector B: 0: 2.76982 1: 0.0880845 2: -4.60629 3: 3.95737 4: 4.25001 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. Sat Aug 27 10:57:29 2022