24 August 2022 09:38:04 PM r8pbl_test(): C version Test r8pbl(). R8PBL_DIF2_TEST R8PBL_DIF2 sets up an R8PBL second difference matrix. Matrix order N = 5 Bandwidth ML = 3 The R8PBL second difference matrix: Col: 0 1 2 3 4 Row --- 0 2 -1 0 0 1 -1 2 -1 0 0 2 0 -1 2 -1 0 3 0 0 -1 2 -1 4 0 0 -1 2 R8PBL_INDICATOR_TEST R8PBL_INDICATOR sets up an R8PBL indicator matrix. Matrix order N = 9 Bandwidth ML = 3 The R8PBL indicator matrix: Col: 0 1 2 3 4 Row --- 0 11 21 31 41 1 21 22 32 42 52 2 31 32 33 43 53 3 41 42 43 44 54 4 52 53 54 55 5 63 64 65 6 74 75 7 85 Col: 5 6 7 8 Row --- 2 63 3 64 74 4 65 75 85 5 66 76 86 96 6 76 77 87 97 7 86 87 88 98 8 96 97 98 99 R8PBL_MV_TEST R8PBL_MV computes A*x, where A is an R8PBL matrix. Matrix order N = 5 Lower bandwidth ML = 2 Matrix A: Col: 0 1 2 3 4 Row --- 0 1.31501 0.218418 0.956318 1 0.218418 1.69061 0.829509 0.561695 2 0.956318 0.829509 3.72095 0.415307 0.0661187 3 0.561695 0.415307 1.3214 0.257578 4 0.0661187 0.257578 0.483706 Vector x: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 Product b=A*x 0 3.664477 1 9.340660 2 15.042566 3 7.544598 4 2.418531 R8PBL_PRINT_TEST R8PBL_PRINT prints an R8PBL matrix. Matrix order N = 9 Bandwidth ML = 3 The R8PBL matrix: Col: 0 1 2 3 4 Row --- 0 11 21 31 41 1 21 22 32 42 52 2 31 32 33 43 53 3 41 42 43 44 54 4 52 53 54 55 5 63 64 65 6 74 75 7 85 Col: 5 6 7 8 Row --- 2 63 3 64 74 4 65 75 85 5 66 76 86 96 6 76 77 87 97 7 86 87 88 98 8 96 97 98 99 R8PBL_PRINT_SOME_TEST R8PBL_PRINT_SOME prints soem of an R8PBL matrix. Matrix order N = 9 Bandwidth ML = 4 Row(3:7), Col(4:8): Col: 4 5 6 7 8 Row --- 3 54 64 74 84 4 55 65 75 85 95 5 65 66 76 86 96 6 75 76 77 87 97 7 85 86 87 88 98 R8PBL_RANDOM_TEST R8PBL_RANDOM randomizes an R8PBL matrix. Matrix order N = 9 Bandwidth ML = 3 The R8PBL random matrix: Col: 0 1 2 3 4 Row --- 0 3.21491 0.218418 0.956318 0.561695 1 0.218418 1.94388 0.829509 0.415307 0.257578 2 0.956318 0.829509 2.62544 0.0661187 0.109957 3 0.561695 0.415307 0.0661187 1.96605 0.043829 4 0.257578 0.109957 0.043829 3.1124 5 0.633966 0.0617272 0.449539 6 0.401306 0.754673 7 0.00183837 Col: 5 6 7 8 Row --- 2 0.633966 3 0.0617272 0.401306 4 0.449539 0.754673 0.00183837 5 3.27926 0.797287 0.897504 0.0945448 6 0.797287 3.1461 0.350752 0.0136169 7 0.897504 0.350752 3.86305 0.859097 8 0.0945448 0.0136169 0.859097 1.23832 R8PBL_TO_R8GE_TEST R8PBL_TO_R8GE converts an R8PBL matrix to R8GE format. Matrix order N = 9 Bandwidth ML = 3 The R8PBL matrix: Col: 0 1 2 3 4 Row --- 0 11 21 31 41 1 21 22 32 42 52 2 31 32 33 43 53 3 41 42 43 44 54 4 52 53 54 55 5 63 64 65 6 74 75 7 85 Col: 5 6 7 8 Row --- 2 63 3 64 74 4 65 75 85 5 66 76 86 96 6 76 77 87 97 7 86 87 88 98 8 96 97 98 99 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 11 21 31 41 0 2 21 22 32 42 52 3 31 32 33 43 53 4 41 42 43 44 54 5 0 52 53 54 55 6 0 0 63 64 65 7 0 0 0 74 75 8 0 0 0 0 85 9 0 0 0 0 0 Col: 6 7 8 9 Row --- 1 0 0 0 0 2 0 0 0 0 3 63 0 0 0 4 64 74 0 0 5 65 75 85 0 6 66 76 86 96 7 76 77 87 97 8 86 87 88 98 9 96 97 98 99 R8PBL_ZEROS_TEST R8PBL_ZEROS zeros an R8PBL matrix. Matrix order N = 5 Bandwidth ML = 3 The R8PBL zero matrix: Col: 0 1 2 3 4 Row --- 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 r8pbl_test(): Normal end of execution. 24 August 2022 09:38:04 PM