23-Aug-2022 07:01:53 r8ge_np_test(): MATLAB version 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 11 Exact value is 11 R8GE_NP_FA_TEST R8GE_NP_FA LU factors an R8GE matrix without pivoting. Matrix order N = 10 Solution: 1: 1 2: 1 3: 1 4: 1 5: 1 6: 1 7: 1 8: 1 9: 1 10: 1 Solution: 1: 1 2: 2 3: 3 4: 4 5: 5 6: 6 7: 7 8: 8 9: 9 10: 10 Solution of transposed system: 1: 1 2: 2 3: 3 4: 4 5: 5 6: 6 7: 7 8: 8 9: 9 10: 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 random matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.0661187 0.0617272 0.00183837 0.859097 2 0.956318 0.257578 0.449539 0.897504 0.840847 3 0.829509 0.109957 0.401306 0.350752 0.123104 4 0.561695 0.043829 0.754673 0.0945448 0.00751236 5 0.415307 0.633966 0.797287 0.0136169 0.260303 B so far Col: 1 2 3 4 5 Row --- 1 4.57837 9.48509 3.16785 -1.05328 -0.0840984 2 -4.37838 -31.3333 -8.97338 1.56317 1.67928 3 -3.7978 -4.42266 -1.59747 2.06331 0.00285593 4 -2.57165 -3.95443 -0.180503 -0.359782 -0.293107 5 -1.90143 15.925 5.64683 -2.19371 -0.107439 The inverse matrix: Col: 1 2 3 4 5 Row --- 1 0.625209 -1.0222 2.84978 -0.868794 -0.0840984 2 -0.471548 -0.149728 0.892032 -2.1207 1.67928 3 -0.274368 0.546937 -1.95265 2.05705 0.00285593 4 -1.39029 1.75848 -1.70625 0.283212 -0.293107 5 1.06404 0.228349 -0.649235 0.235691 -0.107439 The product: Col: 1 2 3 4 5 Row --- 1 1 -6.10623e-16 0 -5.55112e-17 0 2 9.76996e-15 1 1.11022e-16 -2.77556e-17 -4.16334e-17 3 3.747e-15 -9.15934e-16 1 -1.73472e-17 6.07153e-17 4 3.64292e-15 -8.3592e-16 -4.04191e-16 1 6.98226e-17 5 2.4758e-14 -4.20497e-15 -1.22125e-15 -6.17562e-16 1 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 1: 32.6597 32.6597 2: 30.7617 30.7617 3: 20.173 20.173 4: 24.1126 24.1126 5: 30.7597 30.7597 6: 32.7342 32.7342 7: 28.0252 28.0252 8: 28.3091 28.3091 9: 28.7501 28.7501 10: 17.3141 17.3141 A'*x and (PLU)'*x 1: 29.518 29.518 2: 34.2533 34.2533 3: 29.9813 29.9813 4: 25.9818 25.9818 5: 28.7133 28.7133 6: 28.373 28.373 7: 33.4872 33.4872 8: 26.6124 26.6124 9: 24.8414 24.8414 10: 26.539 26.539 R8GE_NP_SL_TEST R8GE_NP_SL solves a linear system factored by R8GE_NP_FA. Matrix order N = 10 Solution: 1 1 2 1 3 1 4 1 5 1 6 1 7 1 8 1 9 1 10 1 Solution: 1: 1 2: 2 3: 3 4: 4 5: 5 6: 6 7: 7 8: 8 9: 9 10: 10 Solution of transposed system: 1: 1 2: 2 3: 3 4: 4 5: 5 6: 6 7: 7 8: 8 9: 9 10: 10 R8GE_NP_TRF_TEST R8GE_NP_TRF factors an R8GE matrix without pivoting, Matrix rows M = 10 Matrix columns N = 10 Solution: 1: 1 2: 1 3: 1 4: 1 5: 1 6: 1 7: 1 8: 1 9: 1 10: 1 Solution: 1: 1 2: 2 3: 3 4: 4 5: 5 6: 6 7: 7 8: 8 9: 9 10: 10 Solution of transposed system: 1: 1 2: 2 3: 3 4: 4 5: 5 6: 6 7: 7 8: 8 9: 9 10: 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: 1: 1 2: 1 3: 1 4: 1 5: 1 6: 1 7: 1 8: 1 9: 1 10: 1 Solution: 1: 1 2: 2 3: 3 4: 4 5: 5 6: 6 7: 7 8: 8 9: 9 10: 10 Solution of transposed system: 1: 1 2: 2 3: 3 4: 4 5: 5 6: 6 7: 7 8: 8 9: 9 10: 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: 1: 1 2: 1 3: 1 4: 1 5: 1 6: 1 7: 1 8: 1 9: 1 10: 1 Solution: 1: 1 2: 2 3: 3 4: 4 5: 5 6: 6 7: 7 8: 8 9: 9 10: 10 Solution of transposed system: 1: 1 2: 2 3: 3 4: 4 5: 5 6: 6 7: 7 8: 8 9: 9 10: 10 r8ge_np_test(): Normal end of execution. 23-Aug-2022 07:01:53