6 October 2025 6:33:38.729 PM r8utt_test(): Fortran90 version Test r8utt(). R8UTT_DET_TEST R8UTT_DET computes the determinant of an R8UTT matrix. Matrix order N = 5 The matrix: Col: 1 2 3 4 5 Row --- 1 0.699719 0.882595 0.425188 0.780196 0.773669 2 0.699719 0.882595 0.425188 0.780196 3 0.699719 0.882595 0.425188 4 0.699719 0.882595 5 0.699719 The determinant = 0.167733 R8UTT_INDICATOR_TEST R8UTT_INDICATOR sets up an indicator matrix in R8UTT format Matrix order N = 5 The indicator matrix: Col: 1 2 3 4 5 Row --- 1 1.00000 2.00000 3.00000 4.00000 5.00000 2 1.00000 2.00000 3.00000 4.00000 3 1.00000 2.00000 3.00000 4 1.00000 2.00000 5 1.00000 R8UTT_INVERSE_TEST R8UTT_INVERSE computes the inverse of an R8UTT matrix. The matrix A: Col: 1 2 3 4 5 Row --- 1 0.397288 0.679453 0.669322 0.412670 0.826119 2 0.397288 0.679453 0.669322 0.412670 3 0.397288 0.679453 0.669322 4 0.397288 0.679453 5 0.397288 The inverse matrix B: Col: 1 2 3 4 5 Row --- 1 2.51706 -4.30475 3.12153 -0.700718 -4.82309 2 2.51706 -4.30475 3.12153 -0.700718 3 2.51706 -4.30475 3.12153 4 2.51706 -4.30475 5 2.51706 The product A * B: Col: 1 2 3 4 5 Row --- 1 1.00000 0.222045E-15 0.222045E-15 -0.222045E-15 -0.444089E-15 2 1.00000 0.222045E-15 0.222045E-15 -0.222045E-15 3 1.00000 0.222045E-15 0.222045E-15 4 1.00000 0.222045E-15 5 1.00000 R8UTT_MM_TEST R8UTT_MM computes C = A * B for R8UTT matrices. Matrix order N = 5 Factor A: Col: 1 2 3 4 5 Row --- 1 0.458990 0.737217 0.339784E-01 0.592138 0.141770 2 0.458990 0.737217 0.339784E-01 0.592138 3 0.458990 0.737217 0.339784E-01 4 0.458990 0.737217 5 0.458990 Factor B: Col: 1 2 3 4 5 Row --- 1 0.982875 0.670549 0.708941 0.408568 0.342234E-01 2 0.982875 0.670549 0.708941 0.408568 3 0.982875 0.670549 0.708941 4 0.982875 0.670549 5 0.982875 The product C = A * B Col: 1 2 3 4 5 Row --- 1 0.451130 1.03237 0.853133 1.31495 0.877400 2 0.451130 1.03237 0.853133 1.31495 3 0.451130 1.03237 0.853133 4 0.451130 1.03237 5 0.451130 The R8GE product C: Col: 1 2 3 4 5 Row --- 1 0.451130 1.03237 0.853133 1.31495 0.877400 2 0.00000 0.451130 1.03237 0.853133 1.31495 3 0.00000 0.00000 0.451130 1.03237 0.853133 4 0.00000 0.00000 0.00000 0.451130 1.03237 5 0.00000 0.00000 0.00000 0.00000 0.451130 R8UTT_MTM_TEST R8UTT_MTM computes C = A' * B for R8UTT matrices. The matrix A: Col: 1 2 3 4 5 Row --- 1 0.492653 0.847935 0.765086 0.621955 0.369518 2 0.492653 0.847935 0.765086 0.621955 3 0.492653 0.847935 0.765086 4 0.492653 0.847935 5 0.492653 The matrix B: Col: 1 2 3 4 5 Row --- 1 0.131402 0.958146 0.642084 0.679178 0.585002 2 0.131402 0.958146 0.642084 0.679178 3 0.131402 0.958146 0.642084 4 0.131402 0.958146 5 0.131402 The product C = A' * B: Col: 1 2 3 4 5 Row --- 1 0.647355E-01 0.472034 0.316325 0.334599 0.288203 2 0.111420 0.877181 1.01648 0.892224 0.830643 3 0.100534 0.844485 1.36843 1.53611 1.33980 4 0.817261E-01 0.696457 1.24383 1.79085 1.89995 5 0.485553E-01 0.435778 0.933719 1.49480 2.00702 The R8GE product C = A' * B: Col: 1 2 3 4 5 Row --- 1 0.647355E-01 0.472034 0.316325 0.334599 0.288203 2 0.111420 0.877181 1.01648 0.892224 0.830643 3 0.100534 0.844485 1.36843 1.53611 1.33980 4 0.817261E-01 0.696457 1.24383 1.79085 1.89995 5 0.485553E-01 0.435778 0.933719 1.49480 2.00702 R8UTT_MTV_TEST R8UTT_MTV computes a matrix product b=A'*x for an R8UTT matrix. Matrix order N = 5 The matrix A: Col: 1 2 3 4 5 Row --- 1 1.00000 2.00000 3.00000 4.00000 5.00000 2 1.00000 2.00000 3.00000 4.00000 3 1.00000 2.00000 3.00000 4 1.00000 2.00000 5 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: 1.0000000 2: 4.0000000 3: 10.000000 4: 20.000000 5: 35.000000 R8UTT_MV_TEST R8UTT_MV computes a product b=A*x for an R8UTT matrix. Matrix order N = 5 The R8UTT matrix A: Col: 1 2 3 4 5 Row --- 1 1.00000 2.00000 3.00000 4.00000 5.00000 2 1.00000 2.00000 3.00000 4.00000 3 1.00000 2.00000 3.00000 4 1.00000 2.00000 5 1.00000 Vector x: 1: 1.0000000 2: 2.0000000 3: 3.0000000 4: 4.0000000 5: 5.0000000 Vector b = A*x: 1: 55.000000 2: 40.000000 3: 26.000000 4: 14.000000 5: 5.0000000 R8UTT_PRINT_TEST R8UTT_PRINT prints an R8UTT matrix. Matrix order N = 5 The matrix: Col: 1 2 3 4 5 Row --- 1 1.00000 2.00000 3.00000 4.00000 5.00000 2 1.00000 2.00000 3.00000 4.00000 3 1.00000 2.00000 3.00000 4 1.00000 2.00000 5 1.00000 R8UTT_PRINT_SOME_TEST R8UTT_PRINT_SOME prints some of an R8UTT matrix. Matrix order N = 6 Some of the matrix: Col: 4 5 6 Row --- 1 4.00000 5.00000 6.00000 2 3.00000 4.00000 5.00000 3 2.00000 3.00000 4.00000 R8UTT_RANDOM_TEST R8UTT_RANDOM randomizes an R8UTT matrix. Matrix order N = 5 Matrix A: Col: 1 2 3 4 5 Row --- 1 0.455444 0.901416 0.419902E-01 0.130348 0.304037 2 0.455444 0.901416 0.419902E-01 0.130348 3 0.455444 0.901416 0.419902E-01 4 0.455444 0.901416 5 0.455444 R8UTT_SL_TEST R8UTT_SL solves a linear system A*x=b with R8UTT matrix Matrix order N = 5 Matrix A: Col: 1 2 3 4 5 Row --- 1 0.263334 0.889342 0.179393 0.151212 0.974111 2 0.263334 0.889342 0.179393 0.151212 3 0.263334 0.889342 0.179393 4 0.263334 0.889342 5 0.263334 Right hand side b: 1: 8.0556034 2: 4.6683285 3: 5.2443378 4: 5.5000471 5: 1.3166690 Solution x: 1: 1.0000000 2: 2.0000000 3: 3.0000000 4: 4.0000000 5: 5.0000000 R8UTT_SLT_TEST R8UTT_SLT solves a linear system A'x=b with R8UTT matrix Matrix order N = 5 Matrix A: Col: 1 2 3 4 5 Row --- 1 0.125564 0.799404 0.873440 0.925432 0.937656 2 0.125564 0.799404 0.873440 0.925432 3 0.125564 0.799404 0.873440 4 0.125564 0.799404 5 0.125564 Right hand side b: 1: 0.12556441 2: 1.0505327 3: 2.8489412 4: 5.5727820 5: 9.2342791 Solution x: 1: 1.0000000 2: 2.0000000 3: 3.0000000 4: 4.0000000 5: 5.0000000 R8UTT_TO_R8GE_TEST R8UTT_TO_R8GE converts an R8UTT matrix to R8GE format. The random R8UTT matrix: Col: 1 2 3 4 5 Row --- 1 0.753941 0.259612E-01 0.531844 0.636332 0.487880 2 0.753941 0.259612E-01 0.531844 0.636332 3 0.753941 0.259612E-01 0.531844 4 0.753941 0.259612E-01 5 0.753941 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 0.753941 0.259612E-01 0.531844 0.636332 0.487880 2 0.00000 0.753941 0.259612E-01 0.531844 0.636332 3 0.00000 0.00000 0.753941 0.259612E-01 0.531844 4 0.00000 0.00000 0.00000 0.753941 0.259612E-01 5 0.00000 0.00000 0.00000 0.00000 0.753941 R8UTT_ZEROS_TEST R8UTT_ZEROS zeros out space for an R8UTT matrix. Matrix order N = 5 The matrix: Col: 1 2 3 4 5 Row --- 1 0.00000 0.00000 0.00000 0.00000 0.00000 2 0.00000 0.00000 0.00000 0.00000 3 0.00000 0.00000 0.00000 4 0.00000 0.00000 5 0.00000 r8utt_test(): Normal end of execution. 6 October 2025 6:33:38.730 PM