11 April 2023 7:10:46.664 AM 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.218418 2 0.956318 0.218418 3 0.829509 0.956318 0.218418 4 0.561695 0.829509 0.956318 0.218418 5 0.415307 0.561695 0.829509 0.956318 0.218418 The determinant = 0.497102E-03 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.218418 2 0.956318 0.218418 3 0.829509 0.956318 0.218418 4 0.561695 0.829509 0.956318 0.218418 5 0.415307 0.561695 0.829509 0.956318 0.218418 The inverse matrix B: Col: 1 2 3 4 5 Row --- 1 4.57837 2 -20.0458 4.57837 3 70.3805 -20.0458 4.57837 4 -243.796 70.3805 -20.0458 4.57837 5 842.985 -243.796 70.3805 -20.0458 4.57837 The product A * B: Col: 1 2 3 4 5 Row --- 1 1.00000 2 0.888178E-15 1.00000 3 0.444089E-14 0.888178E-15 1.00000 4 0.355271E-14 0.444089E-14 0.888178E-15 1.00000 5 -0.393019E-13 0.355271E-14 0.444089E-14 0.888178E-15 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.218418 2 0.956318 0.218418 3 0.829509 0.956318 0.218418 4 0.561695 0.829509 0.956318 0.218418 5 0.415307 0.561695 0.829509 0.956318 0.218418 Factor B: Col: 1 2 3 4 5 Row --- 1 0.661187E-01 2 0.257578 0.661187E-01 3 0.109957 0.257578 0.661187E-01 4 0.438290E-01 0.109957 0.257578 0.661187E-01 5 0.633966 0.438290E-01 0.109957 0.257578 0.661187E-01 The product C = A * B Col: 1 2 3 4 5 Row --- 1 0.144415E-01 2 0.119490 0.144415E-01 3 0.325189 0.119490 0.144415E-01 4 0.365528 0.325189 0.119490 0.144415E-01 5 0.443734 0.365528 0.325189 0.119490 0.144415E-01 The R8GE product C: Col: 1 2 3 4 5 Row --- 1 0.144415E-01 0.00000 0.00000 0.00000 0.00000 2 0.119490 0.144415E-01 0.00000 0.00000 0.00000 3 0.325189 0.119490 0.144415E-01 0.00000 0.00000 4 0.365528 0.325189 0.119490 0.144415E-01 0.00000 5 0.443734 0.365528 0.325189 0.119490 0.144415E-01 R8LTT_MTM_TEST R8LTT_MTM computes C = A' * B for R8LTT matrices. The matrix A: Col: 1 2 3 4 5 Row --- 1 0.218418 2 0.956318 0.218418 3 0.829509 0.956318 0.218418 4 0.561695 0.829509 0.956318 0.218418 5 0.415307 0.561695 0.829509 0.956318 0.218418 The matrix B: Col: 1 2 3 4 5 Row --- 1 0.661187E-01 2 0.257578 0.661187E-01 3 0.109957 0.257578 0.661187E-01 4 0.438290E-01 0.109957 0.257578 0.661187E-01 5 0.633966 0.438290E-01 0.109957 0.257578 0.661187E-01 The product C = A' * B: Col: 1 2 3 4 5 Row --- 1 0.639887 0.356858 0.245192 0.144112 0.274596E-01 2 0.553866 0.376596 0.338656 0.199526 0.371386E-01 3 0.591811 0.197770 0.351978 0.276894 0.548461E-01 4 0.615846 0.659310E-01 0.161413 0.260768 0.632305E-01 5 0.138470 0.957306E-02 0.240166E-01 0.562597E-01 0.144415E-01 The R8GE product C = A' * B: Col: 1 2 3 4 5 Row --- 1 0.639887 0.356858 0.245192 0.144112 0.274596E-01 2 0.553866 0.376596 0.338656 0.199526 0.371386E-01 3 0.591811 0.197770 0.351978 0.276894 0.548461E-01 4 0.615846 0.659310E-01 0.161413 0.260768 0.632305E-01 5 0.138470 0.957306E-02 0.240166E-01 0.562597E-01 0.144415E-01 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.218418 2 0.956318 0.218418 3 0.829509 0.956318 0.218418 4 0.561695 0.829509 0.956318 0.218418 5 0.415307 0.561695 0.829509 0.956318 0.218418 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.218418 2 0.956318 0.218418 3 0.829509 0.956318 0.218418 4 0.561695 0.829509 0.956318 0.218418 5 0.415307 0.561695 0.829509 0.956318 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 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.218418 2 0.956318 0.218418 3 0.829509 0.956318 0.218418 4 0.561695 0.829509 0.956318 0.218418 5 0.415307 0.561695 0.829509 0.956318 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 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.218418 2 0.956318 0.218418 3 0.829509 0.956318 0.218418 4 0.561695 0.829509 0.956318 0.218418 5 0.415307 0.561695 0.829509 0.956318 0.218418 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.00000 0.00000 0.00000 0.00000 2 0.956318 0.218418 0.00000 0.00000 0.00000 3 0.829509 0.956318 0.218418 0.00000 0.00000 4 0.561695 0.829509 0.956318 0.218418 0.00000 5 0.415307 0.561695 0.829509 0.956318 0.218418 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. 11 April 2023 7:10:46.665 AM