10 April 2023 4:18:30.176 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.381025E-02 0.277488E-01 0.247715E-01 0.465839E-01 0.492143E-01 2 0.277488E-01 0.268432 0.208725 0.350545 0.521707 3 0.247715E-01 0.208725 0.488639 0.540951 0.426804 4 0.465839E-01 0.350545 0.540951 1.65848 1.45021 5 0.492143E-01 0.521707 0.426804 1.45021 1.77774 The desired solution: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 The right hand side: 1 0.566030 2 5.20150 3 6.20596 4 16.2555 5 17.0626 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.381025E-02 0.277488E-01 0.247715E-01 0.465839E-01 0.492143E-01 2 0.277488E-01 0.268432 0.208725 0.350545 0.521707 3 0.247715E-01 0.208725 0.488639 0.540951 0.426804 4 0.465839E-01 0.350545 0.540951 1.65848 1.45021 5 0.492143E-01 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.381025E-02 0.277488E-01 0.247715E-01 0.465839E-01 0.492143E-01 2 0.277488E-01 0.268432 0.208725 0.350545 0.521707 3 0.247715E-01 0.208725 0.488639 0.540951 0.426804 4 0.465839E-01 0.350545 0.540951 1.65848 1.45021 5 0.492143E-01 0.521707 0.426804 1.45021 1.77774 The desired solution: 1 1.00000 2 2.00000 3 3.00000 4 4.00000 5 5.00000 The right hand side: 1 0.566030 2 5.20150 3 6.20596 4 16.2555 5 17.0626 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. 10 April 2023 4:18:30.177 PM