17-Aug-2022 08:54:55 r83s_test(): MATLAB version Test r83s(). R83S_CG_TEST R83S_CG applies CG to an R83S matrix. Number of variables N = 10 Norm of residual ||Ax-b|| = 1.60815e-15 Norm of error ||x1-x2|| = 5.41056e-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 3: -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 5: -1 2 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|| = 4.51901e-05 Norm of error ||x1-x2|| = 0.000529646 R83S_INDICATOR_TEST R83S_INDICATOR sets an R83S indicator matrix. We check three cases, MN. a = 3 2 1 R83S indicator matrix: Col: 1 2 3 4 5 Row --- 1: 2 3 2: 1 2 3 3: 1 2 3 a = 3 2 1 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 5: 1 2 a = 3 2 1 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.0023461 Norm of error ||x1-x2|| = 0.0289483 R83S_MTV_TEST R83S_MV computes b=A'*x, where A is an R83S matrix. We check three cases, MN. Product comparison: 1: 2.61534 2.61534 2: 4.61958 4.61958 3: 3.30579 3.30579 4: 0.655255 0.655255 5: 0 0 Product comparison: 1: 2.61534 2.61534 2: 4.61958 4.61958 3: 6.62383 6.62383 4: 8.62807 8.62807 5: 5.65526 5.65526 Product comparison: 1: 2.61534 2.61534 2: 4.61958 4.61958 3: 6.62383 6.62383 R83S_MV_TEST R83S_MV computes b=A*x, where A is an R83S matrix. We check three cases, MN. Product comparison: 1: 1.39315 1.39315 2: 3.3974 3.3974 3: 5.40164 5.40164 Product comparison: 1: 1.39315 1.39315 2: 3.3974 3.3974 3: 5.40164 5.40164 4: 7.40589 7.40589 5: 8.09962 8.09962 Product comparison: 1: 1.39315 1.39315 2: 3.3974 3.3974 3: 4.52797 4.52797 4: 2.48853 2.48853 5: 0 0 R83S_PRINT_TEST R83S_PRINT prints an R83S matrix. a = 3 2 1 R83S matrix: Col: 1 2 3 4 5 Row --- 1: 2 3 2: 1 2 3 3: 1 2 3 4: 1 2 3 5: 1 2 R83S_PRINT_SOME_TEST R83S_PRINT_SOME prints some of an R83S matrix. a = 3 2 1 Rows 2-5, Cols 2-4: Col: 2 3 4 Row --- 2: 2 3 3: 1 2 3 4: 1 2 5: 1 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.956318 0.218418 2: 0.829509 0.956318 0.218418 3: 0.829509 0.956318 0.218418 Random R83S matrix: Col: 1 2 3 4 5 Row --- 1: 0.956318 0.218418 2: 0.829509 0.956318 0.218418 3: 0.829509 0.956318 0.218418 4: 0.829509 0.956318 0.218418 5: 0.829509 0.956318 Random R83S matrix: Col: 1 2 3 Row --- 1: 0.956318 0.218418 2: 0.829509 0.956318 0.218418 3: 0.829509 0.956318 4: 0.829509 R83S_RES_TEST R83S_RES computes b-A*x, where A is an R83S matrix. We check three cases, MN. Residual A*x-b: 1: 0 2: 0 3: 0 Residual A*x-b: 1: 0 2: 0 3: 0 4: 0 5: 0 Residual A*x-b: 1: 0 2: 0 3: 0 4: 0 5: 0 R83S_TO_R8GE_TEST R83S_TO_R8GE converse an R83S matrix to R8GE format. We check three cases, MN. R83s matrix: Col: 1 2 3 4 5 Row --- 1: 0.956318 0.218418 2: 0.829509 0.956318 0.218418 3: 0.829509 0.956318 0.218418 R8GE matrix: Col: 1 2 3 4 5 Row --- 1 0.956318 0.218418 0 0 0 2 0.829509 0.956318 0.218418 0 0 3 0 0.829509 0.956318 0.218418 0 R83s matrix: Col: 1 2 3 4 5 Row --- 1: 0.956318 0.218418 2: 0.829509 0.956318 0.218418 3: 0.829509 0.956318 0.218418 4: 0.829509 0.956318 0.218418 5: 0.829509 0.956318 R8GE matrix: Col: 1 2 3 4 5 Row --- 1 0.956318 0.218418 0 0 0 2 0.829509 0.956318 0.218418 0 0 3 0 0.829509 0.956318 0.218418 0 4 0 0 0.829509 0.956318 0.218418 5 0 0 0 0.829509 0.956318 R83s matrix: Col: 1 2 3 Row --- 1: 0.956318 0.218418 2: 0.829509 0.956318 0.218418 3: 0.829509 0.956318 4: 0.829509 R8GE matrix: Col: 1 2 3 Row --- 1 0.956318 0.218418 0 2 0.829509 0.956318 0.218418 3 0 0.829509 0.956318 4 0 0 0.829509 5 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 3: 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 5: 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. 17-Aug-2022 08:54:56