08 October 2015 08:26:54 AM R8BLT_PRB C++ version Test the R8BLT library. R8BLT_DET_TEST R8BLT_DET gets the determinant of an R8BLT matrix; Matrix order N = 10 Lower bandwidth ML = 3 The R8BLT matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 2 0.956318 0.829509 3 0.561695 0.415307 0.0661187 4 0.257578 0.109957 0.043829 0.633966 5 0.0617272 0.449539 0.401306 0.754673 6 0.797287 0.00183837 0.897504 7 0.0945448 0.0136169 8 0.123104 Col: 6 7 8 9 10 Row --- 6 0.350752 7 0.859097 0.840847 8 0.00751236 0.260303 0.912484 9 0.113664 0.351629 0.822887 0.267132 10 0.692066 0.561662 0.861216 0.453794 Determinant = 0.000186977 R8BLT_INDICATOR_TEST R8BLT_INDICATOR sets up an R8BLT indicator matrix. Matrix order N = 6 Lower bandwidth ML = 2 The R8BLT indicator matrix: Col: 1 2 3 4 5 Row --- 1 11 2 21 22 3 31 32 33 4 42 43 44 5 53 54 55 6 64 65 Col: 6 Row --- 6 66 R8BLT_MTV_TEST R8BLT_MTV computes A'*x, where A is an R8BLT matrix; Matrix order N = 10 Lower bandwidth ML = 3 The R8BLT matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 2 0.956318 0.829509 3 0.561695 0.415307 0.0661187 4 0.257578 0.109957 0.043829 0.633966 5 0.0617272 0.449539 0.401306 0.754673 6 0.797287 0.00183837 0.897504 7 0.0945448 0.0136169 8 0.123104 Col: 6 7 8 9 10 Row --- 6 0.350752 7 0.859097 0.840847 8 0.00751236 0.260303 0.912484 9 0.113664 0.351629 0.822887 0.267132 10 0.692066 0.561662 0.861216 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 4.84645 2 3.6534 3 7.40509 4 5.21524 5 10.2385 6 9.20127 7 18.0537 8 20.3225 9 11.0163 10 4.53794 R8BLT_MV_TEST R8BLT_MV computes A*x, where A is an R8BLT matrix; Matrix order N = 10 Lower bandwidth ML = 3 The R8BLT matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 2 0.956318 0.829509 3 0.561695 0.415307 0.0661187 4 0.257578 0.109957 0.043829 0.633966 5 0.0617272 0.449539 0.401306 0.754673 6 0.797287 0.00183837 0.897504 7 0.0945448 0.0136169 8 0.123104 Col: 6 7 8 9 10 Row --- 6 0.350752 7 0.859097 0.840847 8 0.00751236 0.260303 0.912484 9 0.113664 0.351629 0.822887 0.267132 10 0.692066 0.561662 0.861216 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.218418 2 2.61534 3 1.59067 4 3.14484 5 6.85066 6 8.99125 7 11.4868 8 9.78258 9 12.1307 10 21.6266 R8BLT_PRINT_TEST R8BLT_PRINT prints an R8BLT matrix. Matrix order N = 6 Lower bandwidth ML = 2 The R8BLT matrix: Col: 1 2 3 4 5 Row --- 1 11 2 21 22 3 31 32 33 4 42 43 44 5 53 54 55 6 64 65 Col: 6 Row --- 6 66 R8BLT_PRINT_SOME_TEST R8BLT_PRINT_SOME prints some of an R8BLT matrix. Matrix order N = 10 Lower bandwidth ML = 3 Rows 1:4, Cols 2:4: Col: 2 3 4 Row --- 2 202 3 302 303 4 402 403 404 R8BLT_RANDOM_TEST R8BLT_RANDOM randomizes an R8BLT matrix. Matrix order N = 10 Lower bandwidth ML = 3 The R8BLT matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 2 0.956318 0.829509 3 0.561695 0.415307 0.0661187 4 0.257578 0.109957 0.043829 0.633966 5 0.0617272 0.449539 0.401306 0.754673 6 0.797287 0.00183837 0.897504 7 0.0945448 0.0136169 8 0.123104 Col: 6 7 8 9 10 Row --- 6 0.350752 7 0.859097 0.840847 8 0.00751236 0.260303 0.912484 9 0.113664 0.351629 0.822887 0.267132 10 0.692066 0.561662 0.861216 0.453794 R8BLT_SL_TEST R8BLT_SL solves A*x=b, where A is an R8BLT matrix. Matrix order N = 10 Lower bandwidth ML = 3 The R8BLT matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 2 0.956318 0.829509 3 0.561695 0.415307 0.0661187 4 0.257578 0.109957 0.043829 0.633966 5 0.0617272 0.449539 0.401306 0.754673 6 0.797287 0.00183837 0.897504 7 0.0945448 0.0136169 8 0.123104 Col: 6 7 8 9 10 Row --- 6 0.350752 7 0.859097 0.840847 8 0.00751236 0.260303 0.912484 9 0.113664 0.351629 0.822887 0.267132 10 0.692066 0.561662 0.861216 0.453794 The right hand side: 1 0.218418 2 2.61534 3 1.59067 4 3.14484 5 6.85066 6 8.99125 7 11.4868 8 9.78258 9 12.1307 10 21.6266 Solution to A*x=b: 1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 10 10 R8BLT_SLT_TEST R8BLT_SLT solves A'*x=b, where A is an R8BLT matrix; Matrix order N = 10 Lower bandwidth ML = 3 The R8BLT matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 2 0.956318 0.829509 3 0.561695 0.415307 0.0661187 4 0.257578 0.109957 0.043829 0.633966 5 0.0617272 0.449539 0.401306 0.754673 6 0.797287 0.00183837 0.897504 7 0.0945448 0.0136169 8 0.123104 Col: 6 7 8 9 10 Row --- 6 0.350752 7 0.859097 0.840847 8 0.00751236 0.260303 0.912484 9 0.113664 0.351629 0.822887 0.267132 10 0.692066 0.561662 0.861216 0.453794 The right hand side: 1 4.84645 2 3.6534 3 7.40509 4 5.21524 5 10.2385 6 9.20127 7 18.0537 8 20.3225 9 11.0163 10 4.53794 Solution to A'*x=b: 1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 10 10 R8BLT_TO_R8GE_TEST R8BLT_TO_R8GE converts a matrix from R8BLT to R8GE format; Matrix order N = 10 Lower bandwidth ML = 3 The R8BLT matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 2 0.956318 0.829509 3 0.561695 0.415307 0.0661187 4 0.257578 0.109957 0.043829 0.633966 5 0.0617272 0.449539 0.401306 0.754673 6 0.797287 0.00183837 0.897504 7 0.0945448 0.0136169 8 0.123104 Col: 6 7 8 9 10 Row --- 6 0.350752 7 0.859097 0.840847 8 0.00751236 0.260303 0.912484 9 0.113664 0.351629 0.822887 0.267132 10 0.692066 0.561662 0.861216 0.453794 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0 0 0 0 2 0.956318 0.829509 0 0 0 3 0.561695 0.415307 0.0661187 0 0 4 0.257578 0.109957 0.043829 0.633966 0 5 0 0.0617272 0.449539 0.401306 0.754673 6 0 0 0.797287 0.00183837 0.897504 7 0 0 0 0.0945448 0.0136169 8 0 0 0 0 0.123104 9 0 0 0 0 0 10 0 0 0 0 0 Col: 6 7 8 9 10 Row --- 1 0 0 0 0 0 2 0 0 0 0 0 3 0 0 0 0 0 4 0 0 0 0 0 5 0 0 0 0 0 6 0.350752 0 0 0 0 7 0.859097 0.840847 0 0 0 8 0.00751236 0.260303 0.912484 0 0 9 0.113664 0.351629 0.822887 0.267132 0 10 0 0.692066 0.561662 0.861216 0.453794 R8BLT_ZEROS_TEST R8BLT_ZEROS zeros an R8BLT matrix. Matrix order N = 6 Lower bandwidth ML = 2 The R8BLT matrix: Col: 1 2 3 4 5 Row --- 1 0 2 0 0 3 0 0 0 4 0 0 0 5 0 0 0 6 0 0 Col: 6 Row --- 6 0 R8BLT_PRB Normal end of execution. 08 October 2015 08:26:54 AM