Tue Aug 23 07:39:31 2022 r8ge_np_test(): Python version: 3.6.9 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.52443104 0.55218076 0.02964901 0.54398192 0.39568472] [0.11695021 0.31923684 0.09350822 0.25800002 0.68911447] [0.78863635 0.00718747 0.30313117 0.25656601 0.88139378] [0.80017965 0.3848524 0.67004309 0.01628026 0.49654146] [0.78502001 0.73645826 0.9001442 0.2577137 0.60054704]] The inverse matrix, B [[ 1.08516036 -0.84945564 0.02197216 2.21246083 -1.60179602] [ 1.00737052 1.66767704 -2.31727464 3.02828636 -1.6802344 ] [-1.47916706 -1.20503672 1.19019203 -3.31720069 3.35326251] [ 0.24408757 -2.09654142 2.20633997 -5.57935477 3.61986959] [-0.5415063 1.77118415 0.08222333 0.76064412 -0.76004928]] C = A * B, the product [[ 1.00000000e+00 -1.33226763e-15 5.68989300e-16 -9.43689571e-16 4.44089210e-16] [ 7.77156117e-16 1.00000000e+00 -2.91433544e-16 -1.11022302e-16 1.11022302e-16] [ 7.77156117e-16 -8.88178420e-16 1.00000000e+00 -9.99200722e-16 5.55111512e-16] [ 1.66533454e-15 -2.44249065e-15 2.15105711e-16 1.00000000e+00 3.33066907e-16] [ 2.66453526e-15 -2.44249065e-15 -1.04083409e-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: 24.2406 24.2406 1: 12.9425 12.9425 2: 31.285 31.285 3: 27.901 27.901 4: 23.7398 23.7398 5: 30.245 30.245 6: 24.1935 24.1935 7: 26.0867 26.0867 8: 27.1334 27.1334 9: 28.8117 28.8117 A'*x and (PLU)'*x 0: 23.782 23.782 1: 35.556 35.556 2: 24.7938 24.7938 3: 26.9103 26.9103 4: 39.8748 39.8748 5: 27.5538 27.5538 6: 31.8767 31.8767 7: 26.6884 26.6884 8: 30.7617 30.7617 9: 30.8024 30.8024 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 Aug 23 07:39:31 2022