20 January 2017 11:31:36 AM R8PP_PRB C version Test the R8PP library. R8PP_DET_TEST R8PP_DET computes the determinant of an R8PP matrix factored by R8PP_FA. Matrix order N = 5 The R8PP matrix: Col: 0 1 2 3 4 Row --- 0 2 -1 0 0 0 1 -1 2 -1 0 0 2 0 -1 2 -1 0 3 0 0 -1 2 -1 4 0 0 0 -1 2 Determinant = 6 Exact determinant = 6 R8PP_DIF2_TEST R8PP_DIF2 sets up an R8PP second difference matrix. Matrix order N = 5 The R8PP second difference matrix: Col: 0 1 2 3 4 Row --- 0 2 -1 0 0 0 1 -1 2 -1 0 0 2 0 -1 2 -1 0 3 0 0 -1 2 -1 4 0 0 0 -1 2 R8PP_FA_TEST R8PP_FA factors an R8PP system, Matrix order N = 5 The R8PP matrix: Col: 0 1 2 3 4 Row --- 0 0.00381025 0.0277488 0.0247715 0.0465839 0.0492143 1 0.0277488 0.268432 0.208725 0.350545 0.521707 2 0.0247715 0.208725 0.488639 0.540951 0.426804 3 0.0465839 0.350545 0.540951 1.65848 1.45021 4 0.0492143 0.521707 0.426804 1.45021 1.77774 The desired solution: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 The right hand side: 0 0.566030 1 5.201501 2 6.205960 3 16.255486 4 17.062599 The R8PP matrix has been factored. Solution: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 R8PP_INDICATOR_TEST R8PP_INDICATOR sets up an R8PP indicator matrix. Matrix order N = 5 The R8PP indicator matrix: Col: 0 1 2 3 4 Row --- 0 11 12 13 14 15 1 12 22 23 24 25 2 13 23 33 34 35 3 14 24 34 44 45 4 15 25 35 45 55 R8PP_MV_TEST R8PP_MV computes b=A*x, where A is an R8PP matrix. Matrix order N = 5 The R8PP matrix A: Col: 0 1 2 3 4 Row --- 0 11 12 13 14 15 1 12 22 23 24 25 2 13 23 33 34 35 3 14 24 34 44 45 4 15 25 35 45 55 Vector x: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 Product b=A*x: 0 205.000000 1 346.000000 2 469.000000 3 565.000000 4 625.000000 R8PP_PRINT_TEST R8PP_PRINT prints an R8PP matrix. Matrix order N = 5 The R8PP matrix: Col: 0 1 2 3 4 Row --- 0 2 -1 0 0 0 1 -1 2 -1 0 0 2 0 -1 2 -1 0 3 0 0 -1 2 -1 4 0 0 0 -1 2 R8PP_PRINT_SOME_TEST R8PP_PRINT_SOME prints some of an R8PP matrix. Matrix order N = 10 Rows 1-5, Cols 2-4: Col: 2 3 4 Row --- 1 203 204 205 2 303 304 305 3 304 404 405 4 305 405 505 5 306 406 506 R8PP_RANDOM_TEST R8PP_RANDOM, compute a random positive definite symmetric packed matrix. Matrix order N = 5 The matrix (printed by R8PP_PRINT): Col: 0 1 2 3 4 Row --- 0 0.00381025 0.0277488 0.0247715 0.0465839 0.0492143 1 0.0277488 0.268432 0.208725 0.350545 0.521707 2 0.0247715 0.208725 0.488639 0.540951 0.426804 3 0.0465839 0.350545 0.540951 1.65848 1.45021 4 0.0492143 0.521707 0.426804 1.45021 1.77774 The random R8PP matrix (printed by R8GE_PRINT): Col: 1 2 3 4 5 Row --- 1 0.00381025 0.0277488 0.0247715 0.0465839 0.0492143 2 0.0277488 0.268432 0.208725 0.350545 0.521707 3 0.0247715 0.208725 0.488639 0.540951 0.426804 4 0.0465839 0.350545 0.540951 1.65848 1.45021 5 0.0492143 0.521707 0.426804 1.45021 1.77774 R8PP_SL_TEST R8PP_SL solves a linear system factored by R8PP_FA. Matrix order N = 5 The R8PP matrix: Col: 0 1 2 3 4 Row --- 0 0.00381025 0.0277488 0.0247715 0.0465839 0.0492143 1 0.0277488 0.268432 0.208725 0.350545 0.521707 2 0.0247715 0.208725 0.488639 0.540951 0.426804 3 0.0465839 0.350545 0.540951 1.65848 1.45021 4 0.0492143 0.521707 0.426804 1.45021 1.77774 The desired solution: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 The right hand side: 0 0.566030 1 5.201501 2 6.205960 3 16.255486 4 17.062599 The R8PP matrix has been factored. Solution: 0 1.000000 1 2.000000 2 3.000000 3 4.000000 4 5.000000 R8PP_TO_R8GE_TEST R8PP_TO_R8GE converts an R8PP matrix to R8GE format. Matrix order N = 5 The R8PP indicator matrix: Col: 0 1 2 3 4 Row --- 0 11 12 13 14 15 1 12 22 23 24 25 2 13 23 33 34 35 3 14 24 34 44 45 4 15 25 35 45 55 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 11 12 13 14 15 2 12 22 23 24 25 3 13 23 33 34 35 4 14 24 34 44 45 5 15 25 35 45 55 R8PP_ZEROS_TEST R8PP_ZEROS zeros an R8PP matrix. Matrix order N = 5 The R8PP zero matrix: Col: 0 1 2 3 4 Row --- 0 0 0 0 0 0 1 0 0 0 0 0 2 0 0 0 0 0 3 0 0 0 0 0 4 0 0 0 0 0 R8PP_PRB Normal end of execution. 20 January 2017 11:31:36 AM