14 September 2021 10:27:48.128 AM R8LT_TEST 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.443689 0.875263 0.768891 0.929323 2 0.614670 0.706941 0.188690 0.387637 3 0.281527 0.599347 0.669767 0.954172 4 0.130758 0.780843E-01 0.238984 0.649772 5 0.363581 0.316556 0.667498 0.779434 The R8LT matrix Col: 1 2 3 4 Row --- 1 0.443689 2 0.614670 0.706941 3 0.281527 0.599347 0.669767 4 0.130758 0.780843E-01 0.238984 0.649772 5 0.363581 0.316556 0.667498 0.779434 R8LT_DET_TEST R8LT_DET computes the determinant of an R8LT matrix. The matrix A: Col: 1 2 3 4 5 Row --- 1 0.895532 2 0.849938E-01 0.641384 3 0.535464 0.193334 0.695354 4 0.488869 0.201395E-01 0.204782 0.157723 5 0.146980 0.249765 0.365739 0.311273 0.787002 Determinant is 0.495766E-01 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.122215 2 0.151315 0.434231E-01 3 0.835476 0.676050 0.942171 4 0.451533 0.709349E-01 0.538759 0.549781 5 0.922331 0.760402 0.239080 0.595140 0.317470 The inverse matrix B: Col: 1 2 3 4 5 Row --- 1 8.18227 2 -28.5124 23.0292 3 13.2033 -16.5245 1.06138 4 -15.9798 13.2219 -1.04010 1.81890 5 64.5345 -67.5013 1.15050 -3.40978 3.14990 The product C = A * B: Col: 1 2 3 4 5 Row --- 1 1.00000 2 0.00000 1.00000 3 0.00000 0.177636E-14 1.00000 4 0.00000 0.00000 0.00000 1.00000 5 0.00000 0.355271E-14 0.00000 0.00000 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.519851 2 0.203898 0.284569 3 0.321391 0.638963 0.479853 4 0.708206 0.879292 0.989969 0.567691 5 0.167339 0.205749 0.454090 0.360189 0.992338 6 0.935759 0.384611 0.395892 0.527422 0.497264 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.625768 2 0.445358 0.916776 3 0.957389 0.405093 0.974167 4 0.534812 0.929233 0.450285 0.628260 5 0.828340 0.426793 0.304900 0.761532 The R8GE matrix Col: 1 2 3 4 Row --- 1 0.625768 0.00000 0.00000 0.00000 2 0.445358 0.916776 0.00000 0.00000 3 0.957389 0.405093 0.974167 0.00000 4 0.534812 0.929233 0.450285 0.628260 5 0.828340 0.426793 0.304900 0.761532 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_TEST Normal end of execution. 14 September 2021 10:27:48.129 AM