9 April 2023 6:18:16.325 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.218418 0.956318 0.829509 0.561695 0.415307 2 0.218418 0.956318 0.829509 0.561695 3 0.218418 0.956318 0.829509 4 0.218418 0.956318 5 0.218418 The determinant = 0.497102E-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.218418 0.956318 0.829509 0.561695 0.415307 2 0.218418 0.956318 0.829509 0.561695 3 0.218418 0.956318 0.829509 4 0.218418 0.956318 5 0.218418 The inverse matrix B: Col: 1 2 3 4 5 Row --- 1 4.57837 -20.0458 70.3805 -243.796 842.985 2 4.57837 -20.0458 70.3805 -243.796 3 4.57837 -20.0458 70.3805 4 4.57837 -20.0458 5 4.57837 The product A * B: Col: 1 2 3 4 5 Row --- 1 1.00000 0.888178E-15 0.444089E-14 0.355271E-14 -0.393019E-13 2 1.00000 0.888178E-15 0.444089E-14 0.355271E-14 3 1.00000 0.888178E-15 0.444089E-14 4 1.00000 0.888178E-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.218418 0.956318 0.829509 0.561695 0.415307 2 0.218418 0.956318 0.829509 0.561695 3 0.218418 0.956318 0.829509 4 0.218418 0.956318 5 0.218418 Factor B: Col: 1 2 3 4 5 Row --- 1 0.661187E-01 0.257578 0.109957 0.438290E-01 0.633966 2 0.661187E-01 0.257578 0.109957 0.438290E-01 3 0.661187E-01 0.257578 0.109957 4 0.661187E-01 0.257578 5 0.661187E-01 The product C = A * B Col: 1 2 3 4 5 Row --- 1 0.144415E-01 0.119490 0.325189 0.365528 0.443734 2 0.144415E-01 0.119490 0.325189 0.365528 3 0.144415E-01 0.119490 0.325189 4 0.144415E-01 0.119490 5 0.144415E-01 The R8GE product C: Col: 1 2 3 4 5 Row --- 1 0.144415E-01 0.119490 0.325189 0.365528 0.443734 2 0.00000 0.144415E-01 0.119490 0.325189 0.365528 3 0.00000 0.00000 0.144415E-01 0.119490 0.325189 4 0.00000 0.00000 0.00000 0.144415E-01 0.119490 5 0.00000 0.00000 0.00000 0.00000 0.144415E-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.218418 0.956318 0.829509 0.561695 0.415307 2 0.218418 0.956318 0.829509 0.561695 3 0.218418 0.956318 0.829509 4 0.218418 0.956318 5 0.218418 The matrix B: Col: 1 2 3 4 5 Row --- 1 0.661187E-01 0.257578 0.109957 0.438290E-01 0.633966 2 0.661187E-01 0.257578 0.109957 0.438290E-01 3 0.661187E-01 0.257578 0.109957 4 0.661187E-01 0.257578 5 0.661187E-01 The product C = A' * B: Col: 1 2 3 4 5 Row --- 1 0.144415E-01 0.562597E-01 0.240166E-01 0.957306E-02 0.138470 2 0.632305E-01 0.260768 0.161413 0.659310E-01 0.615846 3 0.548461E-01 0.276894 0.351978 0.197770 0.591811 4 0.371386E-01 0.199526 0.338656 0.376596 0.553866 5 0.274596E-01 0.144112 0.245192 0.356858 0.639887 The R8GE product C = A' * B: Col: 1 2 3 4 5 Row --- 1 0.144415E-01 0.562597E-01 0.240166E-01 0.957306E-02 0.138470 2 0.632305E-01 0.260768 0.161413 0.659310E-01 0.615846 3 0.548461E-01 0.276894 0.351978 0.197770 0.591811 4 0.371386E-01 0.199526 0.338656 0.376596 0.553866 5 0.274596E-01 0.144112 0.245192 0.356858 0.639887 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.218418 0.956318 0.829509 0.561695 0.415307 2 0.218418 0.956318 0.829509 0.561695 3 0.218418 0.956318 0.829509 4 0.218418 0.956318 5 0.218418 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.218418 0.956318 0.829509 0.561695 0.415307 2 0.218418 0.956318 0.829509 0.561695 3 0.218418 0.956318 0.829509 4 0.218418 0.956318 5 0.218418 Right hand side b: 1: 8.9428983 2: 9.4323035 3: 8.6280714 4: 5.6552611 5: 1.0920915 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.218418 0.956318 0.829509 0.561695 0.415307 2 0.218418 0.956318 0.829509 0.561695 3 0.218418 0.956318 0.829509 4 0.218418 0.956318 5 0.218418 Right hand side b: 1: 0.21841830 2: 1.3931542 3: 3.3973993 4: 5.9633398 5: 8.9445875 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.218418 0.956318 0.829509 0.561695 0.415307 2 0.218418 0.956318 0.829509 0.561695 3 0.218418 0.956318 0.829509 4 0.218418 0.956318 5 0.218418 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.956318 0.829509 0.561695 0.415307 2 0.00000 0.218418 0.956318 0.829509 0.561695 3 0.00000 0.00000 0.218418 0.956318 0.829509 4 0.00000 0.00000 0.00000 0.218418 0.956318 5 0.00000 0.00000 0.00000 0.00000 0.218418 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 April 2023 6:18:16.326 PM