09 May 2025 9:29:15.366 PM r8po_test(): FORTRAN90 version: Test r8po(). R8GE_TO_R8PO_TEST R8GE_TO_R8PO converts an R8GE matrix to R8PO format. Matrix order N = 5 The random R8GE matrix A: Col: 1 2 3 4 5 Row --- 1 0.269507 0.767484 0.714255 0.794642E-01 0.580711 2 0.862970E-01 0.470868 0.991080 0.368289 0.190286 3 0.839386 0.674958 0.547155 0.630377 0.224927E-02 4 0.695687 0.575351 0.821229 0.460252E-01 0.380114 5 0.281711 0.923503 0.396213 0.811412 0.748773 The R8PO matrix Col: 1 2 3 4 5 Row --- 1 0.269507 0.767484 0.714255 0.794642E-01 0.580711 2 0.767484 0.470868 0.991080 0.368289 0.190286 3 0.714255 0.991080 0.547155 0.630377 0.224927E-02 4 0.794642E-01 0.368289 0.630377 0.460252E-01 0.380114 5 0.580711 0.190286 0.224927E-02 0.380114 0.748773 R8PO_DET_TEST R8PO_DET computes the determinant of a symmetric positive definite matrix factored by R8PO_FA,. Matrix order N = 4 Matrix A: Col: 1 2 3 4 Row --- 1 2.00000 -1.00000 0.00000 0.00000 2 -1.00000 2.00000 -1.00000 0.00000 3 0.00000 -1.00000 2.00000 -1.00000 4 0.00000 0.00000 -1.00000 2.00000 Determinant of A = 5.00000 R8PO_DIF2_TEST R8PO_DIF2 sets the second difference matrix in R8PO format. The R8PO DIF2 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 R8PO_FA_TEST R8PO_FA factors a positive definite symmetric linear system, Matrix order N = 5 The matrix A: Col: 1 2 3 4 5 Row --- 1 1.00000 1.00000 1.00000 1.00000 1.00000 2 1.00000 2.00000 2.00000 2.00000 2.00000 3 1.00000 2.00000 3.00000 3.00000 3.00000 4 1.00000 2.00000 3.00000 4.00000 4.00000 5 1.00000 2.00000 3.00000 4.00000 5.00000 The factor R (an R8UT matrix): Col: 1 2 3 4 5 Row --- 1 1.00000 1.00000 1.00000 1.00000 1.00000 2 1.00000 1.00000 1.00000 1.00000 3 1.00000 1.00000 1.00000 4 1.00000 1.00000 5 1.00000 The product R' * R: Col: 1 2 3 4 5 Row --- 1 1.00000 1.00000 1.00000 1.00000 1.00000 2 1.00000 2.00000 2.00000 2.00000 2.00000 3 1.00000 2.00000 3.00000 3.00000 3.00000 4 1.00000 2.00000 3.00000 4.00000 4.00000 5 1.00000 2.00000 3.00000 4.00000 5.00000 R8PO_INDICATOR_TEST R8PO_INDICATOR sets up an R8PO indicator matrix. Matrix order N = 5 The R8PO 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 R8PO_INVERSE_TEST R8PO_INVERSE computes the inverse of a symmetric positive definite matrix factored by R8PO_FA, Matrix order N = 4 Matrix A: Col: 1 2 3 4 Row --- 1 1.00000 1.00000 1.00000 1.00000 2 1.00000 2.00000 2.00000 2.00000 3 1.00000 2.00000 3.00000 3.00000 4 1.00000 2.00000 3.00000 4.00000 Inverse matrix B: Col: 1 2 3 4 Row --- 1 2.00000 -1.00000 0.00000 -0.00000 2 -1.00000 2.00000 -1.00000 -0.00000 3 0.00000 -1.00000 2.00000 -1.00000 4 -0.00000 -0.00000 -1.00000 1.00000 Product A * B: Col: 1 2 3 4 Row --- 1 1.00000 0.00000 0.00000 0.00000 2 0.00000 1.00000 0.00000 0.00000 3 0.00000 0.00000 1.00000 0.00000 4 0.00000 0.00000 0.00000 1.00000 R8PO_ML_TEST R8PO_ML can compute A*x for an R8PO matrix A even after it has been factored by R8PO_FA. Solution: 1: 1.0000000 2: 2.0000000 3: 3.0000000 4: 4.0000000 5: 5.0000000 6: 6.0000000 7: 7.0000000 8: 8.0000000 9: 9.0000000 10: 10.000000 Solution: 1: 1.0000000 2: 1.0000000 3: 1.0000000 4: 1.0000000 5: 1.0000000 6: 1.0000000 7: 1.0000000 8: 1.0000000 9: 1.0000000 10: 1.0000000 R8PO_MM_TEST R8PO_MM computes the product of two symmetric positive definite matrices. Matrix order N = 5 Matrix A: 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 Matrix B: Col: 1 2 3 4 5 Row --- 1 1.00000 2.00000 0.00000 0.00000 0.00000 2 2.00000 3.00000 4.00000 0.00000 0.00000 3 0.00000 4.00000 5.00000 6.00000 0.00000 4 0.00000 0.00000 6.00000 7.00000 8.00000 5 0.00000 0.00000 0.00000 8.00000 9.00000 Product matrix C = A * B: Col: 1 2 3 4 5 Row --- 1 0.00000 1.00000 -4.00000 0.00000 0.00000 2 1.00000 0.00000 3.00000 -6.00000 0.00000 3 -4.00000 3.00000 0.00000 5.00000 -8.00000 4 0.00000 -6.00000 5.00000 0.00000 7.00000 5 0.00000 0.00000 -8.00000 7.00000 10.0000 R8PO_MV_TEST R8PO_MV computes the product of a symmetric positive definite matrix times a vector. Matrix order N = 5 Matrix A: 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 Vector V: 1: 1.0000000 2: 2.0000000 3: 3.0000000 4: 4.0000000 5: 5.0000000 Product w = A * v: 1: 0.0000000 2: 0.0000000 3: 0.0000000 4: 0.0000000 5: 6.0000000 R8PO_PRINT_TEST R8PO_PRINT prints an R8PO matrix. The R8PO 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 R8PO_PRINT_SOME_TEST R8PO_PRINT_SOME prints some of an R8PO matrix. Rows 2-6, Columns 5-7: Col: 5 6 7 Row --- 2 205.000 206.000 207.000 3 305.000 306.000 307.000 4 405.000 406.000 407.000 5 505.000 506.000 507.000 6 506.000 606.000 607.000 R8PO_RANDOM_TEST R8PO_RANDOM computes a random positive definite symmetric matrix. Matrix order N = 5 The random R8PO matrix: Col: 1 2 3 4 5 Row --- 1 0.393430 0.606412 0.614513 0.398556 0.538201 2 0.606412 1.30235 1.14297 0.794837 0.930983 3 0.614513 1.14297 1.62834 0.769620 1.59812 4 0.398556 0.794837 0.769620 0.501092 0.685586 5 0.538201 0.930983 1.59812 0.685586 1.99733 R8PO_SL_TEST R8PO_SL solves a linear system with R8PO matrix after it has been factored buy R8PO_FA. Matrix order N = 5 Matrix A: 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 Right hand side b: 1: 0.0000000 2: 0.0000000 3: 0.0000000 4: 0.0000000 5: 6.0000000 Solution vector X: 1: 1.0000000 2: 2.0000000 3: 3.0000000 4: 4.0000000 5: 5.0000000 R8PO_TO_R8GE_TEST R8PO_TO_R8GE converts an R8PO matrix to R8GE format. Matrix order N = 5 The random R8PO matrix A: Col: 1 2 3 4 5 Row --- 1 0.225852E-01 0.142662 0.319902E-01 0.630728E-01 0.976954E-01 2 0.142662 1.42002 0.904708 0.891017 1.09937 3 0.319902E-01 0.904708 1.36453 1.10249 0.819027 4 0.630728E-01 0.891017 1.10249 1.35387 1.18782 5 0.976954E-01 1.09937 0.819027 1.18782 1.45938 The random R8PO matrix (printed by R8GE_PRINT): Col: 1 2 3 4 5 Row --- 1 0.225852E-01 0.142662 0.319902E-01 0.630728E-01 0.976954E-01 2 0.00000 1.42002 0.904708 0.891017 1.09937 3 0.00000 0.00000 1.36453 1.10249 0.819027 4 0.00000 0.00000 0.00000 1.35387 1.18782 5 0.00000 0.00000 0.00000 0.00000 1.45938 The random R8GE matrix (printed by R8GE_PRINT): Col: 1 2 3 4 5 Row --- 1 0.225852E-01 0.142662 0.319902E-01 0.630728E-01 0.976954E-01 2 0.142662 1.42002 0.904708 0.891017 1.09937 3 0.319902E-01 0.904708 1.36453 1.10249 0.819027 4 0.630728E-01 0.891017 1.10249 1.35387 1.18782 5 0.976954E-01 1.09937 0.819027 1.18782 1.45938 R8PO_ZEROS_TEST R8PO_ZEROS zeros out a matrix in R8PO format. Matrix order N = 5 The zero R8PO 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 r8po_test(): Normal end of execution. 09 May 2025 9:29:15.367 PM