14-Aug-2022 18:49:02 r8bb_test(): MATLAB 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: 1 2 3 4 5 Row --- 1 11 0 0 14 15 2 21 22 0 24 25 3 0 32 33 34 35 4 41 42 43 44 45 5 51 52 53 54 55 The R8BB matrix after additions: Col: 1 2 3 4 5 Row --- 1 111 0 0 214 215 2 21 122 0 224 225 3 0 32 133 234 235 4 41 42 43 44 445 5 51 52 53 454 55 R8BB_DIF2_TEST R8BB_DIF2 sets up an R8BB second difference matrix. Matrix order N = 8 Matrix suborder N1 = 6 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The R8BB 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 6 0 0 0 0 -1 7 0 0 0 0 0 8 0 0 0 0 0 Col: 6 7 8 Row --- 1 0 0 0 2 0 0 0 3 0 0 0 4 0 0 0 5 -1 0 0 6 2 -1 0 7 -1 2 -1 8 0 -1 2 R8BB_FA_TEST R8BB_FA factors an R8BB matrix; Matrix order N = 10 Matrix suborder N1 = 8 Matrix suborder N2 = 2 Lower bandwidth ML = 0 Upper bandwidth MU = 0 Solution: 1: 1 2: 2 3: 3 4: 4 5: 5 6: 6 7: 7 8: 8 9: 9 10: 10 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: 1 2 3 4 5 Row --- 1 11 0 0 14 15 2 21 22 0 24 25 3 0 32 33 34 35 4 41 42 43 44 45 5 51 52 53 54 55 A(2,5) = 25 A(5,3) = 53 A(3,1) = 0 A(2,1) = 21 A(1,4) = 14 A(1,3) = 0 A(3,4) = 34 A(4,1) = 41 A(5,2) = 52 A(1,1) = 11 R8BB_INDICATOR_TEST R8BB_INDICATOR sets up an indicator matrix. Matrix order N = 8 Matrix suborder N1 = 6 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The indicator matrix: Col: 1 2 3 4 5 Row --- 1 11 12 0 0 0 2 21 22 23 0 0 3 0 32 33 34 0 4 0 0 43 44 45 5 0 0 0 54 55 6 0 0 0 0 65 7 71 72 73 74 75 8 81 82 83 84 85 Col: 6 7 8 Row --- 1 0 17 18 2 0 27 28 3 0 37 38 4 0 47 48 5 56 57 58 6 66 67 68 7 76 77 78 8 86 87 88 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: 1 2 3 4 5 Row --- 1 11 12 0 0 0 2 21 22 23 0 0 3 0 32 33 34 0 4 0 0 43 44 45 5 0 0 0 54 55 6 0 0 0 0 65 7 71 72 73 74 75 8 81 82 83 84 85 Col: 6 7 8 Row --- 1 0 17 18 2 0 27 28 3 0 37 38 4 0 47 48 5 56 57 58 6 66 67 68 7 76 77 78 8 86 87 88 The vector x: 1: 1 2: 2 3: 3 4: 4 5: 5 6: 6 7: 7 8: 8 The product b=A'*x: 1: 1198 2: 1312 3: 1492 4: 1738 5: 2050 6: 1896 7: 2292 8: 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: 1 2 3 4 5 Row --- 1 11 12 0 0 0 2 21 22 23 0 0 3 0 32 33 34 0 4 0 0 43 44 45 5 0 0 0 54 55 6 0 0 0 0 65 7 71 72 73 74 75 8 81 82 83 84 85 Col: 6 7 8 Row --- 1 0 17 18 2 0 27 28 3 0 37 38 4 0 47 48 5 56 57 58 6 66 67 68 7 76 77 78 8 86 87 88 The vector x: 1: 1 2: 2 3: 3 4: 4 5: 5 6: 6 7: 7 8: 8 The product b=A*x: 1: 298 2: 547 3: 862 4: 1243 5: 1690 6: 1734 7: 2724 8: 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: 1 2 3 4 5 Row --- 1 0.218418 0.829509 0 0 0 2 0.956318 0.561695 0.0661187 0 0 3 0 0.415307 0.257578 0.043829 0 4 0 0 0.109957 0.633966 0.449539 5 0 0 0 0.0617272 0.401306 6 0 0 0 0 0.754673 7 0 0 0 0 0 8 0 0 0 0 0 9 0.396988 0.574366 0.617205 0.21293 0.117707 10 0.185314 0.367027 0.361529 0.714471 0.299329 Col: 6 7 8 9 10 Row --- 1 0 0 0 0.123104 0.692066 2 0 0 0 0.00751236 0.561662 3 0 0 0 0.260303 0.861216 4 0 0 0 0.912484 0.453794 5 0.797287 0 0 0.113664 0.911977 6 0.00183837 0.350752 0 0.351629 0.597917 7 0.897504 0.0945448 0.859097 0.822887 0.188955 8 0 0.0136169 0.840847 0.267132 0.761492 9 0.825003 0.0618618 0.0882833 0.745303 0.899373 10 0.82466 0.710781 0.777994 0.308675 0.763537 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: 7 8 Row --- 7 77 78 8 87 88 R8BB_RANDOM_TEST R8BB_RANDOM returns a random R8BB matrix; Matrix order N = 10 Matrix suborder N1 = 8 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The border-banded matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.829509 0 0 0 2 0.956318 0.561695 0.0661187 0 0 3 0 0.415307 0.257578 0.043829 0 4 0 0 0.109957 0.633966 0.449539 5 0 0 0 0.0617272 0.401306 6 0 0 0 0 0.754673 7 0 0 0 0 0 8 0 0 0 0 0 9 0.396988 0.574366 0.617205 0.21293 0.117707 10 0.185314 0.367027 0.361529 0.714471 0.299329 Col: 6 7 8 9 10 Row --- 1 0 0 0 0.123104 0.692066 2 0 0 0 0.00751236 0.561662 3 0 0 0 0.260303 0.861216 4 0 0 0 0.912484 0.453794 5 0.797287 0 0 0.113664 0.911977 6 0.00183837 0.350752 0 0.351629 0.597917 7 0.897504 0.0945448 0.859097 0.822887 0.188955 8 0 0.0136169 0.840847 0.267132 0.761492 9 0.825003 0.0618618 0.0882833 0.745303 0.899373 10 0.82466 0.710781 0.777994 0.308675 0.763537 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: 1 2 3 4 5 Row --- 1 11 12 0 0 15 2 21 22 23 0 25 3 31 32 33 34 35 4 0 42 43 44 45 5 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: 1: 1 2: 2 3: 3 4: 4 5: 5 6: 6 7: 7 8: 8 9: 9 10: 10 R8BB_TO_R8GE_TEST R8BB_TO_R8GE converts an R8BB matrix to R8GE format. Matrix order N = 8 Matrix suborder N1 = 6 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The R8BB matrix: Col: 1 2 3 4 5 Row --- 1 11 12 0 0 0 2 21 22 23 0 0 3 0 32 33 34 0 4 0 0 43 44 45 5 0 0 0 54 55 6 0 0 0 0 65 7 71 72 73 74 75 8 81 82 83 84 85 Col: 6 7 8 Row --- 1 0 17 18 2 0 27 28 3 0 37 38 4 0 47 48 5 56 57 58 6 66 67 68 7 76 77 78 8 86 87 88 ij = 25 b = 11 12 0 0 0 0 17 0 21 22 23 0 0 0 0 0 0 32 33 34 0 0 0 0 0 0 43 44 45 0 0 0 0 0 0 54 55 56 0 0 0 0 0 0 65 66 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ij = 31 b = 11 12 0 0 0 0 17 18 21 22 23 0 0 0 0 0 0 32 33 34 0 0 0 0 0 0 43 44 45 0 0 0 0 0 0 54 55 56 0 0 0 0 0 0 65 66 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ij = 26 b = 11 12 0 0 0 0 17 18 21 22 23 0 0 0 27 0 0 32 33 34 0 0 0 0 0 0 43 44 45 0 0 0 0 0 0 54 55 56 0 0 0 0 0 0 65 66 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ij = 32 b = 11 12 0 0 0 0 17 18 21 22 23 0 0 0 27 28 0 32 33 34 0 0 0 0 0 0 43 44 45 0 0 0 0 0 0 54 55 56 0 0 0 0 0 0 65 66 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ij = 27 b = 11 12 0 0 0 0 17 18 21 22 23 0 0 0 27 28 0 32 33 34 0 0 37 0 0 0 43 44 45 0 0 0 0 0 0 54 55 56 0 0 0 0 0 0 65 66 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ij = 33 b = 11 12 0 0 0 0 17 18 21 22 23 0 0 0 27 28 0 32 33 34 0 0 37 38 0 0 43 44 45 0 0 0 0 0 0 54 55 56 0 0 0 0 0 0 65 66 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ij = 28 b = 11 12 0 0 0 0 17 18 21 22 23 0 0 0 27 28 0 32 33 34 0 0 37 38 0 0 43 44 45 0 47 0 0 0 0 54 55 56 0 0 0 0 0 0 65 66 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ij = 34 b = 11 12 0 0 0 0 17 18 21 22 23 0 0 0 27 28 0 32 33 34 0 0 37 38 0 0 43 44 45 0 47 48 0 0 0 54 55 56 0 0 0 0 0 0 65 66 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ij = 29 b = 11 12 0 0 0 0 17 18 21 22 23 0 0 0 27 28 0 32 33 34 0 0 37 38 0 0 43 44 45 0 47 48 0 0 0 54 55 56 57 0 0 0 0 0 65 66 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ij = 35 b = 11 12 0 0 0 0 17 18 21 22 23 0 0 0 27 28 0 32 33 34 0 0 37 38 0 0 43 44 45 0 47 48 0 0 0 54 55 56 57 58 0 0 0 0 65 66 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ij = 30 b = 11 12 0 0 0 0 17 18 21 22 23 0 0 0 27 28 0 32 33 34 0 0 37 38 0 0 43 44 45 0 47 48 0 0 0 54 55 56 57 58 0 0 0 0 65 66 67 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ij = 36 b = 11 12 0 0 0 0 17 18 21 22 23 0 0 0 27 28 0 32 33 34 0 0 37 38 0 0 43 44 45 0 47 48 0 0 0 54 55 56 57 58 0 0 0 0 65 66 67 68 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 11 12 0 0 0 2 21 22 23 0 0 3 0 32 33 34 0 4 0 0 43 44 45 5 0 0 0 54 55 6 0 0 0 0 65 7 71 72 73 74 75 8 81 82 83 84 85 Col: 6 7 8 Row --- 1 0 17 18 2 0 27 28 3 0 37 38 4 0 47 48 5 56 57 58 6 66 67 68 7 76 77 78 8 86 87 88 R8BB_ZEROS_TEST R8BB_ZEROS zeros an R8BB matrix. Matrix order N = 8 Matrix suborder N1 = 6 Matrix suborder N2 = 2 Lower bandwidth ML = 1 Upper bandwidth MU = 1 The zero R8BB 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 6 0 0 0 0 0 7 0 0 0 0 0 8 0 0 0 0 0 Col: 6 7 8 Row --- 1 0 0 0 2 0 0 0 3 0 0 0 4 0 0 0 5 0 0 0 6 0 0 0 7 0 0 0 8 0 0 0 r8bb_test(): Normal end of execution. 14-Aug-2022 18:49:02