08 June 2016 7:59:26.531 AM R8GE_NP_PRB FORTRAN90 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.0000 Exact value is 11.0000 R8GE_NP_FA_TEST R8GE_NP_FA computes the LU factors of a general storage matrix without pivoting, Matrix order N = 10 Solution: 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 10 1.000000 Solution: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 Solution of transposed system: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 R8GE_NP_INVERSE R8GE_NP_INVERSE computes the inverse of an R8GE matrix. Matrix order N = 5 The random matrix: Col: 1 2 3 4 5 Row --- 1 0.218418 0.661187E-01 0.617272E-01 0.183837E-02 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.438290E-01 0.754673 0.945448E-01 0.751236E-02 5 0.415307 0.633966 0.797287 0.136169E-01 0.260303 The inverse matrix: Col: 1 2 3 4 5 Row --- 1 0.625209 -1.02220 2.84978 -0.868794 -0.840984E-01 2 -0.471548 -0.149728 0.892032 -2.12070 1.67928 3 -0.274368 0.546937 -1.95265 2.05705 0.285593E-02 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.00000 -0.610623E-15 0.00000 -0.555112E-16 0.00000 2 0.976996E-14 1.00000 0.111022E-15 -0.277556E-16 -0.416334E-16 3 0.374700E-14 -0.915934E-15 1.00000 -0.173472E-16 0.607153E-16 4 0.364292E-14 -0.835920E-15 -0.404191E-15 1.00000 0.698226E-16 5 0.247580E-13 -0.420497E-14 -0.122125E-14 -0.617562E-15 1.00000 R8GE_NP_ML_TEST 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 1 32.6597 32.6597 2 30.7617 30.7617 3 20.1730 20.1730 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.5180 29.5180 2 34.2533 34.2533 3 29.9813 29.9813 4 25.9818 25.9818 5 28.7133 28.7133 6 28.3730 28.3730 7 33.4872 33.4872 8 26.6124 26.6124 9 24.8414 24.8414 10 26.5390 26.5390 R8GE_NP_SL_TEST R8GE_NP_SL solves a linear system that was factored by R8GE_NP_FA. Matrix order N = 10 Solution: 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 10 1.000000 Solution: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 Solution of transposed system: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 R8GE_NP_TRF_TEST R8GE_NP_TRF factors an R8GE matrix without pivoting, Matrix rows M = 10 Matrix columns N = 10 Solution: 1 1.00000 2 1.00000 3 1.00000 4 1.00000 5 1.00000 6 1.00000 7 1.00000 8 1.00000 9 1.00000 10 1.00000 Solution: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 Solution of transposed system: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 R8GE_NP_TRM_TEST R8GE_NP_TRM computes b=A*x after A has been factored by R8GE_NP_TRF. Matrix rows M = 10 Matrix columns N = 10 Solution: 1 1.00000 2 1.00000 3 1.00000 4 1.00000 5 1.00000 6 1.00000 7 1.00000 8 1.00000 9 1.00000 10 1.00000 Solution: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 Solution of transposed system: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 R8GE_NP_TRS_TEST R8GE_NP_TRS solves a linear system A*x=b which has been factored by R8GE_NP_TRF. Matrix rows M = 10 Matrix columns N = 10 Solution: 1 1.00000 2 1.00000 3 1.00000 4 1.00000 5 1.00000 6 1.00000 7 1.00000 8 1.00000 9 1.00000 10 1.00000 Solution: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 Solution of transposed system: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 6 6.00000 7 7.00000 8 8.00000 9 9.00000 10 10.0000 R8GE_NP_PRB Normal end of execution. 08 June 2016 7:59:26.532 AM