03 July 2016 08:41:50 PM R8LT_PRB C++ 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.0661187 0.0617272 0.00183837 2 0.956318 0.257578 0.449539 0.897504 3 0.829509 0.109957 0.401306 0.350752 4 0.561695 0.043829 0.754673 0.0945448 5 0.415307 0.633966 0.797287 0.0136169 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.043829 0.754673 0.0945448 5 0.415307 0.633966 0.797287 0.0136169 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.0661187 3 0.829509 0.257578 0.633966 4 0.561695 0.109957 0.0617272 0.401306 5 0.415307 0.043829 0.449539 0.754673 0.797287 Determinant is 0.00292934 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 2 21 22 3 31 32 33 4 41 42 43 44 5 51 52 53 54 55 6 61 62 63 64 65 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.0661187 3 0.829509 0.257578 0.633966 4 0.561695 0.109957 0.0617272 0.401306 5 0.415307 0.043829 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 A * B: Col: 1 2 3 4 5 Row --- 1 1 2 0 1 3 0 0 1 4 0 0 0 1 5 0 0 0 2.22045e-16 1 R8LT_MM_TEST R8LT_MM computes C=A*B for R8LT matrices. The matrix: Col: 1 2 3 4 5 Row --- 1 1 2 1 1 3 1 1 1 4 1 1 1 1 5 1 1 1 1 1 The product C = A * A: Col: 1 2 3 4 5 Row --- 1 1 2 2 1 3 3 2 1 4 4 3 2 1 5 5 4 3 2 1 R8LT_MTM_TEST R8LT_MTM computes C=A'*B for R8LT matrices. The matrix: Col: 1 2 3 4 5 Row --- 1 1 2 1 1 3 1 1 1 4 1 1 1 1 5 1 1 1 1 1 The product C = A' * A: Col: 1 2 3 4 5 Row --- 1 5 4 3 2 1 2 4 4 3 2 1 3 3 3 3 2 1 4 2 2 2 2 1 5 1 1 1 1 1 R8LT_MTV_TEST R8LT_MTV computes a product b=A'*x where A is in R8LT format. The matrix: Col: 1 2 3 4 Row --- 1 11 2 21 22 3 31 32 33 4 41 42 43 44 5 51 52 53 54 Vector x: 1 1 2 2 3 3 4 4 5 5 Vector b: 1 565 2 568 3 536 4 446 R8LT_MV_TEST R8LT_MV computes a product b=A*x where A is in R8LT format. The matrix: Col: 1 2 3 4 Row --- 1 11 2 21 22 3 31 32 33 4 41 42 43 44 5 51 52 53 54 Vector x: 1 1 2 2 3 3 4 4 Vector b: 1 11 2 65 3 194 4 430 5 530 R8LT_PRINT_TEST R8LT_PRINT prints an R8LT marix. The matrix: Col: 1 2 3 4 Row --- 1 11 2 21 22 3 31 32 33 4 41 42 43 44 5 51 52 53 54 6 61 62 63 64 R8LT_PRINT_SOME_TEST R8LT_PRINT_SOME prints some of an R8LT marix. Rows 1-2, Cols 1-3: Col: 1 2 3 Row --- 1 11 2 21 22 R8LT_RANDOM_TEST R8LT_RANDOM randomizes a matrix in R8LT format; The random R8LT 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.043829 0.401306 0.00183837 5 0.415307 0.633966 0.754673 0.897504 0.0945448 6 0.0661187 0.0617272 0.797287 0.350752 0.0136169 R8LT_SL_TEST R8LT_SL solves a linear system A*x=b with R8LT matrix; The R8LT matrix: Col: 1 2 3 4 5 Row --- 1 1 2 2 2 3 3 3 3 4 4 4 4 4 5 5 5 5 5 5 Solution: 1 1 2 2 3 3 4 4 5 5 R8LT_SLT_TEST R8LT_SLT solves a linear system A'*x=b with R8LT matrix; The R8LT matrix: Col: 1 2 3 4 5 Row --- 1 1 2 2 2 3 3 3 3 4 4 4 4 4 5 5 5 5 5 5 Solution to transposed system: 1 1 2 2 3 3 4 4 5 5 R8LT_TO_R8GE_TEST R8LT_TO_R8GE converts an R8LT matrix to R8GE format. R8LT matrix: Col: 1 2 3 4 Row --- 1 0.218418 2 0.956318 0.0661187 3 0.829509 0.257578 0.633966 4 0.561695 0.109957 0.0617272 0.401306 5 0.415307 0.043829 0.449539 0.754673 R8GE matrix Col: 1 2 3 4 Row --- 1 0.218418 0 0 0 2 0.956318 0.0661187 0 0 3 0.829509 0.257578 0.633966 0 4 0.561695 0.109957 0.0617272 0.401306 5 0.415307 0.043829 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 2 0 0 3 0 0 0 4 0 0 0 0 5 0 0 0 0 0 6 0 0 0 0 0 R8LT_PRB Normal end of execution. 03 July 2016 08:41:50 PM