09 May 2025 9:29:15.497 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.548698E-02 0.570232E-01 0.102330E-01 0.305607E-01 0.823952E-02 2 0.570232E-01 0.790178 0.203050 0.482627 0.209931 3 0.102330E-01 0.203050 0.299671 0.179940 0.187123 4 0.305607E-01 0.482627 0.179940 0.413271 0.479699 5 0.823952E-02 0.209931 0.187123 0.479699 1.13202 The desired solution: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 The right hand side: 1 0.313673 2 5.22669 3 2.97072 4 5.58721 5 8.56838 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.453918E-01 0.503482E-01 0.116001 0.171690 0.103561 2 0.503482E-01 0.903138E-01 0.141531 0.296098 0.202715 3 0.116001 0.141531 0.440306 0.796829 0.543264 4 0.171690 0.296098 0.796829 1.71088 1.23164 5 0.103561 0.202715 0.543264 1.23164 1.68839 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.497852 0.585265 0.103975 0.405252E-01 0.256849 2 0.585265 0.707469 0.137090 0.104150 0.430984 3 0.103975 0.137090 0.131890 0.272216 0.419319 4 0.405252E-01 0.104150 0.272216 1.16504 1.55018 5 0.256849 0.430984 0.419319 1.55018 2.68798 The desired solution: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 The right hand side: 1 3.42665 2 4.98299 3 3.95928 4 13.4765 5 22.0174 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. 09 May 2025 9:29:15.498 PM