Sun Aug 14 19:26:59 2022 r8bb_test(): Python version: Test r8bb(). R8BB_ADD_TEST R8BB_ADD adds a value to elements of an R8BB matrix. Matrix order N = 5 Matrix suborder N1 = 3 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 0 Matrix before additions: Col: 0 1 2 3 4 Row 0 : 11 0 0 14 15 1 : 21 22 0 24 25 2 : 0 32 33 34 35 3 : 41 42 43 44 45 4 : 51 52 53 54 55 The R8BB matrix after additions: Col: 0 1 2 3 4 Row 0 : 111 0 0 214 215 1 : 21 122 0 224 225 2 : 0 32 133 234 235 3 : 41 42 43 44 445 4 : 51 52 53 454 55 R8BB_DIF2_TEST R8BB_DIF2 sets up an R8BB second difference matrix. Matrix order N = 5 Matrix suborder N1 = 3 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The R8BB second difference matrix: Col: 0 1 2 3 4 Row 0 : 2 -1 0 0 0 1 : -1 2 -1 0 0 2 : 0 -1 2 -1 0 3 : 0 0 -1 2 -1 4 : 0 0 0 -1 2 R8BB_FA_TEST R8BB_FA factors an R8BB matrix Matrix order N = 5 Matrix suborder N1 = 3 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 Solution: 0: 12.9481 1: 10.9006 2: 7.81533 3: 1.14833 4: -17.6524 R8BB_GET_TEST R8BB_GET gets a value of an element of an R8BB matrix. Matrix order N = 5 Matrix suborder N1 = 3 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 0 The matrix to be queried: Col: 0 1 2 3 4 Row 0 : 11 0 0 14 15 1 : 21 22 0 24 25 2 : 0 32 33 34 35 3 : 41 42 43 44 45 4 : 51 52 53 54 55 A(3,0) = 41 A(1,3) = 24 A(0,4) = 15 A(0,0) = 11 A(3,2) = 43 A(3,3) = 44 A(4,1) = 52 A(1,0) = 21 A(3,2) = 43 A(3,2) = 43 R8BB_INDICATOR_TEST R8BB_INDICATOR sets up an indicator matrix. Matrix order N = 5 Matrix suborder N1 = 3 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The indicator matrix: Col: 0 1 2 3 4 Row 0 : 11 12 0 14 15 1 : 21 22 23 24 25 2 : 0 32 33 34 35 3 : 41 42 43 44 45 4 : 51 52 53 54 55 R8BB_MTV_TEST R8BB_MTV computes b=A*x, where A is an R8BB matrix. Matrix order N = 8 Matrix suborder N1 = 6 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The R8BB matrix A: Col: 0 1 2 3 4 Row 0 : 11 12 0 0 0 1 : 21 22 23 0 0 2 : 0 32 33 34 0 3 : 0 0 43 44 45 4 : 0 0 0 54 55 5 : 0 0 0 0 65 6 : 71 72 73 74 75 7 : 81 82 83 84 85 Col: 5 6 7 Row 0 : 0 17 18 1 : 0 27 28 2 : 0 37 38 3 : 0 47 48 4 : 56 57 58 5 : 66 67 68 6 : 76 77 78 7 : 86 87 88 The vector x: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 The product b=A*x: 0: 1198 1: 1312 2: 1492 3: 1738 4: 2050 5: 1896 6: 2694 7: 2328 R8BB_MV_TEST R8BB_MV computes b=A*x, where A is an R8BB matrix. Matrix order N = 8 Matrix suborder N1 = 6 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The R8BB matrix A: Col: 0 1 2 3 4 Row 0 : 11 12 0 0 0 1 : 21 22 23 0 0 2 : 0 32 33 34 0 3 : 0 0 43 44 45 4 : 0 0 0 54 55 5 : 0 0 0 0 65 6 : 71 72 73 74 75 7 : 81 82 83 84 85 Col: 5 6 7 Row 0 : 0 17 18 1 : 0 27 28 2 : 0 37 38 3 : 0 47 48 4 : 56 57 58 5 : 66 67 68 6 : 76 77 78 7 : 86 87 88 The vector x: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 The product b=A*x: 0: 298 1: 547 2: 862 3: 1243 4: 1690 5: 2203 6: 2724 7: 3084 R8BB_PRINT_TEST R8BB_PRINT prints an R8BB matrix Matrix order N = 10 Matrix suborder N1 = 8 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The R8BB matrix: Col: 0 1 2 3 4 Row 0 : 0.505443 0.423405 0 0 0 1 : 0.295817 0.684293 0.236621 0 0 2 : 0 0.526706 0.24573 0.569739 0 3 : 0 0 0.389176 0.776168 0.472914 4 : 0 0 0 0.331929 0.0521081 5 : 0 0 0 0 0.732028 6 : 0 0 0 0 0 7 : 0 0 0 0 0 8 : 0.39133 0.286743 0.156305 0.234964 0.278939 9 : 0.176637 0.702952 0.605129 0.573207 0.515506 Col: 5 6 7 8 9 Row 0 : 0 0 0 0.473663 0.864556 1 : 0 0 0 0.359795 0.557291 2 : 0 0 0 0.419137 0.535366 3 : 0 0 0 0.36323 0.603084 4 : 0.100883 0 0 0.602599 0.526312 5 : 0.170148 0.426399 0 0.25303 0.0763787 6 : 0.725945 0.981512 0.953648 0.517566 0.24299 7 : 0 0.989189 0.87947 0.806935 0.88236 8 : 0.909798 0.640285 0.969821 0.461925 0.0336699 9 : 0.754398 0.510318 0.423898 0.0535408 0.0019223 R8BB_PRINT_SOME_TEST R8BB_PRINT_SOME prints some of an R8BB matrix Matrix order N = 8 Matrix suborder N1 = 6 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The Lower Right Block: Col: 6 7 Row 6 : 77 78 7 : 87 88 R8BB_RANDOM_TEST R8BB_RANDOM returns a random R8BB matrix Matrix order N = 5 Matrix suborder N1 = 3 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The border-banded matrix: Col: 0 1 2 3 4 Row 0 : 0.0608006 0.249707 0 0.389206 0.883799 1 : 0.212621 0.611199 0.0535265 0.98992 0.331014 2 : 0 0.962543 0.367561 0.217534 0.495193 3 : 0.923393 0.268867 0.237095 0.586899 0.24873 4 : 0.734528 0.646614 0.0825888 0.100583 0.923361 R8BB_SET_TEST R8BB_SET sets elements of an R8BB matrix. Matrix order N = 5 Matrix suborder N1 = 4 Matrix suborder N2 = 1 Lower bandwidth ML = 2 Upper bandwidth MU = 1 The R8BB matrix: Col: 0 1 2 3 4 Row 0 : 11 12 0 0 15 1 : 21 22 23 0 25 2 : 31 32 33 34 35 3 : 0 42 43 44 45 4 : 51 52 53 54 55 R8BB_SL_TEST R8BB_SL solves a linear system factored by R8BB_FA. Matrix order N = 10 Matrix suborder N1 = 8 Matrix suborder N2 = 2 Lower bandwidth ML = 0 Upper bandwidth MU = 0 Solution: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 8: 9 9: 10 R8BB_TO_R8GE_TEST R8BB_TO_R8GE converts an R8BB matrix to R8GE format. Matrix order N = 5 Matrix suborder N1 = 3 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The R8BB matrix: Col: 0 1 2 3 4 Row 0 : 11 12 0 14 15 1 : 21 22 23 24 25 2 : 0 32 33 34 35 3 : 41 42 43 44 45 4 : 51 52 53 54 55 The R8GE matrix: [[11. 12. 0. 14. 15.] [21. 22. 23. 24. 25.] [ 0. 32. 33. 34. 35.] [41. 42. 43. 44. 45.] [51. 52. 53. 54. 55.]] R8BB_ZEROS_TEST R8BB_ZEROS zeros an R8BB matrix. Matrix order N = 5 Matrix suborder N1 = 3 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The zero R8BB matrix: Col: 0 1 2 3 4 Row 0 : 0 0 0 0 0 1 : 0 0 0 0 0 2 : 0 0 0 0 0 3 : 0 0 0 0 0 4 : 0 0 0 0 0 r8bb_test(): Normal end of execution. Sun Aug 14 19:26:59 2022