06 October 2015 10:26:03.821 AM R8BLT_PRB FORTRAN90 version: Test the R8BLT library. R8BLT_DET_TEST R8BLT_DET computes 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.661187E-01 4 0.257578 0.109957 0.438290E-01 0.633966 5 0.617272E-01 0.449539 0.401306 0.754673 6 0.797287 0.183837E-02 0.897504 7 0.945448E-01 0.136169E-01 8 0.123104 Col: 6 7 8 9 10 Row --- 6 0.350752 7 0.859097 0.840847 8 0.751236E-02 0.260303 0.912484 9 0.113664 0.351629 0.822887 0.267132 10 0.692066 0.561662 0.861216 0.453794 Determinant = 0.186977E-03 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.0000 2 21.0000 22.0000 3 31.0000 32.0000 33.0000 4 42.0000 43.0000 44.0000 5 53.0000 54.0000 55.0000 6 64.0000 65.0000 Col: 6 Row --- 6 66.0000 R8BLT_MTV_TEST R8BLT_MTV computes b=A'*x, where A is an R8BLT matrix. Matrix order N = 5 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.661187E-01 4 0.257578 0.109957 0.438290E-01 0.633966 5 0.617272E-01 0.449539 0.401306 0.754673 x: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 b=A'*x: 1 4.84645 2 3.65340 3 2.62137 4 4.54239 5 3.77337 R8BLT_MV_TEST R8BLT_MV computes b=A*x, where A is an R8BLT matrix. Matrix order N = 5 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.661187E-01 4 0.257578 0.109957 0.438290E-01 0.633966 5 0.617272E-01 0.449539 0.401306 0.754673 x: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 b=A*x: 1 0.218418 2 2.61534 3 1.59067 4 3.14484 5 6.85066 R8BLT_PRINT_TEST R8BLT_PRINT prints 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.661187E-01 4 0.257578 0.109957 0.438290E-01 0.633966 5 0.617272E-01 0.449539 0.401306 0.754673 6 0.797287 0.183837E-02 0.897504 7 0.945448E-01 0.136169E-01 8 0.123104 Col: 6 7 8 9 10 Row --- 6 0.350752 7 0.859097 0.840847 8 0.751236E-02 0.260303 0.912484 9 0.113664 0.351629 0.822887 0.267132 10 0.692066 0.561662 0.861216 0.453794 R8BLT_PRINT_SOME_TEST R8BLT_PRINT_SOME prints some of an R8BLT matrix. Matrix order N = 10 Lower bandwidth ML = 3 Rows 2:5, Cols 3:5: Col: 3 4 5 Row --- 3 303.000 4 403.000 404.000 5 503.000 504.000 505.000 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.661187E-01 4 0.257578 0.109957 0.438290E-01 0.633966 5 0.617272E-01 0.449539 0.401306 0.754673 6 0.797287 0.183837E-02 0.897504 7 0.945448E-01 0.136169E-01 8 0.123104 Col: 6 7 8 9 10 Row --- 6 0.350752 7 0.859097 0.840847 8 0.751236E-02 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.661187E-01 4 0.257578 0.109957 0.438290E-01 0.633966 5 0.617272E-01 0.449539 0.401306 0.754673 6 0.797287 0.183837E-02 0.897504 7 0.945448E-01 0.136169E-01 8 0.123104 Col: 6 7 8 9 10 Row --- 6 0.350752 7 0.859097 0.840847 8 0.751236E-02 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: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 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.661187E-01 4 0.257578 0.109957 0.438290E-01 0.633966 5 0.617272E-01 0.449539 0.401306 0.754673 6 0.797287 0.183837E-02 0.897504 7 0.945448E-01 0.136169E-01 8 0.123104 Col: 6 7 8 9 10 Row --- 6 0.350752 7 0.859097 0.840847 8 0.751236E-02 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.65340 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: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 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.661187E-01 4 0.257578 0.109957 0.438290E-01 0.633966 5 0.617272E-01 0.449539 0.401306 0.754673 6 0.797287 0.183837E-02 0.897504 7 0.945448E-01 0.136169E-01 8 0.123104 Col: 6 7 8 9 10 Row --- 6 0.350752 7 0.859097 0.840847 8 0.751236E-02 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.00000 0.00000 0.00000 0.00000 2 0.956318 0.829509 0.00000 0.00000 0.00000 3 0.561695 0.415307 0.661187E-01 0.00000 0.00000 4 0.257578 0.109957 0.438290E-01 0.633966 0.00000 5 0.00000 0.617272E-01 0.449539 0.401306 0.754673 6 0.00000 0.00000 0.797287 0.183837E-02 0.897504 7 0.00000 0.00000 0.00000 0.945448E-01 0.136169E-01 8 0.00000 0.00000 0.00000 0.00000 0.123104 9 0.00000 0.00000 0.00000 0.00000 0.00000 10 0.00000 0.00000 0.00000 0.00000 0.00000 Col: 6 7 8 9 10 Row --- 1 0.00000 0.00000 0.00000 0.00000 0.00000 2 0.00000 0.00000 0.00000 0.00000 0.00000 3 0.00000 0.00000 0.00000 0.00000 0.00000 4 0.00000 0.00000 0.00000 0.00000 0.00000 5 0.00000 0.00000 0.00000 0.00000 0.00000 6 0.350752 0.00000 0.00000 0.00000 0.00000 7 0.859097 0.840847 0.00000 0.00000 0.00000 8 0.751236E-02 0.260303 0.912484 0.00000 0.00000 9 0.113664 0.351629 0.822887 0.267132 0.00000 10 0.00000 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.00000 2 0.00000 0.00000 3 0.00000 0.00000 0.00000 4 0.00000 0.00000 0.00000 0.00000 5 0.00000 0.00000 0.00000 0.00000 6 0.00000 0.00000 0.00000 7 0.00000 0.00000 8 0.00000 Col: 6 7 8 9 10 Row --- 6 0.00000 7 0.00000 0.00000 8 0.00000 0.00000 0.00000 9 0.00000 0.00000 0.00000 0.00000 10 0.00000 0.00000 0.00000 0.00000 R8BLT_PRB Normal end of execution. 06 October 2015 10:26:03.834 AM