Tue May 20 22:27:32 2025 r83s_test(): python version: 3.10.12 numpy version: 1.26.4 Test r83s(). r83s_cg_test(): r83s_cg() applies CG to an R83S matrix. Number of variables N = 10 Norm of residual ||Ax-b|| = 1.3002420749351734e-15 Norm of error ||x1-x2|| = 7.133574416487974e-16 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.660936391920518e-05 Norm of error ||x1-x2|| = 0.0003118722227750594 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.007964820194528495 Norm of error ||x1-x2|| = 0.02415407546601222 r83s_mtv_test(): r83s_mtv() computes b=A'*x, where A is an R83S matrix. We check three cases, MN. Product comparison: 0: 1.1509 1.1509 1: 2.34085 2.34085 2: 1.51644 1.51644 3: 1.6279 1.6279 4: 0 0 Product comparison: 0: 1.3693 1.3693 1: 2.78302 2.78302 2: 4.19674 4.19674 3: 5.61047 5.61047 4: 3.80476 3.80476 Product comparison: 0: 2.06002 2.06002 1: 3.60086 3.60086 2: 5.14171 5.14171 r83s_mv_test(): r83s_mv() computes b=A*x, where A is an R83S matrix. We check three cases, MN. Product comparison: 0: 2.32734 2.32734 1: 4.35145 4.35145 2: 6.37557 6.37557 Product comparison: 0: 1.21411 1.21411 1: 1.94467 1.94467 2: 2.67523 2.67523 3: 3.40579 3.40579 4: 1.21357 1.21357 Product comparison: 0: 0.46486 0.46486 1: 0.816905 0.816905 2: 0.674183 0.674183 3: 0.0326305 0.0326305 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.439827 0.978762 2: 0.439929 0.439827 0.978762 Random R83S matrix: Col: 1 2 3 4 5 Row --- 1: 0.922904 0.899384 2: 0.021126 0.922904 0.899384 3: 0.021126 0.922904 0.899384 4: 0.021126 0.922904 0.899384 Random R83S matrix: Col: 1 2 3 Row --- 1: 0.506757 0.107631 2: 0.748784 0.506757 0.107631 3: 0.748784 0.506757 4: 0.748784 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.307038 0.968286 2: 0.644583 0.307038 0.968286 The R8GE matrix: [[0.30703761 0.96828603 0. 0. 0. ] [0.64458317 0.30703761 0.96828603 0. 0. ] [0. 0.64458317 0.30703761 0.96828603 0. ]] R83s matrix: Col: 1 2 3 4 5 Row --- 1: 0.538445 0.0181025 2: 0.39307 0.538445 0.0181025 3: 0.39307 0.538445 0.0181025 4: 0.39307 0.538445 0.0181025 The R8GE matrix: [[0.53844535 0.01810246 0. 0. 0. ] [0.39306961 0.53844535 0.01810246 0. 0. ] [0. 0.39306961 0.53844535 0.01810246 0. ] [0. 0. 0.39306961 0.53844535 0.01810246] [0. 0. 0. 0.39306961 0.53844535]] R83s matrix: Col: 1 2 3 Row --- 1: 0.133695 0.868833 2: 0.543531 0.133695 0.868833 3: 0.543531 0.133695 4: 0.543531 The R8GE matrix: [[0.13369494 0.86883348 0. ] [0.54353076 0.13369494 0.86883348] [0. 0.54353076 0.13369494] [0. 0. 0.54353076] [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. Tue May 20 22:27:32 2025