23 August 2022 07:52:06 AM r8ltt_test(): C version Test r8ltt(). R8LTT_DET_TEST R8LTT_DET computes the determinant of an R8LTT matrix. The matrix A: Col: 0 1 2 3 4 Row --- 0 0.218418 1 0.956318 0.218418 2 0.829509 0.956318 0.218418 3 0.561695 0.829509 0.956318 0.218418 4 0.415307 0.561695 0.829509 0.956318 0.218418 Determinant is 0.000497102 R8LTT_INDICATOR_TEST R8LTT_INDICATOR sets up an indicator matrix in R8LTT format Matrix order N = 5 The indicator matrix: Col: 0 1 2 3 4 Row --- 0 1 1 2 1 2 3 2 1 3 4 3 2 1 4 5 4 3 2 1 R8LTT_INVERSE_TEST R8LTT_INVERSE computes the inverse of an R8LTT matrix. The matrix A: Col: 0 1 2 3 4 Row --- 0 0.218418 1 0.956318 0.218418 2 0.829509 0.956318 0.218418 3 0.561695 0.829509 0.956318 0.218418 4 0.415307 0.561695 0.829509 0.956318 0.218418 The inverse matrix B: Col: 0 1 2 3 4 Row --- 0 4.57837 1 -20.0458 4.57837 2 70.3805 -20.0458 4.57837 3 -243.796 70.3805 -20.0458 4.57837 4 842.985 -243.796 70.3805 -20.0458 4.57837 The product A * B: Col: 0 1 2 3 4 Row --- 0 1 1 8.88178e-16 1 2 4.44089e-15 8.88178e-16 1 3 3.55271e-15 4.44089e-15 8.88178e-16 1 4 -3.93019e-14 3.55271e-15 4.44089e-15 8.88178e-16 1 R8LTT_MM_TEST R8LTT_MM computes C = A * B for R8LTT matrices. Matrix order N = 5 Factor A: Col: 0 1 2 3 4 Row --- 0 0.218418 1 0.956318 0.218418 2 0.829509 0.956318 0.218418 3 0.561695 0.829509 0.956318 0.218418 4 0.415307 0.561695 0.829509 0.956318 0.218418 Factor B: Col: 0 1 2 3 4 Row --- 0 0.0661187 1 0.257578 0.0661187 2 0.109957 0.257578 0.0661187 3 0.043829 0.109957 0.257578 0.0661187 4 0.633966 0.043829 0.109957 0.257578 0.0661187 The product C = A * B Col: 0 1 2 3 4 Row --- 0 0.0144415 1 0.11949 0.0144415 2 0.325189 0.11949 0.0144415 3 0.365528 0.325189 0.11949 0.0144415 4 0.443734 0.365528 0.325189 0.11949 0.0144415 The R8GE product C: Col: 1 2 3 4 5 Row --- 1 0.0144415 0 0 0 0 2 0.11949 0.0144415 0 0 0 3 0.325189 0.11949 0.0144415 0 0 4 0.365528 0.325189 0.11949 0.0144415 0 5 0.443734 0.365528 0.325189 0.11949 0.0144415 R8LTT_MTM_TEST R8LTT_MTM computes C = A' * B for R8LTT matrices. The matrix A: Col: 0 1 2 3 4 Row --- 0 0.218418 1 0.956318 0.218418 2 0.829509 0.956318 0.218418 3 0.561695 0.829509 0.956318 0.218418 4 0.415307 0.561695 0.829509 0.956318 0.218418 The matrix B: Col: 0 1 2 3 4 Row --- 0 0.0661187 1 0.257578 0.0661187 2 0.109957 0.257578 0.0661187 3 0.043829 0.109957 0.257578 0.0661187 4 0.633966 0.043829 0.109957 0.257578 0.0661187 The product C = A' * B: Col: 1 2 3 4 5 Row --- 1 0.639887 0.356858 0.245192 0.144112 0.0274596 2 0.553866 0.376596 0.338656 0.199526 0.0371386 3 0.591811 0.19777 0.351978 0.276894 0.0548461 4 0.615846 0.065931 0.161413 0.260768 0.0632305 5 0.13847 0.00957306 0.0240166 0.0562597 0.0144415 The R8GE product C = A' * B: Col: 1 2 3 4 5 Row --- 1 0.639887 0.356858 0.245192 0.144112 0.0274596 2 0.553866 0.376596 0.338656 0.199526 0.0371386 3 0.591811 0.19777 0.351978 0.276894 0.0548461 4 0.615846 0.065931 0.161413 0.260768 0.0632305 5 0.13847 0.00957306 0.0240166 0.0562597 0.0144415 R8LTT_MTV_TEST R8LTT_MTV computes a matrix product b=A'*x for an R8LTT matrix. Matrix order N = 5 The matrix A: Col: 0 1 2 3 4 Row --- 0 1 1 2 1 2 3 2 1 3 4 3 2 1 4 5 4 3 2 1 The vector X: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 The vector b=A'*x: 0 55.000000 1 40.000000 2 26.000000 3 14.000000 4 5.000000 R8LTT_MV_TEST R8LTT_MV computes a product b=A*x for an R8LTT matrix. Matrix order N = 5 The R8LTT matrix A: Col: 0 1 2 3 4 Row --- 0 1 1 2 1 2 3 2 1 3 4 3 2 1 4 5 4 3 2 1 Vector x: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 Vector b = A*x: 0 1.000000 1 4.000000 2 10.000000 3 20.000000 4 35.000000 R8LTT_PRINT_TEST R8LTT_PRINT prints an R8LTT matrix. Matrix order N = 5 The matrix: Col: 0 1 2 3 4 Row --- 0 1 1 2 1 2 3 2 1 3 4 3 2 1 4 5 4 3 2 1 R8LTT_PRINT_SOME_TEST R8LTT_PRINT_SOME prints some of an R8LTT matrix. Matrix order N = 5 Some of the matrix: Col: 0 1 2 3 Row --- 1 2 1 2 3 2 1 3 4 3 2 1 4 5 4 3 2 R8LTT_RANDOM_TEST R8LTT_RANDOM randomizes an R8LTT matrix. Matrix order N = 5 Matrix A: Col: 0 1 2 3 4 Row --- 0 0.218418 1 0.956318 0.218418 2 0.829509 0.956318 0.218418 3 0.561695 0.829509 0.956318 0.218418 4 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: 0 1 2 3 4 Row --- 0 0.218418 1 0.956318 0.218418 2 0.829509 0.956318 0.218418 3 0.561695 0.829509 0.956318 0.218418 4 0.415307 0.561695 0.829509 0.956318 0.218418 Right hand side b: 0 0.218418 1 1.393154 2 3.397399 3 5.963340 4 8.944587 Solution x: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 R8LTT_SLT_TEST R8LTT_SLT solves a linear system A'x=b with R8LTT matrix Matrix order N = 5 Matrix A: Col: 0 1 2 3 4 Row --- 0 0.218418 1 0.956318 0.218418 2 0.829509 0.956318 0.218418 3 0.561695 0.829509 0.956318 0.218418 4 0.415307 0.561695 0.829509 0.956318 0.218418 Right hand side b: 0 8.942898 1 9.432303 2 8.628071 3 5.655261 4 1.092091 Solution x: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 R8LTT_TO_R8GE_TEST R8LTT_TO_R8GE converts an R8LTT matrix to R8GE format. R8LTT matrix: Col: 0 1 2 3 4 Row --- 0 0.218418 1 0.956318 0.218418 2 0.829509 0.956318 0.218418 3 0.561695 0.829509 0.956318 0.218418 4 0.415307 0.561695 0.829509 0.956318 0.218418 R8GE matrix Col: 1 2 3 4 5 Row --- 1 0.218418 0 0 0 0 2 0.956318 0.218418 0 0 0 3 0.829509 0.956318 0.218418 0 0 4 0.561695 0.829509 0.956318 0.218418 0 5 0.415307 0.561695 0.829509 0.956318 0.218418 R8LTT_ZEROS_TEST R8LTT_ZEROS returns a zeroed out R8LTT matrix. Matrix order N = 5 Matrix A: Col: 0 1 2 3 4 Row --- 0 0 1 0 0 2 0 0 0 3 0 0 0 0 4 0 0 0 0 0 r8ltt_test(): Normal end of execution. 23 August 2022 07:52:06 AM