Tue May 20 22:27:32 2025 r85_test(): python version: 3.10.12 numpy version: 1.26.4 Test r85(). R85_DIF2_TEST R85_DIF2 sets up a R85 second difference matrix. Matrix order N = 5 The R85 second difference matrix: Col: 0 1 2 3 4 Row 0 : 2 -1 0 1 : -1 2 -1 0 2 : 0 -1 2 -1 0 3 : 0 -1 2 -1 4 : 0 -1 2 R85_INDICATOR_TEST R85_INDICATOR sets up a R85 indicator matrix. Matrix order N = 5 The R85 indicator matrix: Col: 0 1 2 3 4 Row 0 : 11 12 13 1 : 21 22 23 24 2 : 31 32 33 34 35 3 : 42 43 44 45 4 : 53 54 55 R85_MTV_TEST R85_MTV computes b=A'*x, where A is an R85 matrix. Matrix order N = 5 The R85 matrix: Col: 0 1 2 3 4 Row 0 : 11 12 13 1 : 21 22 23 24 2 : 31 32 33 34 35 3 : 42 43 44 45 4 : 53 54 55 The vector x: 0: 1 1: 2 2: 3 3: 4 4: 5 The product b=A'*x: 0: 74 1: 230 2: 505 3: 614 4: 650 R85_MV_TEST R85_MV computes b=A*x, where A is an R85 matrix. Matrix order N = 5 The R85 matrix: Col: 0 1 2 3 4 Row 0 : 11 12 13 1 : 21 22 23 24 2 : 31 32 33 34 35 3 : 42 43 44 45 4 : 53 54 55 The vector x: 0: 1 1: 2 2: 3 3: 4 4: 5 The product b=A*x: 0: 146 1: 320 2: 595 3: 596 4: 560 R85_NP_FS_TEST R85_NP_FS factors and solves a pentadiagonal linear system, with no pivoting. Matrix order N = 10 The pentadiagonal matrix: Col: 0 1 2 3 4 Row 0 : 0.911739 0.0212587 0.793181 1 : 0.575954 0.395773 0.253246 0.141324 2 : 0.646459 0.220699 0.874481 0.689574 0.943795 3 : 0.828802 0.472467 0.872075 0.303233 4 : 0.934327 0.814036 0.309674 5 : 0.894886 0.882205 6 : 0.195615 Col: 5 6 7 8 9 Row 3 : 0.673755 4 : 0.651791 0.892228 5 : 0.460902 0.0522266 0.261903 6 : 0.688416 0.919644 0.490387 0.0634863 7 : 0.705876 0.760082 0.633002 0.28624 0.175203 8 : 0.362157 0.700525 0.0624774 0.213134 9 : 0.521918 0.988201 0.770979 The right hand side: 0: 4.00302 1: 4.79011 2: 10.4846 3: 15.2792 4: 12.2552 5: 19.1853 6: 20.5521 7: 21.592 8: 13.1786 9: 11.0296 Solution: 0: 1 1: 2 2: 3 3: 4 4: 5 5: 6 6: 7 7: 8 8: 9 9: 10 R85_PRINT_TEST R85_PRINT prints an R85 matrix. Matrix order N = 5 The R85 matrix: Col: 0 1 2 3 4 Row 0 : 11 12 13 1 : 21 22 23 24 2 : 31 32 33 34 35 3 : 42 43 44 45 4 : 53 54 55 R85_PRINT_SOME_TEST R85_PRINT_SOME prints some of an R85 matrix. Matrix order N = 8 Rows 1-5, Cols 2-4: Col: 2 3 4 Row 1 : 23 24 2 : 33 34 35 3 : 43 44 45 4 : 53 54 55 5 : 64 65 R85_RANDOM_TEST R85_RANDOM sets up a random R85 matrix. Matrix order N = 5 The R85 random matrix: Col: 0 1 2 3 4 Row 0 : 0.729483 0.82341 0.965857 1 : 0.209499 0.284322 0.100176 0.0714103 2 : 0.271463 0.35086 0.850101 0.757701 0.440329 3 : 0.466169 0.545937 0.86749 0.0857228 4 : 0.506214 0.531601 0.395379 R85_TO_R8GE_TEST R85_TO_R8GE converts an R85 matrix to R8GE form. Matrix order N = 5 The R85 matrix: Col: 0 1 2 3 4 Row 0 : 11 12 13 1 : 21 22 23 24 2 : 31 32 33 34 35 3 : 42 43 44 45 4 : 53 54 55 The R8GE matrix: [[11. 21. 31. 0. 0.] [12. 22. 32. 42. 0.] [13. 23. 33. 43. 53.] [ 0. 24. 34. 44. 54.] [ 0. 0. 35. 45. 55.]] R85_ZEROS_TEST R85_ZEROS sets up a R85 zero matrix. Matrix order N = 5 The R85 zero matrix: Col: 0 1 2 3 4 Row 0 : 0 0 0 1 : 0 0 0 0 2 : 0 0 0 0 0 3 : 0 0 0 0 4 : 0 0 0 r85_test(): Normal end of execution. Tue May 20 22:27:32 2025