9 May 2025 9:29:17.371 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.158819 0.553677 0.444242E-01 0.122193 0.558862 2 0.158819 0.553677 0.444242E-01 0.122193 3 0.158819 0.553677 0.444242E-01 4 0.158819 0.553677 5 0.158819 The determinant = 0.101045E-03 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.363223 0.870941 0.602289 0.749382 0.751923 2 0.363223 0.870941 0.602289 0.749382 3 0.363223 0.870941 0.602289 4 0.363223 0.870941 5 0.363223 The inverse matrix B: Col: 1 2 3 4 5 Row --- 1 2.75313 -6.60149 11.2640 -21.7425 41.3773 2 2.75313 -6.60149 11.2640 -21.7425 3 2.75313 -6.60149 11.2640 4 2.75313 -6.60149 5 2.75313 The product A * B: Col: 1 2 3 4 5 Row --- 1 1.00000 0.00000 0.133227E-14 0.177636E-14 -0.133227E-14 2 1.00000 0.00000 0.133227E-14 0.177636E-14 3 1.00000 0.00000 0.133227E-14 4 1.00000 0.00000 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.178132 0.365374 0.135288 0.262140 0.607356 2 0.178132 0.365374 0.135288 0.262140 3 0.178132 0.365374 0.135288 4 0.178132 0.365374 5 0.178132 Factor B: Col: 1 2 3 4 5 Row --- 1 0.322378 0.584637 0.297625 0.610535 0.307008 2 0.322378 0.584637 0.297625 0.610535 3 0.322378 0.584637 0.297625 4 0.322378 0.584637 5 0.322378 The product C = A * B Col: 1 2 3 4 5 Row --- 1 0.574257E-01 0.221931 0.310242 0.381103 0.667082 2 0.574257E-01 0.221931 0.310242 0.381103 3 0.574257E-01 0.221931 0.310242 4 0.574257E-01 0.221931 5 0.574257E-01 The R8GE product C: Col: 1 2 3 4 5 Row --- 1 0.574257E-01 0.221931 0.310242 0.381103 0.667082 2 0.00000 0.574257E-01 0.221931 0.310242 0.381103 3 0.00000 0.00000 0.574257E-01 0.221931 0.310242 4 0.00000 0.00000 0.00000 0.574257E-01 0.221931 5 0.00000 0.00000 0.00000 0.00000 0.574257E-01 R8UTT_MTM_TEST R8UTT_MTM computes C = A' * B for R8UTT matrices. The matrix A: Col: 1 2 3 4 5 Row --- 1 0.700750 0.517856 0.882245 0.610530 0.218629 2 0.700750 0.517856 0.882245 0.610530 3 0.700750 0.517856 0.882245 4 0.700750 0.517856 5 0.700750 The matrix B: Col: 1 2 3 4 5 Row --- 1 0.999588 0.628677 0.796685 0.449351 0.468521 2 0.999588 0.628677 0.796685 0.449351 3 0.999588 0.628677 0.796685 4 0.999588 0.628677 5 0.999588 The product C = A' * B: Col: 1 2 3 4 5 Row --- 1 0.700461 0.440545 0.558277 0.314883 0.328316 2 0.517642 1.02602 0.853113 0.790976 0.557509 3 0.881881 1.07229 1.72890 1.24955 1.20433 4 0.610278 1.26571 1.55869 2.00324 1.53560 5 0.218539 0.747725 1.43989 1.65693 2.10567 The R8GE product C = A' * B: Col: 1 2 3 4 5 Row --- 1 0.700461 0.440545 0.558277 0.314883 0.328316 2 0.517642 1.02602 0.853113 0.790976 0.557509 3 0.881881 1.07229 1.72890 1.24955 1.20433 4 0.610278 1.26571 1.55869 2.00324 1.53560 5 0.218539 0.747725 1.43989 1.65693 2.10567 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.437482 0.587146E-01 0.848947 0.799485 0.813103 2 0.437482 0.587146E-01 0.848947 0.799485 3 0.437482 0.587146E-01 0.848947 4 0.437482 0.587146E-01 5 0.437482 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.438527 0.390198 0.869522E-01 0.772901 0.339886 2 0.438527 0.390198 0.869522E-01 0.772901 3 0.438527 0.390198 0.869522E-01 4 0.438527 0.390198 5 0.438527 Right hand side b: 1: 6.2708148 2: 6.2599617 3: 3.3111325 4: 3.7050958 5: 2.1926342 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.569718 0.488307 0.356992 0.945980 0.132869 2 0.569718 0.488307 0.356992 0.945980 3 0.569718 0.488307 0.356992 4 0.569718 0.488307 5 0.569718 Right hand side b: 1: 0.56971843 2: 1.6277434 3: 3.0427600 4: 5.4037566 5: 7.8976224 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.772510 0.113059 0.274241 0.124476 0.729742 2 0.772510 0.113059 0.274241 0.124476 3 0.772510 0.113059 0.274241 4 0.772510 0.113059 5 0.772510 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 0.772510 0.113059 0.274241 0.124476 0.729742 2 0.00000 0.772510 0.113059 0.274241 0.124476 3 0.00000 0.00000 0.772510 0.113059 0.274241 4 0.00000 0.00000 0.00000 0.772510 0.113059 5 0.00000 0.00000 0.00000 0.00000 0.772510 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. 9 May 2025 9:29:17.372 PM