13 October 2015 08:35:19 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: 1 2 3 4 5 Row --- 1 0.633966 0.0661187 0.829509 0.218418 2 0.0617272 0.257578 0.561695 0.956318 3 0.449539 0.109957 0.415307 4 0.401306 0.043829 5 0.754673 Determinant = 0.00532776 R8BUT_INDICATOR_TEST R8BUT_INDICATOR sets an R8BUT indicator matrix. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 11 12 13 14 2 22 23 24 25 3 33 34 35 4 44 45 5 55 R8BUT_MTV_TEST R8BUT_MTV computes b=A'*x, where A is an R8BUT matrix; Matrix order N = 10 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 0.260303 0.00183837 0.109957 0.218418 2 0.912484 0.897504 0.043829 0.956318 3 0.113664 0.350752 0.633966 4 0.351629 0.0945448 5 0.822887 6 7 8 Col: 6 7 8 9 10 Row --- 6 0.267132 0.859097 0.401306 0.0661187 7 0.692066 0.840847 0.754673 0.257578 8 0.561662 0.123104 0.797287 9 0.861216 0.00751236 10 0.453794 x: 1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 10 10 b=A'*x: 1 0.260303 2 1.82681 3 2.24596 4 2.76485 5 8.30715 6 4.40631 7 14.4935 8 14.8636 9 14.4152 10 12.7869 R8BUT_MV_TEST R8BUT_MV computes b=A*x, where A is an R8BUT matrix; Matrix order N = 10 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 0.260303 0.00183837 0.109957 0.218418 2 0.912484 0.897504 0.043829 0.956318 3 0.113664 0.350752 0.633966 4 0.351629 0.0945448 5 0.822887 6 7 8 Col: 6 7 8 9 10 Row --- 6 0.267132 0.859097 0.401306 0.0661187 7 0.692066 0.840847 0.754673 0.257578 8 0.561662 0.123104 0.797287 9 0.861216 0.00751236 10 0.453794 x: 1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 10 10 b=A*x: 1 1.46752 2 9.47438 3 9.89089 4 6.18147 5 10.6654 6 11.422 7 20.9391 8 13.5741 9 7.82607 10 4.53794 R8BUT_PRINT_TEST R8BUT_PRINT prints an R8BUT matrix. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 11 12 13 14 2 22 23 24 25 3 33 34 35 4 44 45 5 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 22 23 24 3 33 34 4 44 R8BUT_RANDOM R8BUT_RANDOM sets a random R8BUT matrix. Matrix order N = 5 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 0.633966 0.0661187 0.829509 0.218418 2 0.0617272 0.257578 0.561695 0.956318 3 0.449539 0.109957 0.415307 4 0.401306 0.043829 5 0.754673 R8BUT_SL_TEST R8BUT_SL solves A*x=b, where A is an R8BUT matrix; Matrix order N = 10 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 0.260303 0.00183837 0.109957 0.218418 2 0.912484 0.897504 0.043829 0.956318 3 0.113664 0.350752 0.633966 4 0.351629 0.0945448 5 0.822887 6 7 8 Col: 6 7 8 9 10 Row --- 6 0.267132 0.859097 0.401306 0.0661187 7 0.692066 0.840847 0.754673 0.257578 8 0.561662 0.123104 0.797287 9 0.861216 0.00751236 10 0.453794 x: 1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 10 10 b=A*x 1 1.46752 2 9.47438 3 9.89089 4 6.18147 5 10.6654 6 11.422 7 20.9391 8 13.5741 9 7.82607 10 4.53794 x: 1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 10 10 R8BUT_SLT_TEST R8BUT_SLT solves A'*x=b, where A is an R8BUT matrix; Matrix order N = 10 Upper bandwidth MU = 3 The R8BUT matrix: Col: 1 2 3 4 5 Row --- 1 0.260303 0.00183837 0.109957 0.218418 2 0.912484 0.897504 0.043829 0.956318 3 0.113664 0.350752 0.633966 4 0.351629 0.0945448 5 0.822887 6 7 8 Col: 6 7 8 9 10 Row --- 6 0.267132 0.859097 0.401306 0.0661187 7 0.692066 0.840847 0.754673 0.257578 8 0.561662 0.123104 0.797287 9 0.861216 0.00751236 10 0.453794 x: 1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 10 10 b=A'*x 1 0.260303 2 1.82681 3 2.24596 4 2.76485 5 8.30715 6 4.40631 7 14.4935 8 14.8636 9 14.4152 10 12.7869 x: 1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 10 10 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: 1 2 3 4 5 Row --- 1 0.633966 0.0661187 0.829509 0.218418 2 0.0617272 0.257578 0.561695 0.956318 3 0.449539 0.109957 0.415307 4 0.401306 0.043829 5 0.754673 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 0.633966 0.0661187 0.829509 0.218418 0 2 0 0.0617272 0.257578 0.561695 0.956318 3 0 0 0.449539 0.109957 0.415307 4 0 0 0 0.401306 0.043829 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: 1 2 3 4 5 Row --- 1 0 0 0 0 2 0 0 0 0 3 0 0 0 4 0 0 5 0 R8BUT_PRB Normal end of execution. 13 October 2015 08:35:19 AM