Tue May 20 22:27:35 2025 r8ltt_test(): python version: 3.10.12 numpy version: 1.26.4 Test r8ltt(). R8LTT_DET_TEST R8LTT_DET computes the determinant of an R8LTT matrix. Matrix order N = 5 The matrix: Col: 0 1 2 3 4 Row 0 : 0.0489715 1 : 0.8731 0.0489715 2 : 0.294886 0.8731 0.0489715 3 : 0.0476487 0.294886 0.8731 0.0489715 4 : 0.260322 0.0476487 0.294886 0.8731 0.0489715 The determinant = 2.8165413832034327e-07 R8LTT_INDICATOR_TEST R8LTT_INDICATOR sets up an indicator matrix in R8LTT format Matrix order n = 5 The indicator matrix: Col: 0 1 2 3 4 Row 0 : 1 1 : 2 1 2 : 3 2 1 3 : 4 3 2 1 4 : 5 4 3 2 1 R8LTT_INVERSE_TEST R8LTT_INVERSE computes the inverse of an R8LTT matrix. The matrix A: Col: 0 1 2 3 4 Row 0 : 0.96824 1 : 0.359925 0.96824 2 : 0.507353 0.359925 0.96824 3 : 0.101886 0.507353 0.359925 0.96824 4 : 0.278646 0.101886 0.507353 0.359925 0.96824 The inverse matrix B: Col: 0 1 2 3 4 Row 0 : 1.0328 1 : -0.383924 1.0328 2 : -0.398466 -0.383924 1.0328 3 : 0.240617 -0.398466 -0.383924 1.0328 4 : -0.137477 0.240617 -0.398466 -0.383924 1.0328 The product A * B: Col: 0 1 2 3 4 Row 0 : 1 1 : 0 1 2 : 1.11022e-16 0 1 3 : 4.16334e-17 1.11022e-16 0 1 4 : 0 4.16334e-17 1.11022e-16 0 1 R8LTT_MM_TEST R8LTT_MM computes C = A * B for R8LTT matrices. Matrix order N = 5 Factor A: Col: 0 1 2 3 4 Row 0 : 0.0200708 1 : 0.378448 0.0200708 2 : 0.52793 0.378448 0.0200708 3 : 0.922448 0.52793 0.378448 0.0200708 4 : 0.574961 0.922448 0.52793 0.378448 0.0200708 Factor B: Col: 0 1 2 3 4 Row 0 : 0.110324 1 : 0.440197 0.110324 2 : 0.171594 0.440197 0.110324 3 : 0.800279 0.171594 0.440197 0.110324 4 : 0.690771 0.800279 0.171594 0.440197 0.110324 The product C = A * B Col: 0 1 2 3 4 Row 0 : 0.0022143 1 : 0.0505871 0.0022143 2 : 0.228279 0.0505871 0.0022143 3 : 0.415163 0.228279 0.0505871 0.0022143 4 : 0.876809 0.415163 0.228279 0.0505871 0.0022143 The product C = A * B using R8GE format: [[0.0022143 0. 0. 0. 0. ] [0.05058708 0.0022143 0. 0. 0. ] [0.22827929 0.05058708 0.0022143 0. 0. ] [0.41516334 0.22827929 0.05058708 0.0022143 0. ] [0.87680889 0.41516334 0.22827929 0.05058708 0.0022143 ]] R8LTT_MTM_TEST R8LTT_MTM computes C = A' * B for R8LTT matrices. The matrix A: Col: 0 1 2 3 4 Row 0 : 0.257903 1 : 0.183121 0.257903 2 : 0.840928 0.183121 0.257903 3 : 0.650711 0.840928 0.183121 0.257903 4 : 0.320575 0.650711 0.840928 0.183121 0.257903 The matrix B: Col: 0 1 2 3 4 Row 0 : 0.569993 1 : 0.0569623 0.569993 2 : 0.348068 0.0569623 0.569993 3 : 0.427402 0.348068 0.0569623 0.569993 4 : 0.707478 0.427402 0.348068 0.0569623 0.569993 The product C = A' * B, using R8LTT format: [[0.95504949 0.51578518 0.62797074 0.38916143 0.18272547] [0.89820783 0.72824955 0.37877068 0.51638871 0.37090072] [0.76297222 0.43784363 0.45013409 0.15227868 0.47932267] [0.23978224 0.16803405 0.07842927 0.1574338 0.10437747] [0.18246077 0.1102283 0.08976785 0.01469075 0.14700281]] The product C = A' * B, using R8GE format: 0.9550494916232706 R8LTT_MTV_TEST R8LTT_MTV computes a matrix product b=A'*x for an R8LTT matrix. Matrix order N = 5 The matrix A: Col: 0 1 2 3 4 Row 0 : 1 1 : 2 1 2 : 3 2 1 3 : 4 3 2 1 4 : 5 4 3 2 1 The vector X: 0: 1 1: 2 2: 3 3: 4 4: 5 The vector b=A'*x: 0: 55 1: 40 2: 26 3: 14 4: 5 R8LTT_MV_TEST R8LTT_MV computes a product b=A*x for an R8LTT matrix. Matrix order N = 5 The R8LTT matrix A: Col: 0 1 2 3 4 Row 0 : 1 1 : 2 1 2 : 3 2 1 3 : 4 3 2 1 4 : 5 4 3 2 1 Vector x: 0: 1 1: 2 2: 3 3: 4 4: 5 Vector b = A*x: 0: 1 1: 4 2: 10 3: 20 4: 35 R8LTT_PRINT_TEST R8LTT_PRINT prints an R8LTT matrix. Matrix order N = 5 The matrix: Col: 0 1 2 3 4 Row 0 : 1 1 : 2 1 2 : 3 2 1 3 : 4 3 2 1 4 : 5 4 3 2 1 R8LTT_PRINT_SOME_TEST R8LTT_PRINT_SOME prints some of an R8LTT matrix. Matrix order N = 6 Some of the matrix: Col: 0 1 2 3 Row 1 : 2 1 2 : 3 2 1 3 : 4 3 2 1 4 : 5 4 3 2 R8LTT_RANDOM_TEST R8LTT_RANDOM randomizes an R8LTT matrix. Matrix order N = 5 Matrix A: Col: 0 1 2 3 4 Row 0 : 0.448308 1 : 0.0922862 0.448308 2 : 0.394121 0.0922862 0.448308 3 : 0.290753 0.394121 0.0922862 0.448308 4 : 0.664363 0.290753 0.394121 0.0922862 0.448308 R8LTT_SL_TEST R8LTT_SL solves a linear system A*x=b with R8LTT matrix Matrix order N = 5 Matrix A: Col: 0 1 2 3 4 Row 0 : 0.648482 1 : 0.358322 0.648482 2 : 0.0300227 0.358322 0.648482 3 : 0.100484 0.0300227 0.358322 0.648482 4 : 0.749946 0.100484 0.0300227 0.358322 0.648482 Right hand side b: 0: 0.648482 1: 1.65529 2: 2.69211 3: 3.82942 4: 5.71668 Solution x: 0: 1 1: 2 2: 3 3: 4 4: 5 R8LTT_SLT_TEST R8LTT_SLT solves a linear system A'x=b with R8LTT matrix Matrix order N = 5 Matrix A: Col: 0 1 2 3 4 Row 0 : 0.491456 1 : 0.758691 0.491456 2 : 0.436226 0.758691 0.491456 3 : 0.0317721 0.436226 0.758691 0.491456 4 : 0.309341 0.0317721 0.436226 0.758691 0.491456 Right hand side b: 0: 4.99131 1: 5.16275 2: 6.69026 3: 5.75928 4: 2.45728 Solution x: 0: 1 1: 2 2: 3 3: 4 4: 5 R8LTT_TO_R8GE_TEST R8LTT_TO_R8GE converts an R8LTT matrix to R8GE format. The random R8LTT matrix: Col: 0 1 2 3 4 Row 0 : 0.43948 1 : 0.0394907 0.43948 2 : 0.270505 0.0394907 0.43948 3 : 0.0547476 0.270505 0.0394907 0.43948 4 : 0.565329 0.0547476 0.270505 0.0394907 0.43948 The R8GE matrix: [[0.43948019 0. 0. 0. 0. ] [0.03949067 0.43948019 0. 0. 0. ] [0.27050522 0.03949067 0.43948019 0. 0. ] [0.05474757 0.27050522 0.03949067 0.43948019 0. ] [0.5653291 0.05474757 0.27050522 0.03949067 0.43948019]] R8LTT_ZEROS_TEST R8LTT_ZEROS zeros out space for an R8LTT matrix. Matrix order N = 5 The matrix: Col: 0 1 2 3 4 Row 0 : 0 1 : 0 0 2 : 0 0 0 3 : 0 0 0 0 4 : 0 0 0 0 0 r8ltt_test(): Normal end of execution. Tue May 20 22:27:35 2025