20 January 2017 10:12:11 AM LINPLUS_PRB C version Test the LINPLUS library. R83P_DET_TEST R83P_DET, determinant of a tridiagonal periodic matrix. Matrix order N = 12 The periodic tridiagonal matrix: Col: 1 2 3 4 5 Row --- 1 0.956318 0.561695 2 0.829509 0.415307 0.257578 3 0.066119 0.109957 0.633966 4 0.043829 0.061727 0.401306 5 0.449539 0.754673 6 0.797287 7 8 9 10 11 12 0.218418 Col: 6 7 8 9 10 Row --- 5 0.001838 6 0.897504 0.094545 7 0.350752 0.013617 0.840847 8 0.859097 0.123104 0.260303 9 0.007512 0.912484 0.351629 10 0.113664 0.822887 11 0.267132 Col: 11 12 Row --- 1 0.597917 2 3 4 5 6 7 8 9 10 0.692066 11 0.561662 0.453794 12 0.861216 0.911977 R83P_DET computes the determinant = -4.19241e-05 R8GE_DET computes the determinant = -4.19241e-05 R83P_FA_TEST R83P_FA factors a tridiagonal periodic system. R83P_SL solves a tridiagonal periodic system. Matrix order N = 10 Solution to A*x=b: 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 to A'*x=b: 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 R83P_INDICATOR_TEST R83P_INDICATOR sets up an R83P indicator matrix. Matrix order N = 5 The R83P indicator matrix: Col: 1 2 3 4 5 Row --- 1 11.000000 12.000000 15.000000 2 21.000000 22.000000 23.000000 3 32.000000 33.000000 34.000000 4 43.000000 44.000000 45.000000 5 51.000000 54.000000 55.000000 R83P_ML_TEST R83P_ML computes A*x or A'*X where A has been factored by R83P_FA. Matrix order N = 10 A*x and PLU*x 0: 4.751031 4.751031 1: 2.432857 2.432857 2: 2.997971 2.997971 3: 2.384927 2.384927 4: 5.582554 5.582554 5: 10.033272 10.033272 6: 8.926612 8.926612 7: 9.341236 9.341236 8: 11.788739 11.788739 9: 9.470268 9.470268 A'*x and (PLU)'*x 0: 9.204759 9.204759 1: 4.071498 4.071498 2: 4.250339 4.250339 3: 4.688539 4.688539 4: 5.554043 5.554043 5: 6.373901 6.373901 6: 11.217532 11.217532 7: 12.683676 12.683676 8: 17.734230 17.734230 9: 14.896107 14.896107 R83P_MTV_TEST R83P_MTV computes A'*x=b for an R83P matrix. Matrix order N = 5 The R83P indicator matrix: Col: 1 2 3 4 5 Row --- 1 11.000000 12.000000 15.000000 2 21.000000 22.000000 23.000000 3 32.000000 33.000000 34.000000 4 43.000000 44.000000 45.000000 5 51.000000 54.000000 55.000000 Vector x: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 b=A'*x: 0 308.000000 1 152.000000 2 317.000000 3 548.000000 4 470.000000 R83P_MV_TEST R83P_MV computes A*x=b for an R83P matrix. Matrix order N = 5 The R83P indicator matrix: Col: 1 2 3 4 5 Row --- 1 11.000000 12.000000 15.000000 2 21.000000 22.000000 23.000000 3 32.000000 33.000000 34.000000 4 43.000000 44.000000 45.000000 5 51.000000 54.000000 55.000000 Vector x: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 b=A*x: 0 110.000000 1 134.000000 2 299.000000 3 530.000000 4 542.000000 R83P_PRINT_TEST R83P_PRINT prints an R83P matrix. Matrix order N = 5 The R83P matrix: Col: 1 2 3 4 5 Row --- 1 11.000000 12.000000 15.000000 2 21.000000 22.000000 23.000000 3 32.000000 33.000000 34.000000 4 43.000000 44.000000 45.000000 5 51.000000 54.000000 55.000000 R83P_PRINT_SOME_TEST R83P_PRINT_SOME prints some of an R83P matrix. Matrix order N = 10 Rows 1:N, Cols 1:2: Col: 1 2 Row --- 1 101.000000 102.000000 2 201.000000 202.000000 3 302.000000 4 5 6 7 8 9 10 1001.000000 R83P_RANDOM_TEST R83P_RANDOM sets up a random R83P matrix. Matrix order N = 5 The random R83P matrix: Col: 1 2 3 4 5 Row --- 1 0.956318 0.561695 0.797287 2 0.829509 0.415307 0.257578 3 0.066119 0.109957 0.633966 4 0.043829 0.061727 0.401306 5 0.218418 0.449539 0.754673 R83P_SL_TEST R83P_SL solves a tridiagonal periodic system after it has been factored by R83P_FA. Matrix order N = 10 Solution to A*x=b: 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 to A'*x=b: 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 R83P_TO_R8GE_TEST R83P_TO_R8GE convert a matrix from R83P to R8GE format. Matrix order N = 5 The R83P matrix: Col: 1 2 3 4 5 Row --- 1 0.956318 0.561695 0.797287 2 0.829509 0.415307 0.257578 3 0.066119 0.109957 0.633966 4 0.043829 0.061727 0.401306 5 0.218418 0.449539 0.754673 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 0.956318 0.561695 0 0 0.797287 2 0.829509 0.415307 0.257578 0 0 3 0 0.0661187 0.109957 0.633966 0 4 0 0 0.043829 0.0617272 0.401306 5 0.218418 0 0 0.449539 0.754673 R83P_ZEROS_TEST R83P_ZEROS sets up a R83P zero matrix. Matrix order N = 5 The R83P zero matrix: Col: 1 2 3 4 5 Row --- 1 0.000000 0.000000 0.000000 2 0.000000 0.000000 0.000000 3 0.000000 0.000000 0.000000 4 0.000000 0.000000 0.000000 5 0.000000 0.000000 0.000000 LINPLUS_PRB Normal end of execution. 20 January 2017 10:12:11 AM