20 January 2017 11:25:39 AM R8GE_NP_PRB C version Test the R8GE_NP library. 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: 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.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 5 6.000000 6 7.000000 7 8.000000 8 9.000000 9 10.000000 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 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 For a matrix in general storage, R8GE_NP_ML computes A*x or A'*X where A has been factored by R8GE_NP_FA. Matrix order N = 10 A*x and PLU*x 0: 32.659651 32.659651 1: 30.761657 30.761657 2: 20.172970 20.172970 3: 24.112617 24.112617 4: 30.759727 30.759727 5: 32.734219 32.734219 6: 28.025246 28.025246 7: 28.309083 28.309083 8: 28.750123 28.750123 9: 17.314086 17.314086 A'*x and (PLU)'*x 0: 29.518040 29.518040 1: 34.253309 34.253309 2: 29.981306 29.981306 3: 25.981850 25.981850 4: 28.713293 28.713293 5: 28.372974 28.372974 6: 33.487236 33.487236 7: 26.612402 26.612402 8: 24.841412 24.841412 9: 26.538970 26.538970 R8GE_NP_SL_TEST R8GE_NP_SL solves a 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.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 5 6.000000 6 7.000000 7 8.000000 8 9.000000 9 10.000000 R8GE_NP_TRF_TEST R8GE_NP_TRF factors an R8GE matrix without pivoting, Matrix rows M = 10 Matrix columns N = 10 Solution: 0 1.000000 1 1.000000 2 1.000000 3 1.000000 4 1.000000 5 1.000000 6 1.000000 7 1.000000 8 1.000000 9 1.000000 Solution: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 5 6.000000 6 7.000000 7 8.000000 8 9.000000 9 10.000000 Solution of transposed system: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 5 6.000000 6 7.000000 7 8.000000 8 9.000000 9 10.000000 R8GE_NP_TRM_TEST R8GE_NP_TRM computes A*x when A was factored by R8GE_NP_TRF. Matrix rows M = 10 Matrix columns N = 10 Solution: 0 1.000000 1 1.000000 2 1.000000 3 1.000000 4 1.000000 5 1.000000 6 1.000000 7 1.000000 8 1.000000 9 1.000000 Solution: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 5 6.000000 6 7.000000 7 8.000000 8 9.000000 9 10.000000 Solution of transposed system: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 5 6.000000 6 7.000000 7 8.000000 8 9.000000 9 10.000000 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.000000 1 1.000000 2 1.000000 3 1.000000 4 1.000000 5 1.000000 6 1.000000 7 1.000000 8 1.000000 9 1.000000 Solution: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 5 6.000000 6 7.000000 7 8.000000 8 9.000000 9 10.000000 Solution of transposed system: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 5 6.000000 6 7.000000 7 8.000000 8 9.000000 9 10.000000 R8GE_NP_PRB Normal end of execution. 20 January 2017 11:25:39 AM