Tue May 20 22:27:33 2025 r8but_test(): python version: 3.10.12 numpy version: 1.26.4 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.564757 0.26892 0.208591 0.00831846 1 0.832394 0.350849 0.647839 0.140016 2 0.728046 0.888992 0.124802 3 0.512583 0.571004 4 0.263163 Determinant = 0.046167660538702075 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.0562604 0.588834 0.00107308 0.123864 1 0.135467 0.0847931 0.156556 0.41184 2 0.621808 0.763597 0.847306 3 0.593295 0.237009 4 0.475166 x: 0: 1 1: 2 2: 3 3: 4 4: 5 b=A'*x: 0: 0.0562604 1: 0.859768 2: 2.03608 3: 5.10095 4: 6.68946 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.384468 0.723332 0.586759 0.440092 1 0.762643 0.771892 0.489023 0.81553 2 0.228782 0.173502 0.361187 3 0.407474 0.580773 4 0.333768 x: 0: 1 1: 2 2: 3 3: 4 4: 5 b=A*x: 0: 5.35178 1: 9.8747 2: 3.18629 3: 4.53376 4: 1.66884 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.483839 0.211789 0.480851 0.682421 1 0.622849 0.0926914 0.070023 0.293195 2 0.661412 0.812909 0.112284 3 0.338961 0.206927 4 0.340598 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.257963 0.46037 0.364884 0.694092 1 0.990152 0.353072 0.277285 0.180874 2 0.228295 0.612594 0.112759 3 0.62159 0.489789 4 0.538611 x: 0: 1 1: 2 2: 3 3: 4 4: 5 b=A*x: 0: 5.04972 1: 5.05303 2: 3.69906 3: 4.9353 4: 2.69305 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.0507712 0.373985 0.610148 0.702158 1 0.555815 0.681253 0.494893 0.0439434 2 0.813305 0.229241 0.882908 3 0.167297 0.226492 4 0.425257 x: 0: 1 1: 2 2: 3 3: 4 4: 5 b=A'*x: 0: 0.0507712 1: 1.48561 2: 4.41257 3: 3.04886 4: 5.76886 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.526455 0.380202 0.499949 0.272868 1 0.254441 0.963502 0.131626 0.268056 2 0.339614 0.737813 0.577398 3 0.845006 0.95937 4 0.0163188 The R8GE matrix: [[0.52645534 0.38020204 0.49994919 0.27286838 0. ] [0. 0.2544415 0.96350158 0.13162646 0.26805628] [0. 0. 0.33961413 0.73781258 0.57739833] [0. 0. 0. 0.84500606 0.95937019] [0. 0. 0. 0. 0.01631879]] 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. Tue May 20 22:27:33 2025