06 October 2025 6:33:35.979 PM r8pbl_test(): Fortran90 version: Test r8pbl(). R8PBL_DIF2_TEST R8PBL_DIF2 sets up an R8PBL indicator matrix. Matrix order N = 5 Bandwidth ML = 1 The R8PBL second difference matrix: Col: 1 2 3 4 5 Row --- 1 2.00000 -1.00000 2 -1.00000 2.00000 -1.00000 3 -1.00000 2.00000 -1.00000 4 -1.00000 2.00000 -1.00000 5 -1.00000 2.00000 R8PBL_INDICATOR_TEST R8PBL_INDICATOR sets up an R8PBL indicator matrix. Matrix order N = 5 Bandwidth ML = 2 The R8PBL indicator matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 21.0000 31.0000 2 21.0000 22.0000 32.0000 42.0000 3 31.0000 32.0000 33.0000 43.0000 53.0000 4 42.0000 43.0000 44.0000 54.0000 5 53.0000 54.0000 55.0000 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: 1 2 3 4 5 Row --- 1 1.36241 0.731579 0.105619 2 0.731579 1.97121 0.389412 0.481089 3 0.105619 0.389412 2.11683 0.400860 0.400982 4 0.481089 0.400860 1.78997 0.499283 5 0.400982 0.499283 1.06597 Vector x: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 Product b=A*x 1 3.14242 2 7.76659 3 10.8433 4 11.8211 5 8.52994 R8PBL_PRINT_TEST R8PBL_PRINT prints an R8PBL matrix. Matrix order N = 5 Bandwidth ML = 2 The R8PBL matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 21.0000 31.0000 2 21.0000 22.0000 32.0000 42.0000 3 31.0000 32.0000 33.0000 43.0000 53.0000 4 42.0000 43.0000 44.0000 54.0000 5 53.0000 54.0000 55.0000 R8PBL_PRINT_SOME_TEST R8PBL_PRINT_SOME prints some of an R8PBL matrix. Matrix order N = 9 Bandwidth ML = 4 Row(4:8), Col(5:9): Col: 5 6 7 8 9 Row --- 4 54.0000 64.0000 74.0000 84.0000 5 55.0000 65.0000 75.0000 85.0000 95.0000 6 65.0000 66.0000 76.0000 86.0000 96.0000 7 75.0000 76.0000 77.0000 87.0000 97.0000 8 85.0000 86.0000 87.0000 88.0000 98.0000 R8PBL_RANDOM_TEST R8PBL_RANDOM sets up a random R8PBL matrix. Matrix order N = 5 Bandwidth ML = 2 The random R8PBL matrix: Col: 1 2 3 4 5 Row --- 1 0.978081 0.453030 0.237169 2 0.453030 2.47199 0.689821 0.224430 3 0.237169 0.689821 3.03676 0.388167 0.471775 4 0.224430 0.388167 2.00956 0.463360 5 0.471775 0.463360 1.58851 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: 1 2 3 4 5 Row --- 1 2.79995 0.671481 0.894549 2 0.671481 1.27258 0.362844 0.185922 3 0.894549 0.362844 3.08373 0.110818 0.745744 4 0.185922 0.110818 0.808491 0.155763 5 0.745744 0.155763 1.26608 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 2.79995 0.671481 0.894549 0.00000 0.00000 2 0.671481 1.27258 0.362844 0.185922 0.00000 3 0.894549 0.362844 3.08373 0.110818 0.745744 4 0.00000 0.185922 0.110818 0.808491 0.155763 5 0.00000 0.00000 0.745744 0.155763 1.26608 R8PBL_ZEROS_TEST R8PBL_ZEROS sets up an R8PBL zero matrix. Matrix order N = 5 Bandwidth ML = 2 The R8PBL indicator matrix: Col: 1 2 3 4 5 Row --- 1 0.00000 0.00000 0.00000 2 0.00000 0.00000 0.00000 0.00000 3 0.00000 0.00000 0.00000 0.00000 0.00000 4 0.00000 0.00000 0.00000 0.00000 5 0.00000 0.00000 0.00000 r8pbl_test(): Normal end of execution. 06 October 2025 6:33:35.979 PM