18 August 2022 11:46:48 AM r8blt_test(): C version Test r8blt(). 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: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 5 6.000000 6 7.000000 7 8.000000 8 9.000000 9 10.000000 b = A'*x: 0 4.846451 1 3.653403 2 7.405089 3 5.215238 4 10.238541 5 9.201267 6 18.053679 7 20.322480 8 11.016348 9 4.537938 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: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 5 6.000000 6 7.000000 7 8.000000 8 9.000000 9 10.000000 b = A*x: 0 0.218418 1 2.615336 2 1.590666 3 3.144841 4 6.850664 5 8.991249 6 11.486777 7 9.782584 8 12.130674 9 21.626645 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 for 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: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 5 6.000000 6 7.000000 7 8.000000 8 9.000000 9 10.000000 The right hand side: 0 0.218418 1 2.615336 2 1.590666 3 3.144841 4 6.850664 5 8.991249 6 11.486777 7 9.782584 8 12.130674 9 21.626645 Solution to A*x=b: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 5 6.000000 6 7.000000 7 8.000000 8 9.000000 9 10.000000 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 x: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 5 6.000000 6 7.000000 7 8.000000 8 9.000000 9 10.000000 The right hand side: 0 4.846451 1 3.653403 2 7.405089 3 5.215238 4 10.238541 5 9.201267 6 18.053679 7 20.322480 8 11.016348 9 4.537938 Solution to A'*x=b: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 5 6.000000 6 7.000000 7 8.000000 8 9.000000 9 10.000000 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 = 10 Lower bandwidth ML = 3 The R8BLT matrix: Col: 1 2 3 4 5 Row --- 1 0 2 0 0 3 0 0 0 4 0 0 0 0 5 0 0 0 0 6 0 0 0 7 0 0 8 0 Col: 6 7 8 9 10 Row --- 6 0 7 0 0 8 0 0 0 9 0 0 0 0 10 0 0 0 0 r8blt(): Normal end of execution. 18 August 2022 11:46:48 AM