Wed Aug 17 10:49:54 2022 r83s_test(): Python version: 3.6.9 Test r83s(). r83s_cg_test(): r83s_cg() applies CG to an R83S matrix. Number of variables N = 10 Norm of residual ||Ax-b|| = 5.789749205971595e-15 Norm of error ||x1-x2|| = 1.8317102817086012e-15 R83S_DIF2_TEST R83S_DIF2 sets an R83S matrix to the second difference. We check three cases, MN. Second difference in R83S format: Col: 1 2 3 4 5 Row --- 1: 2 -1 2: -1 2 -1 Second difference in R83S format: Col: 1 2 3 4 5 Row --- 1: 2 -1 2: -1 2 -1 3: -1 2 -1 4: -1 2 -1 Second difference in R83S format: Col: 1 2 3 Row --- 1: 2 -1 2: -1 2 -1 3: -1 2 4: -1 r83s_gs_sl_test(): r83s_gs_sl() applies Gauss-Seidel iteration with an R83S matrix to solve a linear system A*x=b. Number of variables N = 10 Norm of residual ||Ax-b|| = 2.5628156264697524e-05 Norm of error ||x1-x2|| = 0.0003003720827043647 R83S_INDICATOR_TEST R83S_INDICATOR sets an R83S indicator matrix. We check three cases, MN. R83S indicator matrix: Col: 1 2 3 4 5 Row --- 1: 2 3 2: 1 2 3 R83S indicator matrix: Col: 1 2 3 4 5 Row --- 1: 2 3 2: 1 2 3 3: 1 2 3 4: 1 2 3 R83S indicator matrix: Col: 1 2 3 Row --- 1: 2 3 2: 1 2 3 3: 1 2 4: 1 R83S_JAC_SL_TEST R83S_JAC_SL applies Jacobi iteration with an R83S matrix to solve a linear system A*x=b. Number of variables N = 10 Norm of residual ||Ax-b|| = 0.010102077770625913 Norm of error ||x1-x2|| = 0.026428262147677625 r83s_mtv_test(): r83s_mtv() computes b=A'*x, where A is an R83S matrix. We check three cases, MN. Product comparison: 0: 1.87234 1.87234 1: 3.23824 3.23824 2: 1.80142 1.80142 3: 0.582736 0.582736 4: 0 0 Product comparison: 0: 0.843781 0.843781 1: 2.4545 2.4545 2: 4.06522 4.06522 3: 5.67594 5.67594 4: 6.54692 6.54692 Product comparison: 0: 1.94535 1.94535 1: 4.25871 4.25871 2: 6.57207 6.57207 r83s_mv_test(): r83s_mv() computes b=A*x, where A is an R83S matrix. We check three cases, MN. Product comparison: 0: 2.02367 2.02367 1: 3.9557 3.9557 2: 5.88773 5.88773 Product comparison: 0: 2.22745 2.22745 1: 3.9762 3.9762 2: 5.72495 5.72495 3: 7.47371 7.47371 4: 5.51173 5.51173 Product comparison: 0: 1.29379 1.29379 1: 2.71268 2.71268 2: 2.83783 2.83783 3: 1.34564 1.34564 4: 0 0 R83S_PRINT_TEST R83S_PRINT prints an R83S matrix. R83S matrix: Col: 1 2 3 4 5 Row --- 1: 2 3 2: 1 2 3 3: 1 2 3 4: 1 2 3 R83S_PRINT_SOME_TEST R83S_PRINT_SOME prints some of an R83S matrix. Rows 2-5, Cols 2-4: Col: 2 3 4 Row --- 2: 2 3 3: 1 2 3 4: 1 2 R83S_RANDOM_TEST R83S_RANDOM randomizes an R83S matrix. We check three cases, MN. Random R83S matrix: Col: 1 2 3 4 5 Row --- 1: 0.511471 0.0412355 2: 0.261536 0.511471 0.0412355 Random R83S matrix: Col: 1 2 3 4 5 Row --- 1: 0.31835 0.516847 2: 0.342483 0.31835 0.516847 3: 0.342483 0.31835 0.516847 4: 0.342483 0.31835 0.516847 Random R83S matrix: Col: 1 2 3 Row --- 1: 0.596542 0.75583 2: 0.552658 0.596542 0.75583 3: 0.552658 0.596542 4: 0.552658 r83s_res_test(): r83s_res() computes b-A*x, where A is an R83S matrix. We check three cases, MN. Residual A*x-b: 0: 0 1: 0 2: 0 Residual A*x-b: 0: 0 1: 0 2: 0 3: 0 4: 0 Residual A*x-b: 0: 0 1: 0 2: 0 3: 0 4: 0 R83S_TO_R8GE_TEST R83S_TO_R8GE converts an R83S matrix to R8GE format. We check three cases, MN. R83s matrix: Col: 1 2 3 4 5 Row --- 1: 0.659474 0.347139 2: 0.916638 0.659474 0.347139 The R8GE matrix: [[0.65947414 0.34713898 0. 0. 0. ] [0.91663776 0.65947414 0.34713898 0. 0. ] [0. 0.91663776 0.65947414 0.34713898 0. ]] R83s matrix: Col: 1 2 3 4 5 Row --- 1: 0.213529 0.393573 2: 0.386897 0.213529 0.393573 3: 0.386897 0.213529 0.393573 4: 0.386897 0.213529 0.393573 The R8GE matrix: [[0.21352936 0.39357291 0. 0. 0. ] [0.38689732 0.21352936 0.39357291 0. 0. ] [0. 0.38689732 0.21352936 0.39357291 0. ] [0. 0. 0.38689732 0.21352936 0.39357291] [0. 0. 0. 0.38689732 0.21352936]] R83s matrix: Col: 1 2 3 Row --- 1: 0.238409 0.52991 2: 0.239696 0.238409 0.52991 3: 0.239696 0.238409 4: 0.239696 The R8GE matrix: [[0.23840871 0.52991008 0. ] [0.2396961 0.23840871 0.52991008] [0. 0.2396961 0.23840871] [0. 0. 0.2396961 ] [0. 0. 0. ]] r83s_zeros_test(): r83s_zeros() zeros an R83S matrix. We check three cases, MN. Zeroed R83S matrix: Col: 1 2 3 4 5 Row --- 1: 0 0 2: 0 0 0 Zeroed R83S matrix: Col: 1 2 3 4 5 Row --- 1: 0 0 2: 0 0 0 3: 0 0 0 4: 0 0 0 Zeroed R83S matrix: Col: 1 2 3 Row --- 1: 0 0 2: 0 0 0 3: 0 0 4: 0 r83s_test(): Normal end of execution. Wed Aug 17 10:49:55 2022