Tue May 20 22:27:35 2025 r8pbl_test(): python version: 3.10.12 numpy version: 1.26.4 Test r8pbl(). R8PBL_DIF2_TEST R8PBL_DIF2 sets an R8PBL second difference matrix. Matrix order N = 5 Bandwidth ML = 1 The R8PBL second difference matrix: Col: 0 1 2 3 4 Row 0 : 2 -1 1 : -1 2 -1 2 : -1 2 -1 3 : -1 2 -1 4 : -1 2 R8PBL_INDICATOR_TEST R8PBL_INDICATOR sets up a 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.18302 0.439919 0.276791 1 : 0.439919 2.05887 0.114528 0.655185 2 : 0.276791 0.114528 2.12952 0.608571 0.618362 3 : 0.655185 0.608571 3.00984 0.857902 4 : 0.618362 0.857902 1.98657 Vector x: 0: 1 1: 2 2: 3 3: 4 4: 5 Product b=A*x 0: 2.61644 1: 6.81819 2: 13.8238 3: 20.7577 4: 9.93287 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 some 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 = 5 Bandwidth ML = 2 The R8PBL random matrix: Col: 0 1 2 3 4 Row 0 : 1.56731 0.846558 0.480599 1 : 0.846558 2.37653 0.639083 0.173345 2 : 0.480599 0.639083 2.93114 0.096847 0.864507 3 : 0.173345 0.096847 1.87484 0.964179 4 : 0.864507 0.964179 2.11797 R8PBL_TO_R8GE_TEST R8PBL_TO_R8GE converts an R8PBL matrix to R8GE format. Matrix order N = 5 Bandwidth ML = 2 The R8PBL matrix: Col: 0 1 2 3 4 Row 0 : 1.47053 0.753553 0.536608 1 : 0.753553 4.9329 0.93191 0.949528 2 : 0.536608 0.93191 4.39646 0.494291 0.448852 3 : 0.949528 0.494291 2.97183 0.0529648 4 : 0.448852 0.0529648 0.893882 The R8GE matrix: [[1.47052914 0.75355275 0.53660789 0. 0. ] [0.75355275 4.9328997 0.93191004 0.94952797 0. ] [0.53660789 0.93191004 4.39645956 0.4942913 0.44885198] [0. 0.94952797 0.4942913 2.97183475 0.05296475] [0. 0. 0.44885198 0.05296475 0.89388234]] R8PBL_ZEROS_TEST R8PBL_ZEROS zeros an R8PBL matrix. Matrix order N = 5 Bandwidth ML = 2 The R8PBL zero matrix: Col: 0 1 2 3 4 Row 0 : 0 0 0 1 : 0 0 0 0 2 : 0 0 0 0 0 3 : 0 0 0 0 4 : 0 0 0 r8pbl_test(): Normal end of execution. Tue May 20 22:27:35 2025