Sat Aug 20 09:41:38 2022 r8but_test(): Python version: 3.6.9 Test r8but(). R8BUT_DET_TEST R8BUT_DET computes the determinant of an R8BUT matrix. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 0 1 2 3 4 Row --- 0 0.756827 0.528667 0.820624 0.402478 1 0.706235 0.238708 0.566238 0.232102 2 0.184688 0.973678 0.574877 3 0.565137 0.231686 4 0.760477 Determinant = 0.042425373367055125 R8BUT_INDICATOR_TEST R8BUT_INDICATOR sets up an R8BUT indicator matrix. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 0 1 2 3 4 Row --- 0 11 12 13 14 1 22 23 24 25 2 33 34 35 3 44 45 4 55 R8BUT_MTV_TEST R8BUT_MTV computes b=A'*x, where A is an R8BUT matrix. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 0 1 2 3 4 Row --- 0 0.712814 0.144349 0.486341 0.370178 1 0.917152 0.0896662 0.0644392 0.547506 2 0.846924 0.447638 0.414278 3 0.378999 0.740387 4 0.328929 x: 0: 1 1: 2 2: 3 3: 4 4: 5 b=A'*x: 0: 0.712814 1: 1.97865 2: 3.20645 3: 3.35797 4: 6.94404 R8BUT_MV_TEST R8BUT_MV computes b=A*x, where A is an R8BUT matrix. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 0 1 2 3 4 Row --- 0 0.168357 0.123795 0.696516 0.192694 1 0.590975 0.844963 0.0300269 0.540938 2 0.856371 0.00196726 0.223069 3 0.898586 0.821306 4 0.854136 x: 0: 1 1: 2 2: 3 3: 4 4: 5 b=A*x: 0: 3.27627 1: 6.54164 2: 3.69232 3: 7.70088 4: 4.27068 R8BUT_PRINT_TEST R8BUT_PRINT prints an R8BUT matrix. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 0 1 2 3 4 Row --- 0 11 12 13 14 1 22 23 24 25 2 33 34 35 3 44 45 4 55 R8BUT_PRINT_SOME_TEST R8BUT_PRINT_SOME prints some of an R8BUT matrix. Matrix order N = 10 Upper bandwidth MU = 3 Rows 1:4, Cols 2:4: Col: 2 3 4 Row --- 2 303 304 305 3 404 405 4 505 R8BUT_RANDOM_TEST R8BUT_RANDOM randomizes an R8BUT matrix. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 0 1 2 3 4 Row --- 0 0.964098 0.940507 0.122789 0.0560165 1 0.00771082 0.348098 0.443554 0.233051 2 0.166513 0.798095 0.932516 3 0.0499009 0.0405762 4 0.392951 R8BUT_SL_TEST R8BUT_SL solves A*x=b, where A is an R8BUT matrix. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 0 1 2 3 4 Row --- 0 0.232806 0.958397 0.205816 0.659713 1 0.180575 0.161927 0.30565 0.144377 2 0.814533 0.3044 0.176997 3 0.125607 0.272578 4 0.535913 x: 0: 1 1: 2 2: 3 3: 4 4: 5 b=A*x: 0: 5.4059 1: 2.79142 2: 4.54618 3: 1.86532 4: 2.67957 x: 0: 1 1: 2 2: 3 3: 4 4: 5 R8BUT_SLT_TEST R8BUT_SLT solves A'*x=b, where A is an R8BUT matrix. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 0 1 2 3 4 Row --- 0 0.588384 0.386728 0.419935 0.0459829 1 0.788747 0.0188066 0.254397 0.940593 2 0.817679 0.863401 0.847593 3 0.120099 0.247437 4 0.361706 x: 0: 1 1: 2 2: 3 3: 4 4: 5 b=A'*x: 0: 0.588384 1: 1.96422 2: 2.91058 3: 3.62538 4: 7.22224 x: 0: 1 1: 2 2: 3 3: 4 4: 5 R8BUT_TO_R8GE_TEST R8BUT_TO_R8GE converts a matrix from R8BUT to R8GE format. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 0 1 2 3 4 Row --- 0 0.031289 0.211112 0.406836 0.652772 1 0.702814 0.873518 0.571116 0.0693745 2 0.364992 0.721305 0.121146 3 0.0278004 0.899286 4 0.300349 The R8GE matrix: [[0.03128896 0.21111228 0.40683556 0.65277201 0. ] [0. 0.70281388 0.87351826 0.57111616 0.06937446] [0. 0. 0.36499244 0.72130499 0.12114594] [0. 0. 0. 0.0278004 0.89928595] [0. 0. 0. 0. 0.30034939]] R8BUT_ZEROS_TEST R8BUT_ZEROS zeros an R8BUT matrix. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 0 1 2 3 4 Row --- 0 0 0 0 0 1 0 0 0 0 2 0 0 0 3 0 0 4 0 r8but_test(): Normal end of execution. Sat Aug 20 09:41:38 2022