08-Jan-2022 09:09:11 r8lt_test() MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2 Test r8lt(). 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.814724 0.0975404 0.157613 0.141886 2 0.905792 0.278498 0.970593 0.421761 3 0.126987 0.546882 0.957167 0.915736 4 0.913376 0.957507 0.485376 0.792207 5 0.632359 0.964889 0.80028 0.959492 The R8LT matrix: Col: 1 2 3 4 Row --- 1 0.814724 2 0.905792 0.278498 3 0.126987 0.546882 0.957167 4 0.913376 0.957507 0.485376 0.792207 5 0.632359 0.964889 0.80028 0.959492 R8LT_DET_TEST(): R8LT_DET computes the determinant of an R8LT matrix The matrix A: Col: 1 2 3 4 5 Row --- 1 0.655741 2 0.0357117 0.75774 3 0.849129 0.743132 0.171187 4 0.933993 0.392227 0.706046 0.276923 5 0.678735 0.655478 0.0318328 0.0461714 0.0971318 Determinant is 0.00228793 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.823458 2 0.694829 0.438744 3 0.317099 0.381558 0.186873 4 0.950222 0.765517 0.489764 0.646313 5 0.0344461 0.7952 0.445586 0.709365 0.754687 The inverse matrix B: Col: 1 2 3 4 5 Row --- 1 1.21439 2 -1.9232 2.27923 3 1.86614 -4.65376 5.35124 4 -0.921641 0.82693 -4.05507 1.54724 5 1.73549 -0.431157 0.652043 -1.45432 1.32505 The product C = A * B: Col: 1 2 3 4 5 Row --- 1 1 2 0 1 3 0 0 1 4 0 0 0 1 5 -1.11022e-16 1.11022e-16 0 0 1 R8LT_MM_TEST R8LT_MM computes C = A * B for R8LT matrices. The matrix A: 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 A: 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'*x for an R8LT matrix. The matrix A: 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 The vector x 1: 1 2: 2 3: 3 4: 4 5: 5 b = A'*x: 1: 565 2: 568 3: 536 4: 446 R8LT_MV_TEST R8LT_MV computes A*x for an R8LT matrix. The matrix A: 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 The vector x 1: 1 2: 2 3: 3 4: 4 b = A*x: 1: 11 2: 65 3: 194 4: 430 5: 530 R8LT_PRINT_TEST(): R8LT_PRINT prints an R8LT matrix; The R8LT 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_PRINT_SOME_TEST R8LT_PRINT_SOME prints some of an R8LT matrix; Rows 2-3, Cols 2-4: Col: 2 3 4 Row --- 2 22 3 32 33 R8LT_RANDOM_TEST R8LT_RANDOM randomizes an R8LT matrix; The random R8LT matrix: Col: 1 2 3 4 5 Row --- 1 0.276025 2 0.679703 0.959744 3 0.655098 0.340386 0.255095 4 0.162612 0.585268 0.505957 0.959291 5 0.118998 0.223812 0.699077 0.547216 0.149294 6 0.498364 0.751267 0.890903 0.138624 0.257508 R8LT_SL_TEST(): R8LT_SL solves A*x=b for an R8LT matrix A. The R8LT matrix: Col: 1 2 3 4 5 Row --- 1 1 2 1 2 3 1 2 3 4 1 2 3 4 5 1 2 3 4 5 Solution: 1: 1 2: 2 3: 3 4: 4 5: 5 R8LT_SLT_TEST R8LT_SLT solves A'*x=b for an R8LT matrix A The R8LT matrix: Col: 1 2 3 4 5 Row --- 1 1 2 1 2 3 1 2 3 4 1 2 3 4 5 1 2 3 4 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. The random R8LT matrix: Col: 1 2 3 4 Row --- 1 0.840717 2 0.254282 0.349984 3 0.814285 0.196595 0.473289 4 0.243525 0.251084 0.35166 0.585264 5 0.929264 0.616045 0.830829 0.549724 The R8GE matrix Col: 1 2 3 4 Row --- 1 0.840717 0 0 0 2 0.254282 0.349984 0 0 3 0.814285 0.196595 0.473289 0 4 0.243525 0.251084 0.35166 0.585264 5 0.929264 0.616045 0.830829 0.549724 R8LT_ZEROS_TEST(): R8LT_ZEROS zeros out space for an R8LT matrix. Matrix order M, N = 5, 4 Matrix A: Col: 1 2 3 4 Row --- 1 0 2 0 0 3 0 0 0 4 0 0 0 0 5 0 0 0 0 r8lt_test(): Normal end of execution. 08-Jan-2022 09:09:11