Thu Aug 25 07:16:13 2022 r8pbl_test(): Python version: 3.6.9 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.41771 0.653641 0.41536 1 : 0.653641 3.30208 0.529159 0.930535 2 : 0.41536 0.529159 2.78188 0.792502 0.0176733 3 : 0.930535 0.792502 4.58675 0.882762 4 : 0.0176733 0.882762 1.09623 Vector x: 0: 1 1: 2 2: 3 3: 4 4: 5 Product b=A*x 0: 3.55571 1: 11.9278 2: 15.0324 3: 25.4621 4: 5.48115 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.28886 0.38107 0.882014 1 : 0.38107 2.11055 0.148462 0.922832 2 : 0.882014 0.148462 4.38434 0.922334 0.96016 3 : 0.922832 0.922334 3.71059 0.218879 4 : 0.96016 0.218879 2.15854 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.37616 0.105445 0.875834 1 : 0.105445 2.38071 0.74911 0.57691 2 : 0.875834 0.74911 2.56773 0.107149 0.073588 3 : 0.57691 0.107149 1.63075 0.544985 4 : 0.073588 0.544985 0.923162 The R8GE matrix: [[1.37616025 0.105445 0.87583359 0. 0. ] [0.105445 2.38071339 0.74911008 0.57690951 0. ] [0.87583359 0.74911008 2.56772904 0.10714906 0.07358801] [0. 0.57690951 0.10714906 1.63074568 0.54498525] [0. 0. 0.07358801 0.54498525 0.92316189]] 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. Thu Aug 25 07:16:13 2022