22-Aug-2022 19:16:22 r85_test(): MATLAB version Test r85(). R85_DIF2_TEST R85_DIF2 sets up a R85 second difference matrix. Matrix order N = 5 The R85 second difference matrix: Col: 1 2 3 4 5 Row --- 1 2 -1 0 2 -1 2 -1 0 3 0 -1 2 -1 0 4 0 -1 2 -1 5 0 -1 2 R85_INDICATOR_TEST R85_INDICATOR sets up a R85 indicator matrix. Matrix order N = 10 The R85 indicator matrix: Col: 1 2 3 4 5 Row --- 1 101 102 103 2 201 202 203 204 3 301 302 303 304 305 4 402 403 404 405 5 503 504 505 6 604 605 7 705 Col: 6 7 8 9 10 Row --- 4 406 5 506 507 6 606 607 608 7 706 707 708 709 8 806 807 808 809 810 9 907 908 909 910 10 1008 1009 1010 R85_MTV_TEST R85_MTV computes b=A'*x, where A is an R85 matrix. Matrix order N = 5 The R85 matrix: Col: 1 2 3 4 5 Row --- 1 11 12 13 2 21 22 23 24 3 31 32 33 34 35 4 42 43 44 45 5 53 54 55 The vector x: 1: 1 2: 2 3: 3 4: 4 5: 5 The product b=A'*x: 1: 1369 2: 1799 3: 2073 4: 1763 5: 1309 R85_MV_TEST R85_MV computes b=A*x, where A is an R85 matrix. Matrix order N = 5 The R85 matrix: Col: 1 2 3 4 5 Row --- 1 11 12 13 2 21 22 23 24 3 31 32 33 34 35 4 42 43 44 45 5 53 54 55 The vector x: 1: 1 2: 2 3: 3 4: 4 5: 5 The product b=A*x: 1: 146 2: 320 3: 595 4: 596 5: 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: 1 2 3 4 5 Row --- 1 0.350752 0.043829 0.218418 2 0.351629 0.0945448 0.633966 0.956318 3 0.188955 0.822887 0.0136169 0.0617272 0.829509 4 0.761492 0.267132 0.859097 0.449539 5 0.396988 0.692066 0.840847 6 0.185314 0.561662 7 0.574366 Col: 6 7 8 9 10 Row --- 4 0.561695 5 0.401306 0.415307 6 0.123104 0.754673 0.0661187 7 0.861216 0.00751236 0.797287 0.257578 8 0.367027 0.453794 0.260303 0.00183837 0.109957 9 0.617205 0.911977 0.912484 0.897504 10 0.361529 0.597917 0.113664 The right hand side: 1: 1.62087 2: 5.74755 3: 4.58067 4: 10.1064 5: 15.8815 6: 13.9567 7: 15.8424 8: 19.8832 9: 16.0093 10: 10.0938 Solution: 1: 1 2: 2 3: 3 4: 4 5: 5 6: 6 7: 7 8: 8 9: 9 10: 10 R85_PRINT_TEST R85_PRINT prints an R85 matrix. Matrix order N = 5 The R85 matrix: Col: 1 2 3 4 5 Row --- 1 11 12 13 2 21 22 23 24 3 31 32 33 34 35 4 42 43 44 45 5 53 54 55 R85_PRINT_SOME_TEST R85_PRINT_SOME prints some of an R85 matrix. Matrix order N = 8 Rows 2-6, Cols 3-5: Col: 3 4 5 Row --- 2 23 24 3 33 34 35 4 43 44 45 5 53 54 55 6 64 65 R85_RANDOM_TEST R85_RANDOM sets up a random R85 matrix. Matrix order N = 5 The R85 random matrix: Col: 1 2 3 4 5 Row --- 1 0.109957 0.561695 0.218418 2 0.401306 0.043829 0.415307 0.956318 3 0.897504 0.754673 0.633966 0.0661187 0.829509 4 0.350752 0.797287 0.0617272 0.257578 5 0.0945448 0.00183837 0.449539 R85_TO_R8GE_TEST R85_TO_R8GE converts an R85 matrix to R8GE form. Matrix order N = 5 The R85 matrix: Col: 1 2 3 4 5 Row --- 1 11 12 13 2 21 22 23 24 3 31 32 33 34 35 4 42 43 44 45 5 53 54 55 The R8GE matrix: Col: 1 2 3 4 5 Row --- 1 11 21 31 0 0 2 12 22 32 42 0 3 13 23 33 43 53 4 0 24 34 44 54 5 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: 1 2 3 4 5 Row --- 1 0 0 0 2 0 0 0 0 3 0 0 0 0 0 4 0 0 0 0 5 0 0 0 r85_test(): Normal end of execution. 22-Aug-2022 19:16:22