02 July 2016 9:27:30.658 AM R8LT_PRB FORTRAN90 version: Test the R8LT library. R8GE_TO_R8LT_TEST R8GE_TO_R8LT converts an R8GE matrix to R8LT format. The random R8GE matrix: Col: 1 2 3 4 Row --- 1 0.218418 0.661187E-01 0.617272E-01 0.183837E-02 2 0.956318 0.257578 0.449539 0.897504 3 0.829509 0.109957 0.401306 0.350752 4 0.561695 0.438290E-01 0.754673 0.945448E-01 5 0.415307 0.633966 0.797287 0.136169E-01 The R8LT matrix Col: 1 2 3 4 Row --- 1 0.218418 2 0.956318 0.257578 3 0.829509 0.109957 0.401306 4 0.561695 0.438290E-01 0.754673 0.945448E-01 5 0.415307 0.633966 0.797287 0.136169E-01 R8LT_DET_TEST R8LT_DET computes the determinant of an R8LT matrix. The matrix A: Col: 1 2 3 4 5 Row --- 1 0.218418 2 0.956318 0.661187E-01 3 0.829509 0.257578 0.633966 4 0.561695 0.109957 0.617272E-01 0.401306 5 0.415307 0.438290E-01 0.449539 0.754673 0.797287 Determinant is 0.292934E-02 R8LT_INDICATOR_TEST R8LT_INDICATOR sets up an indicator matrix in R8LT format; The R8LT 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 41.0000 42.0000 43.0000 44.0000 5 51.0000 52.0000 53.0000 54.0000 55.0000 6 61.0000 62.0000 63.0000 64.0000 65.0000 R8LT_INVERSE_TEST R8LT_INVERSE computes the inverse of an R8LT matrix. The matrix A: Col: 1 2 3 4 5 Row --- 1 0.218418 2 0.956318 0.661187E-01 3 0.829509 0.257578 0.633966 4 0.561695 0.109957 0.617272E-01 0.401306 5 0.415307 0.438290E-01 0.449539 0.754673 0.797287 The inverse matrix B: Col: 1 2 3 4 5 Row --- 1 4.57837 2 -66.2199 15.1243 3 20.9143 -6.14495 1.57737 4 8.51892 -3.19883 -0.242625 2.49186 5 -18.6004 5.66117 -0.659722 -2.35868 1.25425 The product C = A * B: Col: 1 2 3 4 5 Row --- 1 1.00000 2 0.00000 1.00000 3 0.00000 0.00000 1.00000 4 0.00000 0.00000 0.00000 1.00000 5 0.00000 0.00000 0.00000 0.222045E-15 1.00000 R8LT_MM_TEST R8LT_MM computes C = A * B for R8LT matrices. The matrix A: Col: 1 2 3 4 5 Row --- 1 1.00000 2 1.00000 1.00000 3 1.00000 1.00000 1.00000 4 1.00000 1.00000 1.00000 1.00000 5 1.00000 1.00000 1.00000 1.00000 1.00000 The product C = A * A Col: 1 2 3 4 5 Row --- 1 1.00000 2 2.00000 1.00000 3 3.00000 2.00000 1.00000 4 4.00000 3.00000 2.00000 1.00000 5 5.00000 4.00000 3.00000 2.00000 1.00000 R8LT_MTM_TEST R8LT_MTM computes C = A' * B for R8LT matrices. The matrix A: Col: 1 2 3 4 5 Row --- 1 1.00000 2 1.00000 1.00000 3 1.00000 1.00000 1.00000 4 1.00000 1.00000 1.00000 1.00000 5 1.00000 1.00000 1.00000 1.00000 1.00000 The product C = A' * A Col: 1 2 3 4 5 Row --- 1 5.00000 4.00000 3.00000 2.00000 1.00000 2 4.00000 4.00000 3.00000 2.00000 1.00000 3 3.00000 3.00000 3.00000 2.00000 1.00000 4 2.00000 2.00000 2.00000 2.00000 1.00000 5 1.00000 1.00000 1.00000 1.00000 1.00000 R8LT_MTV_TEST R8LT_MTV computes A'*x for an R8LT matrix. The matrix A: Col: 1 2 3 4 Row --- 1 11.0000 2 21.0000 22.0000 3 31.0000 32.0000 33.0000 4 41.0000 42.0000 43.0000 44.0000 5 51.0000 52.0000 53.0000 54.0000 The vector x 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 b = A'*x: 1 565.000 2 568.000 3 536.000 4 446.000 R8LT_MV_TEST R8LT_MV computes A*x for an R8LT matrix. The matrix A: Col: 1 2 3 4 Row --- 1 11.0000 2 21.0000 22.0000 3 31.0000 32.0000 33.0000 4 41.0000 42.0000 43.0000 44.0000 5 51.0000 52.0000 53.0000 54.0000 The vector x 1 1.00000 2 2.00000 3 3.00000 4 4.00000 b = A*x: 1 11.0000 2 65.0000 3 194.000 4 430.000 5 530.000 R8LT_PRINT_TEST R8UT_PRINT prints an R8LT matrix. Matrix rows M = 6 Matrix columns N = 4 The R8UT matrix: Col: 1 2 3 4 Row --- 1 11.0000 2 21.0000 22.0000 3 31.0000 32.0000 33.0000 4 41.0000 42.0000 43.0000 44.0000 5 51.0000 52.0000 53.0000 54.0000 6 61.0000 62.0000 63.0000 64.0000 R8LT_PRINT_SOME_TEST R8LT_PRINT_SOME prints some of an R8LT matrix. Matrix rows M = 4 Matrix columns N = 6 Rows 2-3, Cols 2-4: Col: 2 3 4 Row --- 2 22.0000 3 32.0000 33.0000 R8LT_RANDOM_TEST R8LT_RANDOM randomizes an R8LT matrix. The random RLT matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 2 0.956318 0.257578 3 0.829509 0.109957 0.449539 4 0.561695 0.438290E-01 0.401306 0.183837E-02 5 0.415307 0.633966 0.754673 0.897504 0.945448E-01 6 0.661187E-01 0.617272E-01 0.797287 0.350752 0.136169E-01 R8LT_SL_TEST R8LT_SL solves A*x=b for an R8LT matrix A. Matrix order N = 5 The R8LT matrix: Col: 1 2 3 4 5 Row --- 1 1.00000 2 1.00000 2.00000 3 1.00000 2.00000 3.00000 4 1.00000 2.00000 3.00000 4.00000 5 1.00000 2.00000 3.00000 4.00000 5.00000 Solution: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 R8LT_SLT_TEST R8LT_SLT solves A'*x=b for an R8LT matrix A Matrix order N = 5 The R8LT matrix: Col: 1 2 3 4 5 Row --- 1 1.00000 2 1.00000 2.00000 3 1.00000 2.00000 3.00000 4 1.00000 2.00000 3.00000 4.00000 5 1.00000 2.00000 3.00000 4.00000 5.00000 Solution to transposed system: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 R8LT_TO_R8GE_TEST R8LT_TO_R8GE converts an R8LT matrix to R8GE format. The random R8LT matrix: Col: 1 2 3 4 Row --- 1 0.218418 2 0.956318 0.661187E-01 3 0.829509 0.257578 0.633966 4 0.561695 0.109957 0.617272E-01 0.401306 5 0.415307 0.438290E-01 0.449539 0.754673 The R8GE matrix Col: 1 2 3 4 Row --- 1 0.218418 0.00000 0.00000 0.00000 2 0.956318 0.661187E-01 0.00000 0.00000 3 0.829509 0.257578 0.633966 0.00000 4 0.561695 0.109957 0.617272E-01 0.401306 5 0.415307 0.438290E-01 0.449539 0.754673 R8LT_ZEROS_TEST R8LT_ZEROS zeros an R8LT matrix; The R8LT 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 0.00000 6 0.00000 0.00000 0.00000 0.00000 0.00000 R8LT_PRB Normal end of execution. 02 July 2016 9:27:30.659 AM