Sun Aug 28 17:17:12 2022 r8utt_test(): Python version: 3.6.9 Test r8utt(). R8UTT_DET_TEST R8UTT_DET computes the determinant of an R8UTT matrix. Matrix order N = 5 The matrix: Col: 0 1 2 3 4 Row 0 : 0.143871 0.874915 0.880852 0.753068 0.815347 1 : 0.143871 0.874915 0.880852 0.753068 2 : 0.143871 0.874915 0.880852 3 : 0.143871 0.874915 4 : 0.143871 The determinant = 6.164142590535286e-05 R8UTT_INDICATOR_TEST R8UTT_INDICATOR sets up an indicator matrix in R8UTT format Matrix order N = 5 The indicator matrix: Col: 0 1 2 3 4 Row 0 : 1 2 3 4 5 1 : 1 2 3 4 2 : 1 2 3 3 : 1 2 4 : 1 R8UTT_INVERSE_TEST R8UTT_INVERSE computes the inverse of an R8UTT matrix. The matrix A: Col: 0 1 2 3 4 Row 0 : 0.204299 0.937634 0.821366 0.820244 0.883228 1 : 0.204299 0.937634 0.821366 0.820244 2 : 0.204299 0.937634 0.821366 3 : 0.204299 0.937634 4 : 0.204299 The inverse matrix B: Col: 0 1 2 3 4 Row 0 : 4.89478 -22.4646 83.4228 -312.205 1166.51 1 : 4.89478 -22.4646 83.4228 -312.205 2 : 4.89478 -22.4646 83.4228 3 : 4.89478 -22.4646 4 : 4.89478 The product A * B: Col: 0 1 2 3 4 Row 0 : 1 0 8.88178e-16 -1.77636e-15 -6.75016e-14 1 : 1 0 8.88178e-16 -1.77636e-15 2 : 1 0 8.88178e-16 3 : 1 0 4 : 1 R8UTT_MM_TEST R8UTT_MM computes C = A * B for R8UTT matrices. Matrix order N = 5 Factor A: Col: 0 1 2 3 4 Row 0 : 0.183886 0.74717 0.0982378 0.0951454 0.565996 1 : 0.183886 0.74717 0.0982378 0.0951454 2 : 0.183886 0.74717 0.0982378 3 : 0.183886 0.74717 4 : 0.183886 Factor B: Col: 0 1 2 3 4 Row 0 : 0.579224 0.751624 0.177652 0.897127 0.341144 1 : 0.579224 0.751624 0.177652 0.897127 2 : 0.579224 0.751624 0.177652 3 : 0.579224 0.751624 4 : 0.579224 The product C = A * B Col: 0 1 2 3 4 Row 0 : 0.106511 0.570992 0.65116 0.426654 1.14984 1 : 0.106511 0.570992 0.65116 0.426654 2 : 0.106511 0.570992 0.65116 3 : 0.106511 0.570992 4 : 0.106511 The R8GE product C: [[0.10651138 0.57099225 0.65116017 0.42665373 1.14984243] [0. 0.10651138 0.57099225 0.65116017 0.42665373] [0. 0. 0.10651138 0.57099225 0.65116017] [0. 0. 0. 0.10651138 0.57099225] [0. 0. 0. 0. 0.10651138]] R8UTT_MTM_TEST R8UTT_MTM computes C = A' * B for R8UTT matrices. The matrix A: Col: 0 1 2 3 4 Row 0 : 0.728414 0.168842 0.850197 0.524708 0.763983 1 : 0.728414 0.168842 0.850197 0.524708 2 : 0.728414 0.168842 0.850197 3 : 0.728414 0.168842 4 : 0.728414 The matrix B: Col: 0 1 2 3 4 Row 0 : 0.341229 0.442786 0.819831 0.370099 0.392842 1 : 0.341229 0.442786 0.819831 0.370099 2 : 0.341229 0.442786 0.819831 3 : 0.341229 0.442786 4 : 0.341229 The product C = A' * B: [[0.24855587 0.32253178 0.59717642 0.26958527 0.28615193] [0.05761387 0.32331694 0.46095395 0.65966478 0.33591369] [0.29011152 0.43406921 1.0203343 0.77561073 0.99365798] [0.17904559 0.52244516 0.86424131 1.21452827 0.98173839] [0.26069295 0.51732672 1.14878193 1.14699052 1.51465311]] The R8GE product C = A' * B: [[0.24855587 0.32253178 0.59717642 0.26958527 0.28615193] [0.05761387 0.32331694 0.46095395 0.65966478 0.33591369] [0.29011152 0.43406921 1.0203343 0.77561073 0.99365798] [0.17904559 0.52244516 0.86424131 1.21452827 0.98173839] [0.26069295 0.51732672 1.14878193 1.14699052 1.51465311]] R8UTT_MTV_TEST R8UTT_MTV computes a matrix product b=A'*x for an R8UTT matrix. Matrix order N = 5 The matrix A: Col: 0 1 2 3 4 Row 0 : 1 2 3 4 5 1 : 1 2 3 4 2 : 1 2 3 3 : 1 2 4 : 1 The vector X: 0: 1 1: 2 2: 3 3: 4 4: 5 The vector b=A'*x: 0: 1 1: 4 2: 10 3: 20 4: 35 R8UTT_MV_TEST R8UTT_MV computes a product b=A*x for an R8UTT matrix. Matrix order N = 5 The R8UTT matrix A: Col: 0 1 2 3 4 Row 0 : 1 2 3 4 5 1 : 1 2 3 4 2 : 1 2 3 3 : 1 2 4 : 1 Vector x: 0: 1 1: 2 2: 3 3: 4 4: 5 Vector b = A*x: 0: 55 1: 40 2: 26 3: 14 4: 5 R8UTT_PRINT_TEST R8UTT_PRINT prints an R8UTT matrix. Matrix order N = 5 The matrix: Col: 0 1 2 3 4 Row 0 : 1 2 3 4 5 1 : 1 2 3 4 2 : 1 2 3 3 : 1 2 4 : 1 R8UTT_PRINT_SOME_TEST R8UTT_PRINT_SOME prints some of an R8UTT matrix. Matrix order N = 6 Some of the matrix: Col: 4 5 Row 1 : 4 5 2 : 3 4 3 : 2 3 R8UTT_RANDOM_TEST R8UTT_RANDOM randomizes an R8UTT matrix. Matrix order N = 5 Matrix A: Col: 0 1 2 3 4 Row 0 : 0.516423 0.900901 0.96679 0.414407 0.175927 1 : 0.516423 0.900901 0.96679 0.414407 2 : 0.516423 0.900901 0.96679 3 : 0.516423 0.900901 4 : 0.516423 R8UTT_SL_TEST R8UTT_SL solves a linear system A*x=b with R8UTT matrix Matrix order N = 5 Matrix A: Col: 0 1 2 3 4 Row 0 : 0.430306 0.0138733 0.753348 0.542658 0.709931 1 : 0.430306 0.0138733 0.753348 0.542658 2 : 0.430306 0.0138733 0.753348 3 : 0.430306 0.0138733 4 : 0.430306 Right hand side b: 0: 8.43839 1: 6.62892 2: 5.11315 3: 1.79059 4: 2.15153 Solution x: 0: 1 1: 2 2: 3 3: 4 4: 5 R8UTT_SLT_TEST R8UTT_SLT solves a linear system A'x=b with R8UTT matrix Matrix order N = 5 Matrix A: Col: 0 1 2 3 4 Row 0 : 0.156754 0.272812 0.496669 0.403092 0.363841 1 : 0.156754 0.272812 0.496669 0.403092 2 : 0.156754 0.272812 0.496669 3 : 0.156754 0.272812 4 : 0.156754 Right hand side b: 0: 0.156754 1: 0.58632 2: 1.51255 3: 2.84188 4: 4.53505 Solution x: 0: 1 1: 2 2: 3 3: 4 4: 5 R8UTT_TO_R8GE_TEST R8UTT_TO_R8GE converts an R8UTT matrix to R8GE format. The random R8UTT matrix: Col: 0 1 2 3 4 Row 0 : 0.143453 0.697558 0.0251148 0.409025 0.0918546 1 : 0.143453 0.697558 0.0251148 0.409025 2 : 0.143453 0.697558 0.0251148 3 : 0.143453 0.697558 4 : 0.143453 The R8GE matrix: [[0.14345314 0.69755819 0.02511483 0.40902538 0.09185461] [0. 0.14345314 0.69755819 0.02511483 0.40902538] [0. 0. 0.14345314 0.69755819 0.02511483] [0. 0. 0. 0.14345314 0.69755819] [0. 0. 0. 0. 0.14345314]] R8UTT_ZEROS_TEST R8UTT_ZEROS zeros out space for an R8UTT matrix. Matrix order N = 5 The matrix: Col: 0 1 2 3 4 Row 0 : 0 0 0 0 0 1 : 0 0 0 0 2 : 0 0 0 3 : 0 0 4 : 0 r8utt_test(): Normal end of execution. Sun Aug 28 17:17:12 2022