Tue May 20 22:27:34 2025 r8ge_np_test(): python version: 3.10.12 numpy version: 1.26.4 Test r8ge_np(). R8GE_NP_DET_TEST R8GE_NP_DET computes the determinant of a matrix that was factored by R8GE_NP_FA, Matrix order N = 10 Determinant of -1, 2, -1 matrix is 10.999999999999998 Exact value is 11 R8GE_NP_FA_TEST R8GE_NP_FA LU factors an R8GE matrix without pivoting. Matrix order N = 10 Solution: 0: 1 1: 1 2: 1 3: 1 4: 1 5: 1 6: 1 7: 1 8: 1 9: 1 Solution: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 8: 9 9: 10 Solution of transposed system: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 8: 9 9: 10 R8GE_NP_INVERSE_TEST R8GE_NP_INVERSE computes the inverse of a matrix that was factored by R8GE_NP_FA. Matrix order N = 5 The R8GE random matrix, A [[0.96997682 0.17482496 0.10629471 0.22542441 0.38523807] [0.0588354 0.76046772 0.52419275 0.35331712 0.99235679] [0.98011228 0.12408712 0.95871245 0.11408017 0.11371645] [0.95161469 0.57466436 0.56620751 0.62610072 0.97239035] [0.63277575 0.64064988 0.74001459 0.24275853 0.43112123]] The inverse matrix, B [[ 1.34393676 -0.12950002 0.04036737 -0.50985099 0.23649582] [ 0.38507663 -0.6707882 -1.74480151 -0.68448882 3.20401293] [-1.15700156 0.33172141 1.25773672 0.19153362 -0.49344761] [-3.84061626 -3.55821886 -1.42800869 4.74830124 1.28910377] [ 1.60379604 2.62105895 1.17874303 -1.23698223 -2.66765296]] C = A * B, the product [[ 1.00000000e+00 -2.22044605e-16 -1.11022302e-16 5.55111512e-17 0.00000000e+00] [-2.22044605e-16 1.00000000e+00 -4.44089210e-16 4.44089210e-16 4.44089210e-16] [ 1.38777878e-16 -2.77555756e-16 1.00000000e+00 1.66533454e-16 2.22044605e-16] [ 2.22044605e-16 -4.44089210e-16 0.00000000e+00 1.00000000e+00 0.00000000e+00] [-2.22044605e-16 -8.88178420e-16 -4.44089210e-16 4.44089210e-16 1.00000000e+00]] R8GE_NP_ML_TEST R8GE_NP_ML computes A*x or A'*x for a matrix A already factored by R8GE_NP_FA. Matrix order N = 10 A*x and PLU*x 0: 26.8506 26.8506 1: 22.8737 22.8737 2: 18.2407 18.2407 3: 17.27 17.27 4: 28.6638 28.6638 5: 22.4529 22.4529 6: 24.47 24.47 7: 23.4776 23.4776 8: 22.4971 22.4971 9: 25.1852 25.1852 A'*x and (PLU)'*x 0: 25.939 25.939 1: 27.2595 27.2595 2: 31.3161 31.3161 3: 33.7721 33.7721 4: 20.694 20.694 5: 36.4207 36.4207 6: 25.5406 25.5406 7: 32.4011 32.4011 8: 24.0153 24.0153 9: 32.4613 32.4613 R8GE_NP_SL_TEST R8GE_NP_SL solves a linear system factored by R8GE_NP_FA. Matrix order N = 10 Solution: 0: 1 1: 1 2: 1 3: 1 4: 1 5: 1 6: 1 7: 1 8: 1 9: 1 Solution: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 8: 9 9: 10 Solution of transposed system: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 8: 9 9: 10 R8GE_NP_TRF_TEST R8GE_NP_TRF factors an R8GE matrix without pivoting, Matrix rows M = 10 Matrix columns N = 10 Solution: 0: 1 1: 1 2: 1 3: 1 4: 1 5: 1 6: 1 7: 1 8: 1 9: 1 Solution: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 8: 9 9: 10 Solution of transposed system: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 8: 9 9: 10 R8GE_NP_TRM_TEST R8GE_NP_TRM computes A*x after A has been factored by R8GE_NP_TRF. Matrix rows M = 10 Matrix columns N = 10 Solution: 0: 1 1: 1 2: 1 3: 1 4: 1 5: 1 6: 1 7: 1 8: 1 9: 1 Solution: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 8: 9 9: 10 Solution of transposed system: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 8: 9 9: 10 R8GE_NP_TRS_TEST R8GE_NP_TRS solves a linear system factored by R8GE_NP_TRF. Matrix rows M = 10 Matrix columns N = 10 Solution: 0: 1 1: 1 2: 1 3: 1 4: 1 5: 1 6: 1 7: 1 8: 1 9: 1 Solution: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 8: 9 9: 10 Solution of transposed system: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 8: 9 9: 10 r8ge_np_test(): Normal end of execution. Tue May 20 22:27:34 2025