Wed Oct 8 08:50:31 2025 r8lt_test(): python version: 3.10.12 numpy version: 1.26.4 Test r8lt(). i4_log_10_test(): i4_log_10(): whole part of log base 10, X, i4_log_10 0 0 1 0 2 0 3 0 9 0 10 1 11 1 99 1 101 2 -1 0 -2 0 -3 0 -9 0 r8ge_print_test(): r8ge_print prints an R8GE matrix. Here is an R8GE: Col: 0 1 2 3 4 Row 0 : 11 12 13 14 15 1 : 21 22 23 24 25 2 : 31 32 33 34 35 3 : 41 42 43 44 45 Col: 5 Row 0 : 16 1 : 26 2 : 36 3 : 46 r8ge_print_some_test(): r8ge_print_some prints some of an R8GE matrix. Rows 0:2, Cols 3:5: Col: 3 4 5 Row 0 : 14 15 16 1 : 24 25 26 2 : 34 35 36 r8ge_to_r8lt_test(): r8ge_to_r8lt() converts an R8GE matrix to R8LT format. The random R8GE matrix: Col: 0 1 2 3 Row 0 : 0.829292 0.140055 0.480683 0.824443 1 : 0.532189 0.0609981 0.0180261 0.963423 2 : 0.27128 0.779971 0.210222 0.880824 3 : 0.480097 0.0463718 0.406392 0.0760207 4 : 0.972526 0.640908 0.00539382 0.38357 The R8LT matrix: Col: 0 1 2 3 Row 0 : 0.829292 0 0 0 1 : 0.532189 0.0609981 0 0 2 : 0.27128 0.779971 0.210222 0 3 : 0.480097 0.0463718 0.406392 0.0760207 4 : 0.972526 0.640908 0.00539382 0.38357 r8lt_det_test(): r8lt_det() computes the determinant of an R8LT matrix. The matrix A: Col: 0 1 2 3 4 Row 0 : 0.27177 0 0 0 0 1 : 0.510211 0.061003 0 0 0 2 : 0.456983 0.303918 0.293242 0 0 3 : 0.879491 0.0195996 0.0768896 0.708435 0 4 : 0.957845 0.206257 0.597339 0.0703215 0.227847 Determinant is 0.000784735 r8lt_indicator_test(): r8lt_indicator sets up an indicator matrix in R8LT format Matrix rows M = 5 Matrix columns N = 4 The indicator matrix: Col: 0 1 2 3 Row 0 : 11 0 0 0 1 : 21 22 0 0 2 : 31 32 33 0 3 : 41 42 43 44 4 : 51 52 53 54 r8lt_inverse_test(): r8lt_inverse() computes the inverse of an R8LT matrix. The matrix A: Col: 0 1 2 3 4 Row 0 : 0.0217959 0 0 0 0 1 : 0.914999 0.0485959 0 0 0 2 : 0.899466 0.980065 0.134991 0 0 3 : 0.393227 0.000197011 0.701979 0.243752 0 4 : 0.211846 0.145048 0.763186 0.333724 0.951183 The inverse matrix B: Col: 0 1 2 3 4 Row 0 : 45.8802 0 0 0 0 1 : -863.864 20.5778 0 0 0 2 : 5966.14 -149.4 7.4079 0 0 3 : -17255.1 430.238 -21.334 4.10253 0 4 : 1388.55 -34.2162 1.54129 -1.43938 1.05132 The product A * B: Col: 0 1 2 3 4 Row 0 : 1 0 0 0 0 1 : 0 1 0 0 0 2 : 0 0 1 0 0 3 : 0 0 0 1 0 4 : 0 0 0 0 1 r8lt_mm_test(): r8lt_mm computes C = A * B for R8LT matrices. The matrix A: Col: 0 1 2 3 4 Row 0 : 1 0 0 0 0 1 : 1 1 0 0 0 2 : 1 1 1 0 0 3 : 1 1 1 1 0 4 : 1 1 1 1 1 The product C = A * A Col: 0 1 2 3 4 Row 0 : 1 0 0 0 0 1 : 2 1 0 0 0 2 : 3 2 1 0 0 3 : 4 3 2 1 0 4 : 5 4 3 2 1 r8lt_mtm_test(): r8lt_mtm computes C = A' * B for R8LT matrices. The matrix A: Col: 0 1 2 3 4 Row 0 : 1 0 0 0 0 1 : 1 1 0 0 0 2 : 1 1 1 0 0 3 : 1 1 1 1 0 4 : 1 1 1 1 1 The product C = A' * A Col: 0 1 2 3 4 Row 0 : 5 4 3 2 1 1 : 4 4 3 2 1 2 : 3 3 3 2 1 3 : 2 2 2 2 1 4 : 1 1 1 1 1 r8lt_mtv_test(): r8lt_mtv computes A'*x for an R8LT matrix. The matrix A: Col: 0 1 2 3 Row 0 : 11 0 0 0 1 : 21 22 0 0 2 : 31 32 33 0 3 : 41 42 43 44 4 : 51 52 53 54 The vector x 0 1 1 2 2 3 3 4 4 5 b = A'*x: 0 565 1 568 2 536 3 446 r8lt_mv_test(): r8lt_mv computes A*x for an R8LT matrix. The matrix A: Col: 0 1 2 3 Row 0 : 11 0 0 0 1 : 21 22 0 0 2 : 31 32 33 0 3 : 41 42 43 44 4 : 51 52 53 54 The vector x 0 1 1 2 2 3 3 4 b = A*x: 0 11 1 65 2 194 3 430 4 530 r8lt_print_test(): r8lt_print prints an R8LT matrix. The R8LT matrix: Col: 0 1 2 3 4 Row 0 : 11 0 0 0 0 1 : 21 22 0 0 0 2 : 31 32 33 0 0 3 : 41 42 43 44 0 Col: 5 Row r8lt_print_some_test(): r8lt_print_some prints some of an R8LT matrix. Rows 1-3, Columns 1-2: Col: 1 2 Row 1 : 22 0 2 : 32 33 3 : 42 43 r8lt_random_test(): r8lt_random() randomizes an R8LT matrix. Matrix order M, N = 5, 4 Matrix A: Col: 0 1 2 3 Row 0 : 0.0306368 0 0 0 1 : 0.0033415 0.336926 0 0 2 : 0.206675 0.575376 0.821141 0 3 : 0.0705248 0.845985 0.552644 0.95173 4 : 0.700418 0.931105 0.443667 0.849726 r8lt_sl_test(): r8lt_sl() solves A*x=b for an R8LT matrix A. The R8LT matrix: Col: 0 1 2 3 4 Row 0 : 1 0 0 0 0 1 : 1 2 0 0 0 2 : 1 2 3 0 0 3 : 1 2 3 4 0 4 : 1 2 3 4 5 Solution: 0 1 1 2 2 3 3 4 4 5 r8lt_slt_test(): r8lt_slt() solves A'*x=b for an R8LT matrix A The R8LT matrix: Col: 0 1 2 3 4 Row 0 : 1 0 0 0 0 1 : 1 2 0 0 0 2 : 1 2 3 0 0 3 : 1 2 3 4 0 4 : 1 2 3 4 5 Solution to transposed system: 0 1 1 2 2 3 3 4 4 5 r8lt_to_r8ge_test(): r8lt_to_r8ge() converts an R8LT matrix to R8GE format. The R8LT matrix: Col: 0 1 2 3 4 Row 0 : 11 0 0 0 0 1 : 21 22 0 0 0 2 : 31 32 33 0 0 3 : 41 42 43 44 0 Col: 5 Row The R8GE matrix: Col: 0 1 2 3 4 Row 0 : 11 0 0 0 0 1 : 21 22 0 0 0 2 : 31 32 33 0 0 3 : 41 42 43 44 0 Col: 5 Row 0 : 0 1 : 0 2 : 0 3 : 0 r8lt_zeros_test(): r8lt_zeros zeros out space for an R8LT matrix. Matrix A: Col: 0 1 2 3 Row 0 : 0 0 0 0 1 : 0 0 0 0 2 : 0 0 0 0 3 : 0 0 0 0 4 : 0 0 0 0 r8vec_indicator1_test(): r8vec_indicator1 returns the 1-based indicator matrix. The 1-based indicator vector: 0 1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 10 r8lt_test(): Normal end of execution. Wed Oct 8 08:50:31 2025