06 October 2025 6:33:36.843 PM r8pp_test(): Fortran90 version: Test r8pp(). R8PP_DET_TEST R8PP_DET computes the determinant of an R8PP matrix factored by R8PP_FA. Matrix order N = 5 The R8PP matrix: Col: 1 2 3 4 5 Row --- 1 2.00000 -1.00000 0.00000 0.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 0.00000 4 0.00000 0.00000 -1.00000 2.00000 -1.00000 5 0.00000 0.00000 0.00000 -1.00000 2.00000 Determinant = 6.00000 Exact determinant = 6.00000 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.00000 -1.00000 0.00000 0.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 0.00000 4 0.00000 0.00000 -1.00000 2.00000 -1.00000 5 0.00000 0.00000 0.00000 -1.00000 2.00000 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.553362E-02 0.162003E-01 0.735886E-01 0.663626E-01 0.107801E-01 2 0.162003E-01 0.609863 0.596068 0.221301 0.603582 3 0.735886E-01 0.596068 1.34944 1.02412 0.617155 4 0.663626E-01 0.221301 1.02412 0.981978 0.370499 5 0.107801E-01 0.603582 0.617155 0.370499 1.82873 The desired solution: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 The right hand side: 1 0.578051 2 6.92725 3 12.4963 4 9.36175 5 13.6951 The R8PP matrix has been factored. Solution: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 R8PP_INDICATOR_TEST R8PP_INDICATOR sets up an R8PP indicator matrix. Matrix order N = 5 The R8PP indicator matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 13.0000 14.0000 15.0000 2 12.0000 22.0000 23.0000 24.0000 25.0000 3 13.0000 23.0000 33.0000 34.0000 35.0000 4 14.0000 24.0000 34.0000 44.0000 45.0000 5 15.0000 25.0000 35.0000 45.0000 55.0000 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.0000 12.0000 13.0000 14.0000 15.0000 2 12.0000 22.0000 23.0000 24.0000 25.0000 3 13.0000 23.0000 33.0000 34.0000 35.0000 4 14.0000 24.0000 34.0000 44.0000 45.0000 5 15.0000 25.0000 35.0000 45.0000 55.0000 Vector x: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 Result b=A*x: 1 205.000 2 346.000 3 469.000 4 565.000 5 625.000 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.0000 12.0000 13.0000 14.0000 15.0000 2 12.0000 22.0000 23.0000 24.0000 25.0000 3 13.0000 23.0000 33.0000 34.0000 35.0000 4 14.0000 24.0000 34.0000 44.0000 45.0000 5 15.0000 25.0000 35.0000 45.0000 55.0000 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.000 204.000 205.000 3 303.000 304.000 305.000 4 304.000 404.000 405.000 5 305.000 405.000 505.000 6 306.000 406.000 506.000 R8PP_RANDOM_TEST R8PP_RANDOM, compute a random positive definite symmetric packed matrix. Matrix order N = 5 The random R8PP matrix: Col: 1 2 3 4 5 Row --- 1 0.162036 0.193289 0.151943 0.397190 0.392372 2 0.193289 1.04390 0.549082 0.740521 0.960574 3 0.151943 0.549082 1.09930 1.26361 1.45871 4 0.397190 0.740521 1.26361 2.23188 2.00570 5 0.392372 0.960574 1.45871 2.00570 2.23392 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.850978 0.501587 0.226494 0.577121 0.853680 2 0.501587 1.12502 0.795400 0.367061 0.648582 3 0.226494 0.795400 1.41958 0.921149 0.660964 4 0.577121 0.367061 0.921149 1.98380 1.37450 5 0.853680 0.648582 0.660964 1.37450 1.29428 The desired solution: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 The right hand side: 1 9.11052 2 9.84897 3 13.0655 4 18.8824 5 16.1031 The R8PP matrix has been factored. Solution: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 R8PP_TO_R8GE_TEST R8PP_TO_R8GE converts an R8PP matrix to R8GE format. Matrix order N = 5 The R8PP matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 13.0000 14.0000 15.0000 2 12.0000 22.0000 23.0000 24.0000 25.0000 3 13.0000 23.0000 33.0000 34.0000 35.0000 4 14.0000 24.0000 34.0000 44.0000 45.0000 5 15.0000 25.0000 35.0000 45.0000 55.0000 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 11.0000 12.0000 13.0000 14.0000 15.0000 2 12.0000 22.0000 23.0000 24.0000 25.0000 3 13.0000 23.0000 33.0000 34.0000 35.0000 4 14.0000 24.0000 34.0000 44.0000 45.0000 5 15.0000 25.0000 35.0000 45.0000 55.0000 R8PP_ZEROS_TEST R8PP_ZEROS sets a zero R8PP matrix. Matrix order N = 5 The R8PP matrix: Col: 1 2 3 4 5 Row --- 1 0.00000 0.00000 0.00000 0.00000 0.00000 2 0.00000 0.00000 0.00000 0.00000 0.00000 3 0.00000 0.00000 0.00000 0.00000 0.00000 4 0.00000 0.00000 0.00000 0.00000 0.00000 5 0.00000 0.00000 0.00000 0.00000 0.00000 r8pp_test(): Normal end of execution. 06 October 2025 6:33:36.844 PM