18 May 2016 08:38:28 PM R83P_PRB C++ version Test the R83P 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.0661187 0.109957 0.633966 4 0.043829 0.0617272 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.00183837 6 0.897504 0.0945448 7 0.350752 0.0136169 0.840847 8 0.859097 0.123104 0.260303 9 0.00751236 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. which can then be solved by R83P_SL. Matrix order N = 10 Solution to A*x=b: 1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 10 10 Solution to A'*x=b: 1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 10 10 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 12 15 2 21 22 23 3 32 33 34 4 43 44 45 5 51 54 55 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 1 4.75103 4.75103 2 2.43286 2.43286 3 2.99797 2.99797 4 2.38493 2.38493 5 5.58255 5.58255 6 10.0333 10.0333 7 8.92661 8.92661 8 9.34124 9.34124 9 11.7887 11.7887 10 9.47027 9.47027 A'*x and (PLU)'*x 1 9.20476 9.20476 2 4.0715 4.0715 3 4.25034 4.25034 4 4.68854 4.68854 5 5.55404 5.55404 6 6.3739 6.3739 7 11.2175 11.2175 8 12.6837 12.6837 9 17.7342 17.7342 10 14.8961 14.8961 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 12 15 2 21 22 23 3 32 33 34 4 43 44 45 5 51 54 55 Vector x: 1 1 2 2 3 3 4 4 5 5 b=A'*x: 1 308 2 152 3 317 4 548 5 470 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 12 15 2 21 22 23 3 32 33 34 4 43 44 45 5 51 54 55 Vector x: 1 1 2 2 3 3 4 4 5 5 b=A*x: 1 110 2 134 3 299 4 530 5 542 R83P_PRINT_TEST R83P_PRINT prints an R83P matrix. Matrix order N = 5 The R83P indicator matrix: Col: 1 2 3 4 5 Row --- 1 11 12 15 2 21 22 23 3 32 33 34 4 43 44 45 5 51 54 55 R83P_PRINT_SOME_TEST R83P_PRINT_SOME prints some of an R83P matrix. Matrix order N = 5 Row(1:N),Col(1:2) Col: 1 2 Row --- 1 11 12 2 21 22 3 32 4 5 51 R83P_RANDOM_TEST R83P_RANDOM returns a random R83P matrix. 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.0661187 0.109957 0.633966 4 0.043829 0.0617272 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: 1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 10 10 Solution to A'*x=b: 1 1 2 2 3 3 4 4 5 5 6 6 7 7 8 8 9 9 10 10 R83P_TO_R8GE_TEST R83P_R8GE converts a matrix from R83P to R8GE format. Matrix order N = 5 The R83P matrix: Col: 1 2 3 4 5 Row --- 1 11 12 15 2 21 22 23 3 32 33 34 4 43 44 45 5 51 54 55 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 11 12 0 0 15 2 21 22 23 0 0 3 0 32 33 34 0 4 0 0 43 44 45 5 51 0 0 54 55 R83P_ZEROS_TEST R83P_ZEROS returns a zero R83P matrix. Matrix order N = 5 The R83P matrix: Col: 1 2 3 4 5 Row --- 1 0 0 0 2 0 0 0 3 0 0 0 4 0 0 0 5 0 0 0 R83P_PRB Normal end of execution. 18 May 2016 08:38:28 PM