Tue May 20 22:27:35 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.773574 0.470371 0.76458 0.746315 1 : 0.58081 0.994446 0.890119 0.36487 2 : 0.52526 0.0798752 0.615918 0.933733 3 : 0.960011 0.299993 0.887656 0.310668 4 : 0.23615 0.668688 0.627229 0.824524 The R8LT matrix: Col: 0 1 2 3 Row 0 : 0.773574 0 0 0 1 : 0.58081 0.994446 0 0 2 : 0.52526 0.0798752 0.615918 0 3 : 0.960011 0.299993 0.887656 0.310668 4 : 0.23615 0.668688 0.627229 0.824524 r8lt_det_test(): r8lt_det() computes the determinant of an R8LT matrix. The matrix A: Col: 0 1 2 3 4 Row 0 : 0.640367 0 0 0 0 1 : 0.667924 0.404331 0 0 0 2 : 0.431634 0.122797 0.374373 0 0 3 : 0.257127 0.223728 0.45208 0.529866 0 4 : 0.213952 0.850326 0.450022 0.771964 0.133002 Determinant is 0.00683116 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.913011 0 0 0 0 1 : 0.654678 0.59116 0 0 0 2 : 0.520858 0.014354 0.447093 0 0 3 : 0.141622 0.706531 0.107706 0.911799 0 4 : 0.504844 0.898382 0.296133 0.201846 0.505062 The inverse matrix B: Col: 0 1 2 3 4 Row 0 : 1.09528 0 0 0 0 1 : -1.21296 1.69159 0 0 0 2 : -1.23704 -0.0543086 2.23667 0 0 3 : 0.915898 -1.30436 -0.264206 1.09673 0 4 : 1.42204 -2.4558 -1.20584 -0.438305 1.97996 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 2.77556e-17 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.404828 0 0 0 1 : 0.795999 0.812205 0 0 2 : 0.438351 0.00957068 0.0955577 0 3 : 0.162719 0.914107 0.227243 0.39891 4 : 0.493003 0.0113282 0.388166 0.0241406 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. Tue May 20 22:27:35 2025