Wed Oct 8 08:50:36 2025 r8utt_test(): python version: 3.10.12 numpy version: 1.26.4 Test r8utt(). r8utt_det_test(): r8utt_det() computes the determinant of an R8UTT matrix. Matrix order N = 5 The matrix: Col: 0 1 2 3 4 Row 0 : 0.0886058 0.263433 0.1343 0.746896 0.469622 1 : 0.0886058 0.263433 0.1343 0.746896 2 : 0.0886058 0.263433 0.1343 3 : 0.0886058 0.263433 4 : 0.0886058 The determinant = 5.461471842877888e-06 r8utt_indicator_test(): r8utt_indicator() sets up an indicator matrix in R8UTT format Matrix order N = 5 The indicator matrix: Col: 0 1 2 3 4 Row 0 : 1 2 3 4 5 1 : 1 2 3 4 2 : 1 2 3 3 : 1 2 4 : 1 R8UTT_INVERSE_TEST R8UTT_INVERSE computes the inverse of an R8UTT matrix. The matrix A: Col: 0 1 2 3 4 Row 0 : 0.778041 0.748485 0.762939 0.916877 0.189899 1 : 0.778041 0.748485 0.762939 0.916877 2 : 0.778041 0.748485 0.762939 3 : 0.778041 0.748485 4 : 0.778041 The inverse matrix B: Col: 0 1 2 3 4 Row 0 : 1.28528 -1.23646 -0.0708458 -0.234019 1.43799 1 : 1.28528 -1.23646 -0.0708458 -0.234019 2 : 1.28528 -1.23646 -0.0708458 3 : 1.28528 -1.23646 4 : 1.28528 The product A * B: Col: 0 1 2 3 4 Row 0 : 1 0 -2.22045e-16 -2.22045e-16 -6.66134e-16 1 : 1 0 -2.22045e-16 -2.22045e-16 2 : 1 0 -2.22045e-16 3 : 1 0 4 : 1 R8UTT_MM_TEST R8UTT_MM computes C = A * B for R8UTT matrices. Matrix order N = 5 Factor A: Col: 0 1 2 3 4 Row 0 : 0.678069 0.924522 0.634637 0.37276 0.00530727 1 : 0.678069 0.924522 0.634637 0.37276 2 : 0.678069 0.924522 0.634637 3 : 0.678069 0.924522 4 : 0.678069 Factor B: Col: 0 1 2 3 4 Row 0 : 0.810627 0.645034 0.794417 0.440081 0.911655 1 : 0.810627 0.645034 0.794417 0.440081 2 : 0.810627 0.645034 0.794417 3 : 0.810627 0.645034 4 : 0.810627 The product C = A * B Col: 0 1 2 3 4 Row 0 : 0.549661 1.18682 1.64947 1.74439 1.77394 1 : 0.549661 1.18682 1.64947 1.74439 2 : 0.549661 1.18682 1.64947 3 : 0.549661 1.18682 4 : 0.549661 The R8GE product C: [[0.54966135 1.18682042 1.64947155 1.74439269 1.77394096] [0. 0.54966135 1.18682042 1.64947155 1.74439269] [0. 0. 0.54966135 1.18682042 1.64947155] [0. 0. 0. 0.54966135 1.18682042] [0. 0. 0. 0. 0.54966135]] R8UTT_MTM_TEST R8UTT_MTM computes C = A' * B for R8UTT matrices. The matrix A: Col: 0 1 2 3 4 Row 0 : 0.722083 0.27556 0.026669 0.0823927 0.459709 1 : 0.722083 0.27556 0.026669 0.0823927 2 : 0.722083 0.27556 0.026669 3 : 0.722083 0.27556 4 : 0.722083 The matrix B: Col: 0 1 2 3 4 Row 0 : 0.40691 0.275098 0.49897 0.0869746 0.0369715 1 : 0.40691 0.275098 0.49897 0.0869746 2 : 0.40691 0.275098 0.49897 3 : 0.40691 0.275098 4 : 0.40691 The product C = A' * B: [[0.29382274 0.19864329 0.36029764 0.06280287 0.02669649] [0.1121281 0.36962863 0.33613942 0.38426436 0.07299073] [0.0108519 0.11946469 0.38293568 0.33845895 0.38525035] [0.03352639 0.03351792 0.16057614 0.39010174 0.34150513] [0.18706027 0.1599913 0.26289899 0.20055916 0.40709788]] The R8GE product C = A' * B: [[0.29382274 0.19864329 0.36029764 0.06280287 0.02669649] [0.1121281 0.36962863 0.33613942 0.38426436 0.07299073] [0.0108519 0.11946469 0.38293568 0.33845895 0.38525035] [0.03352639 0.03351792 0.16057614 0.39010174 0.34150513] [0.18706027 0.1599913 0.26289899 0.20055916 0.40709788]] R8UTT_MTV_TEST R8UTT_MTV computes a matrix product b=A'*x for an R8UTT matrix. Matrix order N = 5 The matrix A: Col: 0 1 2 3 4 Row 0 : 1 2 3 4 5 1 : 1 2 3 4 2 : 1 2 3 3 : 1 2 4 : 1 The vector X: 0: 1 1: 2 2: 3 3: 4 4: 5 The vector b=A'*x: 0: 1 1: 4 2: 10 3: 20 4: 35 R8UTT_MV_TEST R8UTT_MV computes a product b=A*x for an R8UTT matrix. Matrix order N = 5 The R8UTT matrix A: Col: 0 1 2 3 4 Row 0 : 1 2 3 4 5 1 : 1 2 3 4 2 : 1 2 3 3 : 1 2 4 : 1 Vector x: 0: 1 1: 2 2: 3 3: 4 4: 5 Vector b = A*x: 0: 55 1: 40 2: 26 3: 14 4: 5 R8UTT_PRINT_TEST R8UTT_PRINT prints an R8UTT matrix. Matrix order N = 5 The matrix: Col: 0 1 2 3 4 Row 0 : 1 2 3 4 5 1 : 1 2 3 4 2 : 1 2 3 3 : 1 2 4 : 1 R8UTT_PRINT_SOME_TEST R8UTT_PRINT_SOME prints some of an R8UTT matrix. Matrix order N = 6 Some of the matrix: Col: 4 5 Row 1 : 4 5 2 : 3 4 3 : 2 3 R8UTT_RANDOM_TEST R8UTT_RANDOM randomizes an R8UTT matrix. Matrix order N = 5 Matrix A: Col: 0 1 2 3 4 Row 0 : 0.692892 0.0900891 0.746049 0.887633 0.420128 1 : 0.692892 0.0900891 0.746049 0.887633 2 : 0.692892 0.0900891 0.746049 3 : 0.692892 0.0900891 4 : 0.692892 R8UTT_SL_TEST R8UTT_SL solves a linear system A*x=b with R8UTT matrix Matrix order N = 5 Matrix A: Col: 0 1 2 3 4 Row 0 : 0.431569 0.938854 0.373905 0.407138 0.76494 1 : 0.431569 0.938854 0.373905 0.407138 2 : 0.431569 0.938854 0.373905 3 : 0.431569 0.938854 4 : 0.431569 Right hand side b: 0: 8.88424 1: 7.21101 2: 6.91965 3: 6.42054 4: 2.15784 Solution x: 0: 1 1: 2 2: 3 3: 4 4: 5 R8UTT_SLT_TEST R8UTT_SLT solves a linear system A'x=b with R8UTT matrix Matrix order N = 5 Matrix A: Col: 0 1 2 3 4 Row 0 : 0.235384 0.904817 0.793565 0.605982 0.732122 1 : 0.235384 0.904817 0.793565 0.605982 2 : 0.235384 0.904817 0.793565 3 : 0.235384 0.904817 4 : 0.235384 Right hand side b: 0: 0.235384 1: 1.37559 2: 3.30935 3: 5.8491 4: 9.12097 Solution x: 0: 1 1: 2 2: 3 3: 4 4: 5 R8UTT_TO_R8GE_TEST R8UTT_TO_R8GE converts an R8UTT matrix to R8GE format. The random R8UTT matrix: Col: 0 1 2 3 4 Row 0 : 0.704708 0.168038 0.999657 0.296149 0.197752 1 : 0.704708 0.168038 0.999657 0.296149 2 : 0.704708 0.168038 0.999657 3 : 0.704708 0.168038 4 : 0.704708 The R8GE matrix: [[0.70470769 0.16803833 0.99965656 0.29614937 0.19775186] [0. 0.70470769 0.16803833 0.99965656 0.29614937] [0. 0. 0.70470769 0.16803833 0.99965656] [0. 0. 0. 0.70470769 0.16803833] [0. 0. 0. 0. 0.70470769]] R8UTT_ZEROS_TEST R8UTT_ZEROS zeros out space for an R8UTT matrix. Matrix order N = 5 The matrix: Col: 0 1 2 3 4 Row 0 : 0 0 0 0 0 1 : 0 0 0 0 2 : 0 0 0 3 : 0 0 4 : 0 r8utt_test(): Normal end of execution. Wed Oct 8 08:50:36 2025