14 September 2021 10:27:48.959 AM R8UT_TEST 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.958914 0.225164 0.631466 0.332089 2 0.362857 0.568306 0.310183 0.699005 3 0.570575 0.107566 0.661726 0.828833 4 0.340801 0.171543 0.257229 0.423243 5 0.426795 0.986481E-01 0.165669 0.987513 The R8UT matrix Col: 1 2 3 4 Row --- 1 0.958914 0.225164 0.631466 0.332089 2 0.568306 0.310183 0.699005 3 0.661726 0.828833 4 0.423243 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.504189 0.698771 0.528215 0.267680 0.867082 2 0.752646 0.739072 0.361744 0.729205 3 0.910733 0.634640 0.875955 4 0.511469 0.809051 5 0.613078 Determinant is 0.108370 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.476672 0.874522 0.198183 0.331999 0.132330 2 0.979299E-02 0.216208 0.867753 0.545007 3 0.823802 0.561119 0.973480 4 0.934018 0.463430 5 0.180276 The inverse matrix B: Col: 1 2 3 4 5 Row --- 1 2.09788 -187.342 48.6635 144.071 -68.3082 2 102.114 -26.7999 -78.7691 38.4986 3 1.21388 -0.729251 -4.68025 4 1.07064 -2.75227 5 5.54706 The product C = A * B: Col: 1 2 3 4 5 Row --- 1 1.00000 0.00000 0.166533E-14 -0.166533E-15 -0.133227E-14 2 1.00000 0.00000 0.00000 0.00000 3 1.00000 0.00000 0.888178E-15 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.957779 0.771805 0.344211 0.437216 2 0.790199 0.757877 0.856236 3 0.179961E-01 0.169322 4 0.940756 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.596518 0.916628 0.304036 0.738729 2 0.522322 0.144318E-01 0.994174 3 0.497878 0.444549 4 0.371286 5 The R8GE matrix Col: 1 2 3 4 Row --- 1 0.596518 0.916628 0.304036 0.738729 2 0.00000 0.522322 0.144318E-01 0.994174 3 0.00000 0.00000 0.497878 0.444549 4 0.00000 0.00000 0.00000 0.371286 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_TEST Normal end of execution. 14 September 2021 10:27:48.960 AM