02 July 2016 9:27:48.349 AM R8UT_PRB FORTRAN90 version: Test the R8UT library. R8GE_TO_R8UT_TEST R8GE_TO_R8UT converts an R8GE matrix to R8UT 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 R8UT matrix Col: 1 2 3 4 Row --- 1 0.218418 0.661187E-01 0.617272E-01 0.183837E-02 2 0.257578 0.449539 0.897504 3 0.401306 0.350752 4 0.945448E-01 5 R8UT_DET_TEST R8UT_DET computes the determinant of an R8UT matrix. The matrix A: Col: 1 2 3 4 5 Row --- 1 0.218418 0.956318 0.561695 0.257578 0.617272E-01 2 0.829509 0.415307 0.109957 0.449539 3 0.661187E-01 0.438290E-01 0.401306 4 0.633966 0.754673 5 0.797287 Determinant is 0.605501E-02 R8UT_INDICATOR_TEST R8UT_INDICATOR sets up an indicator matrix. Matrix rows M = 5 Matrix columns N = 4 The R8UT indicator matrix: Col: 1 2 3 4 Row --- 1 11.0000 12.0000 13.0000 14.0000 2 22.0000 23.0000 24.0000 3 33.0000 34.0000 4 44.0000 5 R8UT_INVERSE_TEST R8UT_INVERSE computes the inverse of an R8UT matrix. The matrix A: Col: 1 2 3 4 5 Row --- 1 0.218418 0.956318 0.561695 0.257578 0.617272E-01 2 0.829509 0.415307 0.109957 0.449539 3 0.661187E-01 0.438290E-01 0.401306 4 0.633966 0.754673 5 0.797287 The inverse matrix B: Col: 1 2 3 4 5 Row --- 1 4.57837 -5.27827 -5.74037 -0.547837 6.02953 2 1.20553 -7.57223 0.314412 2.83407 3 15.1243 -1.04561 -6.62294 4 1.57737 -1.49306 5 1.25425 The product C = A * B: Col: 1 2 3 4 5 Row --- 1 1.00000 0.00000 0.00000 0.00000 -0.208167E-15 2 1.00000 0.00000 0.277556E-16 -0.111022E-15 3 1.00000 0.00000 0.00000 4 1.00000 0.00000 5 1.00000 R8UT_MM_TEST R8UT_MM computes C = A * B for R8UT matrices. The matrix A: Col: 1 2 3 4 5 Row --- 1 1.00000 1.00000 1.00000 1.00000 1.00000 2 1.00000 1.00000 1.00000 1.00000 3 1.00000 1.00000 1.00000 4 1.00000 1.00000 5 1.00000 The product C = A * A Col: 1 2 3 4 5 Row --- 1 1.00000 2.00000 3.00000 4.00000 5.00000 2 1.00000 2.00000 3.00000 4.00000 3 1.00000 2.00000 3.00000 4 1.00000 2.00000 5 1.00000 R8UT_MTM_TEST R8UT_MTM computes C = A' * B for R8UT matrices. The matrix A: Col: 1 2 3 4 5 Row --- 1 1.00000 1.00000 1.00000 1.00000 1.00000 2 1.00000 1.00000 1.00000 1.00000 3 1.00000 1.00000 1.00000 4 1.00000 1.00000 5 1.00000 The product C = A' * A Col: 1 2 3 4 5 Row --- 1 1.00000 1.00000 1.00000 1.00000 1.00000 2 1.00000 2.00000 2.00000 2.00000 2.00000 3 1.00000 2.00000 3.00000 3.00000 3.00000 4 1.00000 2.00000 3.00000 4.00000 4.00000 5 1.00000 2.00000 3.00000 4.00000 5.00000 R8UT_MTV_TEST R8UT_MTV computes A'*x for an R8UT matrix. The matrix A: Col: 1 2 3 4 Row --- 1 11.0000 12.0000 13.0000 14.0000 2 22.0000 23.0000 24.0000 3 33.0000 34.0000 4 44.0000 5 The vector x 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 b = A'*x: 1 11.0000 2 56.0000 3 158.000 4 340.000 R8UT_MV_TEST R8UT_MV computes A*x for an R8UT matrix. The matrix A: Col: 1 2 3 4 Row --- 1 11.0000 12.0000 13.0000 14.0000 2 22.0000 23.0000 24.0000 3 33.0000 34.0000 4 44.0000 5 The vector x 1 1.00000 2 2.00000 3 3.00000 4 4.00000 b = A*x: 1 130.000 2 209.000 3 235.000 4 176.000 5 0.00000 R8UT_PRINT_TEST R8UT_PRINT prints an R8UT matrix. Matrix rows M = 6 Matrix columns N = 4 The R8UT matrix: Col: 1 2 3 4 Row --- 1 11.0000 12.0000 13.0000 14.0000 2 22.0000 23.0000 24.0000 3 33.0000 34.0000 4 44.0000 5 6 R8UT_PRINT_SOME_TEST R8UT_PRINT_SOME prints some of an R8UT matrix. Matrix rows M = 4 Matrix columns N = 6 Some of the matrix: Col: 4 5 6 Row --- 1 14.0000 15.0000 16.0000 2 24.0000 25.0000 26.0000 3 34.0000 35.0000 36.0000 R8UT_RANDOM_TEST R8UT_RANDOM randomizes an R8UT matrix. Matrix order M, N = 5 4 The matrix: Col: 1 2 3 4 Row --- 1 0.218418 0.956318 0.561695 0.257578 2 0.829509 0.415307 0.109957 3 0.661187E-01 0.438290E-01 4 0.633966 5 R8UT_SL_TEST R8UT_SL solves A*x=b for an upper triangular A. Matrix order N = 5 The upper triangular matrix: Col: 1 2 3 4 5 Row --- 1 1.00000 2.00000 3.00000 4.00000 5.00000 2 2.00000 3.00000 4.00000 5.00000 3 3.00000 4.00000 5.00000 4 4.00000 5.00000 5 5.00000 Solution: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 R8UT_SLT_TEST R8UT_SLT solves A'*x=b for an upper triangular A. Matrix order N = 5 The upper triangular matrix: Col: 1 2 3 4 5 Row --- 1 1.00000 2.00000 3.00000 4.00000 5.00000 2 2.00000 3.00000 4.00000 5.00000 3 3.00000 4.00000 5.00000 4 4.00000 5.00000 5 5.00000 Solution to transposed system: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 R8UT_TO_R8GE_TEST R8UT_TO_R8GE converts an R8UT matrix to R8GE format. The random R8UT matrix: Col: 1 2 3 4 Row --- 1 0.218418 0.956318 0.561695 0.257578 2 0.829509 0.415307 0.109957 3 0.661187E-01 0.438290E-01 4 0.633966 5 The R8GE matrix Col: 1 2 3 4 Row --- 1 0.218418 0.956318 0.561695 0.257578 2 0.00000 0.829509 0.415307 0.109957 3 0.00000 0.00000 0.661187E-01 0.438290E-01 4 0.00000 0.00000 0.00000 0.633966 5 0.00000 0.00000 0.00000 0.00000 R8UT_ZEROS_TEST R8UT_ZEROS zeros out a matrix in R8UT format. Matrix order M, N = 5 4 The matrix: Col: 1 2 3 4 Row --- 1 0.00000 0.00000 0.00000 0.00000 2 0.00000 0.00000 0.00000 3 0.00000 0.00000 4 0.00000 5 R8UT_PRB Normal end of execution. 02 July 2016 9:27:48.350 AM