Tue May 20 22:27:39 2025 r8utt_test(): python version: 3.10.12 numpy version: 1.26.4 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.0655857 0.52999 0.398779 0.176295 0.674613 1 : 0.0655857 0.52999 0.398779 0.176295 2 : 0.0655857 0.52999 0.398779 3 : 0.0655857 0.52999 4 : 0.0655857 The determinant = 1.2135175586941965e-06 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.651469 0.296479 0.200729 0.775437 0.461497 1 : 0.651469 0.296479 0.200729 0.775437 2 : 0.651469 0.296479 0.200729 3 : 0.651469 0.296479 4 : 0.651469 The inverse matrix B: Col: 0 1 2 3 4 Row 0 : 1.53499 -0.698564 -0.155046 -1.54129 0.493315 1 : 1.53499 -0.698564 -0.155046 -1.54129 2 : 1.53499 -0.698564 -0.155046 3 : 1.53499 -0.698564 4 : 1.53499 The product A * B: Col: 0 1 2 3 4 Row 0 : 1 5.55112e-17 5.55112e-17 2.22045e-16 -2.22045e-16 1 : 1 5.55112e-17 5.55112e-17 2.22045e-16 2 : 1 5.55112e-17 5.55112e-17 3 : 1 5.55112e-17 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.59018 0.253405 0.46169 0.332474 0.28249 1 : 0.59018 0.253405 0.46169 0.332474 2 : 0.59018 0.253405 0.46169 3 : 0.59018 0.253405 4 : 0.59018 Factor B: Col: 0 1 2 3 4 Row 0 : 0.0524534 0.857983 0.416354 0.00756683 0.474943 1 : 0.0524534 0.857983 0.416354 0.00756683 2 : 0.0524534 0.857983 0.416354 3 : 0.0524534 0.857983 4 : 0.0524534 The product C = A * B Col: 0 1 2 3 4 Row 0 : 0.030957 0.519656 0.487359 0.523534 0.77452 1 : 0.030957 0.519656 0.487359 0.523534 2 : 0.030957 0.519656 0.487359 3 : 0.030957 0.519656 4 : 0.030957 The R8GE product C: [[0.03095696 0.51965603 0.48735856 0.52353371 0.77452045] [0. 0.03095696 0.51965603 0.48735856 0.52353371] [0. 0. 0.03095696 0.51965603 0.48735856] [0. 0. 0. 0.03095696 0.51965603] [0. 0. 0. 0. 0.03095696]] R8UTT_MTM_TEST R8UTT_MTM computes C = A' * B for R8UTT matrices. The matrix A: Col: 0 1 2 3 4 Row 0 : 0.901849 0.727265 0.0564127 0.921918 0.795869 1 : 0.901849 0.727265 0.0564127 0.921918 2 : 0.901849 0.727265 0.0564127 3 : 0.901849 0.727265 4 : 0.901849 The matrix B: Col: 0 1 2 3 4 Row 0 : 0.16879 0.130732 0.0813415 0.105991 0.119026 1 : 0.16879 0.130732 0.0813415 0.105991 2 : 0.16879 0.130732 0.0813415 3 : 0.16879 0.130732 4 : 0.16879 The product C = A' * B: [[0.15222346 0.11790032 0.07335773 0.09558767 0.10734312] [0.12275542 0.24730016 0.17705717 0.15044115 0.18215087] [0.00952191 0.13013035 0.25188885 0.18303639 0.1571557 ] [0.1556109 0.13004588 0.20512051 0.34960364 0.29276822] [0.13433512 0.25965632 0.19478307 0.28947532 0.44433247]] The R8GE product C = A' * B: [[0.15222346 0.11790032 0.07335773 0.09558767 0.10734312] [0.12275542 0.24730016 0.17705717 0.15044115 0.18215087] [0.00952191 0.13013035 0.25188885 0.18303639 0.1571557 ] [0.1556109 0.13004588 0.20512051 0.34960364 0.29276822] [0.13433512 0.25965632 0.19478307 0.28947532 0.44433247]] 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.507333 0.771027 0.679609 0.115838 0.487847 1 : 0.507333 0.771027 0.679609 0.115838 2 : 0.507333 0.771027 0.679609 3 : 0.507333 0.771027 4 : 0.507333 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.582944 0.633759 0.922475 0.837083 0.104019 1 : 0.582944 0.633759 0.922475 0.837083 2 : 0.582944 0.633759 0.922475 3 : 0.582944 0.633759 4 : 0.582944 Right hand side b: 0: 8.48631 1: 10.9425 2: 8.89624 3: 5.50057 4: 2.91472 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.662645 0.476842 0.294005 0.822177 0.232117 1 : 0.662645 0.476842 0.294005 0.822177 2 : 0.662645 0.476842 0.294005 3 : 0.662645 0.476842 4 : 0.662645 Right hand side b: 0: 0.662645 1: 1.80213 2: 3.23562 3: 5.49129 4: 7.97908 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.580427 0.102304 0.951196 0.0125088 0.651481 1 : 0.580427 0.102304 0.951196 0.0125088 2 : 0.580427 0.102304 0.951196 3 : 0.580427 0.102304 4 : 0.580427 The R8GE matrix: [[0.58042732 0.1023037 0.95119597 0.01250885 0.6514811 ] [0. 0.58042732 0.1023037 0.95119597 0.01250885] [0. 0. 0.58042732 0.1023037 0.95119597] [0. 0. 0. 0.58042732 0.1023037 ] [0. 0. 0. 0. 0.58042732]] 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. Tue May 20 22:27:39 2025