26-Aug-2022 09:41:14 r8pp_test(): MATLAB version Test r8pp(). R8PP_DET_TEST R8PP_DET computes the determinant of an R8PP matrix which has been factored by R8PP_FA. Matrix order N = 5 The R8PP matrix: Col: 1 2 3 4 5 Row --- 1 2 -1 0 0 0 2 -1 2 -1 0 0 3 0 -1 2 -1 0 4 0 0 -1 2 -1 5 0 0 0 -1 2 Computed 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: 1 2 3 4 5 Row --- 1 2 -1 0 0 0 2 -1 2 -1 0 0 3 0 -1 2 -1 0 4 0 0 -1 2 -1 5 0 0 0 -1 2 R8PP_FA_TEST R8PP_FA factors an R8PP system, Matrix order N = 5 The R8PP matrix: 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 The desired solution: 1: 1 2: 2 3: 3 4: 4 5: 5 The right hand side: 1: 0.56603 2: 5.2015 3: 6.20596 4: 16.2555 5: 17.0626 The R8PP matrix has been factored. Solution: 1: 1 2: 2 3: 3 4: 4 5: 5 R8PP_INDICATOR_TEST R8PP_INDICATOR sets up a R8PP indicator matrix. Matrix order N = 5 The R8PP indicator 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_MV_TEST R8PP_MV computes b=A*x, where A is an R8PP matrix Matrix order N = 5 The R8PP indicator 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 Vector X: 1: 1 2: 2 3: 3 4: 4 5: 5 Product b=A*x 1: 205 2: 346 3: 469 4: 565 5: 625 R8PP_PRINT_TEST R8PP_PRINT prints an R8PP matrix. Matrix order N = 5 The R8PP 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_PRINT_SOME_TEST R8PP_PRINT_SOME prints some of an R8PP matrix. Matrix order N = 10 Rows 2-6, Cols 3-5: Col: 3 4 5 Row --- 2 203 204 205 3 303 304 305 4 304 404 405 5 305 405 505 6 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: 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 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: 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 The desired solution: 1: 1 2: 2 3: 3 4: 4 5: 5 The right hand side: 1: 0.56603 2: 5.2015 3: 6.20596 4: 16.2555 5: 17.0626 The R8PP matrix has been factored. Solution: 1: 1 2: 2 3: 3 4: 4 5: 5 R8PP_TO_R8GE_TEST R8PP_TO_R8GE converts an R8PP matrix to R8GE format. Matrix order N = 5 The R8PP indicator 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 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 sets an R8PP zero matrix. Matrix order N = 5 The R8PP zero matrix: Col: 1 2 3 4 5 Row --- 1 0 0 0 0 0 2 0 0 0 0 0 3 0 0 0 0 0 4 0 0 0 0 0 5 0 0 0 0 0 r8pp_test(): Normal end of execution. 26-Aug-2022 09:41:14