Fri Aug 12 23:49:05 2022 r8pp_test(): Python version: 3.6.9 Test r8pp(). r8ge_to_r8pp_test(): r8ge_to_r8pp() converts an R8GE matrix to R8PP format. Matrix order N = 5 The positive definite symmetric R8GE matrix: [[-1. -1. -1. -1. -1.] [-1. 1. 1. 1. 1.] [-1. 1. 3. 3. 3.] [-1. 1. 3. 5. 5.] [-1. 1. 3. 5. 7.]] The RPP matrix: [-1. -1. 1. -1. 1. 3. -1. 1. 3. 5. -1. 1. 3. 5. 7.] R8PP_DET_TEST R8PP_DET computes the determinant of an R8PP matrix. Matrix order N = 5 The R8PP matrix: Col: 0 1 2 3 4 Row 0 : 2 -1 0 0 0 1 : -1 2 -1 0 0 2 : 0 -1 2 -1 0 3 : 0 0 -1 2 -1 4 : 0 0 0 -1 2 Computed determinant = 5.999999999999999 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: 0 1 2 3 4 Row 0 : 2 -1 0 0 0 1 : -1 2 -1 0 0 2 : 0 -1 2 -1 0 3 : 0 0 -1 2 -1 4 : 0 0 0 -1 2 R8PP_FA_TEST R8PP_FA factors an R8PP system, Matrix order N = 5 The R8PP matrix: Col: 0 1 2 3 4 Row 0 : 0.00381025 0.0277488 0.0247715 0.0465839 0.0492143 1 : 0.0277488 0.268432 0.208725 0.350545 0.521707 2 : 0.0247715 0.208725 0.488639 0.540951 0.426804 3 : 0.0465839 0.350545 0.540951 1.65848 1.45021 4 : 0.0492143 0.521707 0.426804 1.45021 1.77774 The desired solution: 0: 1 1: 2 2: 3 3: 4 4: 5 The right hand side: 0: 0.56603 1: 5.2015 2: 6.20596 3: 16.2555 4: 17.0626 The R8PP matrix has been factored. Solution: 0: 1 1: 2 2: 3 3: 4 4: 5 R8PP_INDICATOR_TEST R8PP_INDICATOR sets up a R8PP indicator matrix. Matrix order N = 5 The R8PP indicator matrix: Col: 0 1 2 3 4 Row 0 : 11 12 13 14 15 1 : 12 22 23 24 25 2 : 13 23 33 34 35 3 : 14 24 34 44 45 4 : 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: 0 1 2 3 4 Row 0 : 11 12 13 14 15 1 : 12 22 23 24 25 2 : 13 23 33 34 35 3 : 14 24 34 44 45 4 : 15 25 35 45 55 Vector X: 0: 1 1: 2 2: 3 3: 4 4: 5 Product b=A*x 0: 205 1: 346 2: 469 3: 565 4: 625 r8pp_print_test(): r8pp_print() prints an R8PP matrix. Matrix order N = 5 The R8PP matrix: Col: 0 1 2 3 4 Row 0 : 11 12 13 14 15 1 : 12 22 23 24 25 2 : 13 23 33 34 35 3 : 14 24 34 44 45 4 : 15 25 35 45 55 r8pp_print_some_test(): r8pp_print_some() prints some of an R8PP matrix. Matrix order N = 10 Rows 1-5, Cols 2-4: Col: 2 3 4 Row 1 : 203 204 205 2 : 303 304 305 3 : 304 404 405 4 : 305 405 505 5 : 306 406 506 R8PP_RANDOM_TEST R8PP_RANDOM, compute a random symmetric positive definite packed matrix. Matrix order N = 5 The random R8PP matrix Col: 0 1 2 3 4 Row 0 : 0.00381025 0.0277488 0.0247715 0.0465839 0.0492143 1 : 0.0277488 0.268432 0.208725 0.350545 0.521707 2 : 0.0247715 0.208725 0.488639 0.540951 0.426804 3 : 0.0465839 0.350545 0.540951 1.65848 1.45021 4 : 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: 0 1 2 3 4 Row 0 : 0.00381025 0.0277488 0.0247715 0.0465839 0.0492143 1 : 0.0277488 0.268432 0.208725 0.350545 0.521707 2 : 0.0247715 0.208725 0.488639 0.540951 0.426804 3 : 0.0465839 0.350545 0.540951 1.65848 1.45021 4 : 0.0492143 0.521707 0.426804 1.45021 1.77774 The desired solution: 0: 1 1: 2 2: 3 3: 4 4: 5 The right hand side: 0: 0.56603 1: 5.2015 2: 6.20596 3: 16.2555 4: 17.0626 The R8PP matrix has been factored. Solution: 0: 1 1: 2 2: 3 3: 4 4: 5 r8pp_to_r8ge_test() r8pp_to_r8ge() converts an R8PP matrix to R8GE format. Matrix order N = 5 The R8PP indicator matrix: [11. 12. 22. 13. 23. 33. 14. 24. 34. 44. 15. 25. 35. 45. 55.] The R8GE matrix: [[11. 12. 13. 14. 15.] [12. 22. 23. 24. 25.] [13. 23. 33. 34. 35.] [14. 24. 34. 44. 45.] [15. 25. 35. 45. 55.]] R8PP_ZEROS_TEST R8PP_ZEROS sets an R8PP zero matrix. Matrix order N = 5 The R8PP zero matrix: Col: 0 1 2 3 4 Row 0 : 0 0 0 0 0 1 : 0 0 0 0 0 2 : 0 0 0 0 0 3 : 0 0 0 0 0 4 : 0 0 0 0 0 r8pp_test(): Normal end of execution. Fri Aug 12 23:49:06 2022