17-Aug-2022 15:46:01 R83T_TEST MATLAB version Test the R83T library. R83T_CG_TEST R83T_CG applies CG to an R83T matrix. x1 = 0.2184 0.9563 0.8295 0.5617 0.4153 0.0661 0.2576 0.1100 0.0438 0.6340 Number of variables N = 10 Norm of residual ||Ax-b|| = 1.60815e-15 Norm of error ||x1-x2|| = 5.41056e-16 R83T_DIF2_TEST R83T_DIF2 sets an R83T matrix to the second difference. R83T matrix: Col: 1 2 3 4 5 Row --- 1: 2 -1 2: -1 2 -1 3: -1 2 -1 4: -1 2 -1 5: -1 2 R83T_GS_SL_TEST R83T_GS_SL solves a linear system using Gauss-Seidel iteration, with R83T matrix storage. Matrix order N = 10 Iterations per call = 25 Current solution estimate: 1: 0.634606 2: 1.3265 3: 2.09566 4: 2.95451 5: 3.90741 6: 4.95088 7: 6.07442 8: 7.2619 9: 8.4933 10: 9.74665 Current solution estimate: 1: 0.953567 2: 1.9145 3: 2.88533 4: 3.86757 5: 4.86173 6: 5.86733 7: 6.88302 8: 7.90675 9: 8.93599 10: 9.968 Current solution estimate: 1: 0.994126 2: 1.98918 3: 2.98549 4: 3.98325 5: 4.98251 6: 5.98322 7: 6.9852 8: 7.9882 9: 8.9919 10: 9.99595 R83T_INDICATOR_TEST R83T_INDICATOR sets an R83T indicator matrix. We check three cases, MN. R83T indicator matrix: Col: 1 2 3 4 5 Row --- 1: 11 12 2: 21 22 23 3: 32 33 34 R83T indicator matrix: Col: 1 2 3 4 5 Row --- 1: 11 12 2: 21 22 23 3: 32 33 34 4: 43 44 45 5: 54 55 R83T indicator matrix: Col: 1 2 3 Row --- 1: 11 12 2: 21 22 23 3: 32 33 4: 43 R83T_JAC_SL_TEST R83T_JAC_SL solves a linear system using Jacobi iteration, with R83T matrix storage. Matrix order N = 10 Iterations per call = 25 Current solution estimate: 1: 0.315171 2: 0.727797 3: 1.14042 4: 1.82758 5: 2.51474 6: 3.59047 7: 4.6662 8: 6.1282 9: 7.5902 10: 9.2951 Current solution estimate: 1: 0.757545 2: 1.51509 3: 2.34936 4: 3.18363 5: 4.14729 6: 5.11094 7: 6.21581 8: 7.32068 9: 8.53366 10: 9.74665 Current solution estimate: 1: 0.910021 2: 1.83432 3: 2.75863 4: 3.72124 5: 4.68386 6: 5.69666 7: 6.70946 8: 7.76839 9: 8.82731 10: 9.91366 R83T_MTV_TEST R83T_MTV multiplies an R83T matrix transposed times a vector. M = 5 N = 6 The R83T matrix A: Col: 1 2 3 4 5 Row --- 1: 11 12 2: 21 22 23 3: 32 33 34 4: 43 44 45 5: 54 55 Col: 6 Row --- 5: 56 The vector x: 1: 1 2: 2 3: 3 4: 4 5: 5 The product b = A'*x: 1: 53 2: 152 3: 317 4: 548 5: 455 6: 280 R83T_MV_TEST R83T_MV multiplies an R83T matrix times a vector. M = 5 N = 6 The R83T matrix A: Col: 1 2 3 4 5 Row --- 1: 11 12 2: 21 22 23 3: 32 33 34 4: 43 44 45 5: 54 55 Col: 6 Row --- 5: 56 The vector x: 1: 1 2: 2 3: 3 4: 4 5: 5 6: 6 The product b = A*x: 1: 35 2: 134 3: 299 4: 530 5: 827 R83T_PRINT_TEST R83T_PRINT prints an R83T matrix. We check three cases, MN. R83T matrix: Col: 1 2 3 4 5 Row --- 1: 11 12 2: 21 22 23 3: 32 33 34 R83T matrix: Col: 1 2 3 4 5 Row --- 1: 11 12 2: 21 22 23 3: 32 33 34 4: 43 44 45 5: 54 55 R83T matrix: Col: 1 2 3 Row --- 1: 11 12 2: 21 22 23 3: 32 33 4: 43 R83T_PRINT_SOME_TEST R83T_PRINT_SOME prints some of an R83T matrix. M = 9 N = 9 Rows 3:6, Cols 5:8: Col: 5 6 7 8 Row --- 4: 45 5: 55 56 6: 65 66 67 R83T_RANDOM_TEST R83T_RANDOM returns a random R83T matrix. R83T matrix: Col: 1 2 3 4 5 Row --- 1: 0.218418 0.956318 2: 0.829509 0.561695 0.415307 3: 0.0661187 0.257578 0.109957 4: 0.043829 0.633966 0.0617272 5: 0.449539 0.401306 R83T_RES_TEST R83T_RES evaluates the residual given an approximate solution of a linear system A*x=b. M = 5 N = 5 The R83T matrix A: Col: 1 2 3 4 5 Row --- 1: 2 -1 2: -1 2 -1 3: -1 2 -1 4: -1 2 -1 5: -1 2 The right hand side B: 1: 0 2: 0 3: 0 4: 0 5: 6 The solution X: 1: 1 2: 2 3: 3 4: 4 5: 5 The residual b-A*x: 1: 2.22045e-16 2: -4.44089e-16 3: 8.88178e-16 4: 0 5: 0 R83T_TO_R8GE_TEST R83T_TO_R8GE converts an R83T matrix to R8GE format. M = 5 N = 5 The R83T indicator matrix: Col: 1 2 3 4 5 Row --- 1: 11 12 2: 21 22 23 3: 32 33 34 4: 43 44 45 5: 54 55 The R8GE format matrix: Col: 1 2 3 4 5 Row --- 1 11 12 0 0 0 2 21 22 23 0 0 3 0 32 33 34 0 4 0 0 43 44 45 5 0 0 0 54 55 R83T_ZEROS_TEST R83T_ZEROS sets an R83T matrix to zero. R83T matrix: Col: 1 2 3 4 5 Row --- 1: 0 0 2: 0 0 0 3: 0 0 0 4: 0 0 0 5: 0 0 R83T_TEST Normal end of execution. 17-Aug-2022 15:46:02