Tue Aug 23 09:17:29 2022 r8ltt_test(): Python version: 3.6.9 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.688033 1 : 0.479553 0.688033 2 : 0.326593 0.479553 0.688033 3 : 0.619433 0.326593 0.479553 0.688033 4 : 0.617516 0.619433 0.326593 0.479553 0.688033 The determinant = 0.15418637659574033 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.983506 1 : 0.0371989 0.983506 2 : 0.777408 0.0371989 0.983506 3 : 0.0466105 0.777408 0.0371989 0.983506 4 : 0.276265 0.0466105 0.777408 0.0371989 0.983506 The inverse matrix B: Col: 0 1 2 3 4 Row 0 : 1.01677 1 : -0.0384571 1.01677 2 : -0.802248 -0.0384571 1.01677 3 : 0.0125545 -0.802248 -0.0384571 1.01677 4 : 0.349872 0.0125545 -0.802248 -0.0384571 1.01677 The product A * B: Col: 0 1 2 3 4 Row 0 : 1 1 : 0 1 2 : 0 0 1 3 : 0 0 0 1 4 : 0 0 0 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.690123 1 : 0.0208953 0.690123 2 : 0.559807 0.0208953 0.690123 3 : 0.732591 0.559807 0.0208953 0.690123 4 : 0.506908 0.732591 0.559807 0.0208953 0.690123 Factor B: Col: 0 1 2 3 4 Row 0 : 0.724376 1 : 0.89946 0.724376 2 : 0.896363 0.89946 0.724376 3 : 0.468542 0.896363 0.89946 0.724376 4 : 0.165603 0.468542 0.896363 0.89946 0.724376 The product C = A * B Col: 0 1 2 3 4 Row 0 : 0.499908 1 : 0.635874 0.499908 2 : 1.04291 0.635874 0.499908 3 : 1.37628 1.04291 0.635874 0.499908 4 : 1.652 1.37628 1.04291 0.635874 0.499908 The product C = A * B using R8GE format: [[0.49990813 0. 0. 0. 0. ] [0.63587381 0.49990813 0. 0. 0. ] [1.04290548 0.63587381 0.49990813 0. 0. ] [1.37627651 1.04290548 0.63587381 0.49990813 0. ] [1.65199514 1.37627651 1.04290548 0.63587381 0.49990813]] R8LTT_MTM_TEST R8LTT_MTM computes C = A' * B for R8LTT matrices. The matrix A: Col: 0 1 2 3 4 Row 0 : 0.80597 1 : 0.306875 0.80597 2 : 0.950485 0.306875 0.80597 3 : 0.0401585 0.950485 0.306875 0.80597 4 : 0.535564 0.0401585 0.950485 0.306875 0.80597 The matrix B: Col: 0 1 2 3 4 Row 0 : 0.328024 1 : 0.132589 0.328024 2 : 0.365165 0.132589 0.328024 3 : 0.0166234 0.365165 0.132589 0.328024 4 : 0.22358 0.0166234 0.365165 0.132589 0.328024 The product C = A' * B, using R8LTT format: [[0.77255895 0.25025383 0.51267584 0.08418294 0.17567788] [0.24370192 0.65281726 0.24135093 0.31710647 0.01317296] [0.51192322 0.23472326 0.65214969 0.22668644 0.31178188] [0.08200928 0.29941331 0.21892297 0.30506565 0.10066246] [0.18019905 0.01339795 0.294312 0.10686279 0.26437732]] The product C = A' * B, using R8GE format: 0.7725589502697685 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.349782 1 : 0.699161 0.349782 2 : 0.422579 0.699161 0.349782 3 : 0.197564 0.422579 0.699161 0.349782 4 : 0.643717 0.197564 0.422579 0.699161 0.349782 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.723019 1 : 0.676601 0.723019 2 : 0.306733 0.676601 0.723019 3 : 0.810299 0.306733 0.676601 0.723019 4 : 0.978001 0.810299 0.306733 0.676601 0.723019 Right hand side b: 0: 0.723019 1: 2.12264 2: 3.82899 3: 6.34565 4: 9.8403 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.747637 1 : 0.308544 0.747637 2 : 0.0254415 0.308544 0.747637 3 : 0.606748 0.0254415 0.308544 0.747637 4 : 0.231197 0.606748 0.0254415 0.308544 0.747637 Right hand side b: 0: 5.02402 1: 5.55641 2: 3.6043 3: 4.53327 4: 3.73819 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.941612 1 : 0.679327 0.941612 2 : 0.8074 0.679327 0.941612 3 : 0.857315 0.8074 0.679327 0.941612 4 : 0.912674 0.857315 0.8074 0.679327 0.941612 The R8GE matrix: [[0.94161201 0. 0. 0. 0. ] [0.67932669 0.94161201 0. 0. 0. ] [0.80740038 0.67932669 0.94161201 0. 0. ] [0.85731502 0.80740038 0.67932669 0.94161201 0. ] [0.91267372 0.85731502 0.80740038 0.67932669 0.94161201]] 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 Aug 23 09:17:29 2022