20 January 2017 11:12:03 AM R8BUT_PRB C version Test the R8BUT library. R8BUT_DET 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.218418 0.956318 0.561695 0.257578 1 0.829509 0.415307 0.109957 0.0617272 2 0.0661187 0.043829 0.449539 3 0.633966 0.401306 4 0.754673 Determinant = 0.00573139 R8BUT_INDICATOR_TEST R8BUT_INDICATOR sets 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.218418 0.956318 0.561695 0.257578 1 0.829509 0.415307 0.109957 0.0617272 2 0.0661187 0.043829 0.449539 3 0.633966 0.401306 4 0.754673 x: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 b=A'*x: 0 0.218418 1 2.615336 2 1.590666 3 3.144841 4 6.850664 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.218418 0.956318 0.561695 0.257578 1 0.829509 0.415307 0.109957 0.0617272 2 0.0661187 0.043829 0.449539 3 0.633966 0.401306 4 0.754673 x: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 b=A*x: 0 4.846451 1 3.653403 2 2.621367 3 4.542394 4 3.773367 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 = 5 Upper bandwidth MU = 3 Rows 1:4, Cols 2:4: Col: 2 3 4 Row --- 2 33 34 35 3 44 45 4 55 R8BUT_RANDOM R8BUT_RANDOM sets a random R8BUT matrix. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 0 1 2 3 4 Row --- 0 0.218418 0.956318 0.561695 0.257578 1 0.829509 0.415307 0.109957 0.0617272 2 0.0661187 0.043829 0.449539 3 0.633966 0.401306 4 0.754673 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.218418 0.956318 0.561695 0.257578 1 0.829509 0.415307 0.109957 0.0617272 2 0.0661187 0.043829 0.449539 3 0.633966 0.401306 4 0.754673 x: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 b=A*x 0 4.846451 1 3.653403 2 2.621367 3 4.542394 4 3.773367 x: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 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.218418 0.956318 0.561695 0.257578 1 0.829509 0.415307 0.109957 0.0617272 2 0.0661187 0.043829 0.449539 3 0.633966 0.401306 4 0.754673 x: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 b=A'*x 0 0.218418 1 2.615336 2 1.590666 3 3.144841 4 6.850664 x: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 R8BUT_TO_R8GE 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.218418 0.956318 0.561695 0.257578 1 0.829509 0.415307 0.109957 0.0617272 2 0.0661187 0.043829 0.449539 3 0.633966 0.401306 4 0.754673 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.956318 0.561695 0.257578 0 2 0 0.829509 0.415307 0.109957 0.0617272 3 0 0 0.0661187 0.043829 0.449539 4 0 0 0 0.633966 0.401306 5 0 0 0 0 0.754673 R8BUT_ZEROS 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_PRB Normal end of execution. 20 January 2017 11:12:03 AM