6 October 2025 6:33:35.716 PM r8ltt_test(): Fortran90 version Test r8ltt(). R8LTT_DET_TEST R8LTT_DET computes the determinant of an R8LTT matrix. Matrix order N = 5 The matrix: Col: 1 2 3 4 5 Row --- 1 0.994214E-01 2 0.452507 0.994214E-01 3 0.714086 0.452507 0.994214E-01 4 0.454479 0.714086 0.452507 0.994214E-01 5 0.711186 0.454479 0.714086 0.452507 0.994214E-01 The determinant = 0.971402E-05 R8LTT_INDICATOR_TEST R8LTT_INDICATOR sets up an indicator matrix in R8LTT format Matrix order N = 5 The indicator matrix: Col: 1 2 3 4 5 Row --- 1 1.00000 2 2.00000 1.00000 3 3.00000 2.00000 1.00000 4 4.00000 3.00000 2.00000 1.00000 5 5.00000 4.00000 3.00000 2.00000 1.00000 R8LTT_INVERSE_TEST R8LTT_INVERSE computes the inverse of an R8LTT matrix. The matrix A: Col: 1 2 3 4 5 Row --- 1 0.550544 2 0.275171E-01 0.550544 3 0.443652 0.275171E-01 0.550544 4 0.179011 0.443652 0.275171E-01 0.550544 5 0.778561 0.179011 0.443652 0.275171E-01 0.550544 The inverse matrix B: Col: 1 2 3 4 5 Row --- 1 1.81638 2 -0.907859E-01 1.81638 3 -1.45918 -0.907859E-01 1.81638 4 -0.444510 -1.45918 -0.907859E-01 1.81638 5 -1.34106 -0.444510 -1.45918 -0.907859E-01 1.81638 The product A * B: Col: 1 2 3 4 5 Row --- 1 1.00000 2 0.00000 1.00000 3 0.111022E-15 0.00000 1.00000 4 0.00000 0.111022E-15 0.00000 1.00000 5 0.00000 0.00000 0.111022E-15 0.00000 1.00000 R8LTT_MM_TEST R8LTT_MM computes C = A * B for R8LTT matrices. Matrix order N = 5 Factor A: Col: 1 2 3 4 5 Row --- 1 0.937025 2 0.123253 0.937025 3 0.519022 0.123253 0.937025 4 0.489260 0.519022 0.123253 0.937025 5 0.490990 0.489260 0.519022 0.123253 0.937025 Factor B: Col: 1 2 3 4 5 Row --- 1 0.488331 2 0.641604 0.488331 3 0.208078 0.641604 0.488331 4 0.235795 0.208078 0.641604 0.488331 5 0.369308 0.235795 0.208078 0.641604 0.488331 The product C = A * B Col: 1 2 3 4 5 Row --- 1 0.457579 2 0.661387 0.457579 3 0.527508 0.661387 0.457579 4 0.818519 0.527508 0.661387 0.457579 5 1.03679 0.818519 0.527508 0.661387 0.457579 The R8GE product C: Col: 1 2 3 4 5 Row --- 1 0.457579 0.00000 0.00000 0.00000 0.00000 2 0.661387 0.457579 0.00000 0.00000 0.00000 3 0.527508 0.661387 0.457579 0.00000 0.00000 4 0.818519 0.527508 0.661387 0.457579 0.00000 5 1.03679 0.818519 0.527508 0.661387 0.457579 R8LTT_MTM_TEST R8LTT_MTM computes C = A' * B for R8LTT matrices. The matrix A: Col: 1 2 3 4 5 Row --- 1 0.435507 2 0.579651 0.435507 3 0.168956 0.579651 0.435507 4 0.470515 0.168956 0.579651 0.435507 5 0.320246 0.470515 0.168956 0.579651 0.435507 The matrix B: Col: 1 2 3 4 5 Row --- 1 0.471262 2 0.318567 0.471262 3 0.254876 0.318567 0.471262 4 0.520977 0.254876 0.318567 0.471262 5 0.667191 0.520977 0.254876 0.318567 0.471262 The product C = A' * B: Col: 1 2 3 4 5 Row --- 1 0.891751 0.613755 0.311136 0.323756 0.150920 2 0.688423 0.678085 0.446914 0.229513 0.221736 3 0.525711 0.374499 0.432958 0.326991 0.796226E-01 4 0.613626 0.412984 0.286477 0.389895 0.273167 5 0.290566 0.226889 0.111000 0.138738 0.205238 The R8GE product C = A' * B: Col: 1 2 3 4 5 Row --- 1 0.891751 0.613755 0.311136 0.323756 0.150920 2 0.688423 0.678085 0.446914 0.229513 0.221736 3 0.525711 0.374499 0.432958 0.326991 0.796226E-01 4 0.613626 0.412984 0.286477 0.389895 0.273167 5 0.290566 0.226889 0.111000 0.138738 0.205238 R8LTT_MTV_TEST R8LTT_MTV computes a matrix product b=A'*x for an R8LTT matrix. Matrix order N = 5 The matrix A: Col: 1 2 3 4 5 Row --- 1 1.00000 2 2.00000 1.00000 3 3.00000 2.00000 1.00000 4 4.00000 3.00000 2.00000 1.00000 5 5.00000 4.00000 3.00000 2.00000 1.00000 The vector X: 1: 1.0000000 2: 2.0000000 3: 3.0000000 4: 4.0000000 5: 5.0000000 The vector b=A'*x: 1: 55.000000 2: 40.000000 3: 26.000000 4: 14.000000 5: 5.0000000 R8LTT_MV_TEST R8LTT_MV computes a product b=A*x for an R8LTT matrix. Matrix order N = 5 The R8LTT matrix A: Col: 1 2 3 4 5 Row --- 1 1.00000 2 2.00000 1.00000 3 3.00000 2.00000 1.00000 4 4.00000 3.00000 2.00000 1.00000 5 5.00000 4.00000 3.00000 2.00000 1.00000 Vector x: 1: 1.0000000 2: 2.0000000 3: 3.0000000 4: 4.0000000 5: 5.0000000 Vector b = A*x: 1: 1.0000000 2: 4.0000000 3: 10.000000 4: 20.000000 5: 35.000000 R8LTT_PRINT_TEST R8LTT_PRINT prints an R8LTT matrix. Matrix order N = 5 The matrix: Col: 1 2 3 4 5 Row --- 1 1.00000 2 2.00000 1.00000 3 3.00000 2.00000 1.00000 4 4.00000 3.00000 2.00000 1.00000 5 5.00000 4.00000 3.00000 2.00000 1.00000 R8LTT_PRINT_SOME_TEST R8LTT_PRINT_SOME prints some of an R8LTT matrix. Matrix order N = 6 Some of the matrix: Col: 1 2 3 4 Row --- 2 2.00000 1.00000 3 3.00000 2.00000 1.00000 4 4.00000 3.00000 2.00000 1.00000 5 5.00000 4.00000 3.00000 2.00000 R8LTT_RANDOM_TEST R8LTT_RANDOM randomizes an R8LTT matrix. Matrix order N = 5 Matrix A: Col: 1 2 3 4 5 Row --- 1 0.230514 2 0.103244E-01 0.230514 3 0.423611 0.103244E-01 0.230514 4 0.830360 0.423611 0.103244E-01 0.230514 5 0.884197 0.830360 0.423611 0.103244E-01 0.230514 R8LTT_SL_TEST R8LTT_SL solves a linear system A*x=b with R8LTT matrix Matrix order N = 5 Matrix A: Col: 1 2 3 4 5 Row --- 1 0.855505 2 0.403880 0.855505 3 0.694087 0.403880 0.855505 4 0.630873 0.694087 0.403880 0.855505 5 0.988934 0.630873 0.694087 0.403880 0.855505 Right hand side b: 1: 0.85550455 2: 2.1148894 3: 4.0683610 4: 6.6527052 5: 10.225983 Solution x: 1: 1.0000000 2: 2.0000000 3: 3.0000000 4: 4.0000000 5: 5.0000000 R8LTT_SLT_TEST R8LTT_SLT solves a linear system A'x=b with R8LTT matrix Matrix order N = 5 Matrix A: Col: 1 2 3 4 5 Row --- 1 0.991536 2 0.937027 0.991536 3 0.381563 0.937027 0.991536 4 0.958254 0.381563 0.937027 0.991536 5 0.178934 0.958254 0.381563 0.937027 0.991536 Right hand side b: 1: 8.7379631 2: 11.111675 3: 8.6305336 4: 8.6512797 5: 4.9576817 Solution x: 1: 1.0000000 2: 2.0000000 3: 3.0000000 4: 4.0000000 5: 5.0000000 R8LTT_TO_R8GE_TEST R8LTT_TO_R8GE converts an R8LTT matrix to R8GE format. The random R8LTT matrix: Col: 1 2 3 4 5 Row --- 1 0.204845 2 0.910343 0.204845 3 0.330728 0.910343 0.204845 4 0.422534E-02 0.330728 0.910343 0.204845 5 0.744027 0.422534E-02 0.330728 0.910343 0.204845 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 0.204845 0.00000 0.00000 0.00000 0.00000 2 0.910343 0.204845 0.00000 0.00000 0.00000 3 0.330728 0.910343 0.204845 0.00000 0.00000 4 0.422534E-02 0.330728 0.910343 0.204845 0.00000 5 0.744027 0.422534E-02 0.330728 0.910343 0.204845 R8LTT_ZEROS_TEST R8LTT_ZEROS zeros out space for an R8LTT matrix. Matrix order N = 5 The matrix: Col: 1 2 3 4 5 Row --- 1 0.00000 2 0.00000 0.00000 3 0.00000 0.00000 0.00000 4 0.00000 0.00000 0.00000 0.00000 5 0.00000 0.00000 0.00000 0.00000 0.00000 r8ltt_test(): Normal end of execution. 6 October 2025 6:33:35.717 PM