17 August 2022 07:26:01 AM r83s_test(): C version Test r83s(). R83S_CG_TEST R83S_CG applies CG to an R83S matrix. Number of variables N = 10 Norm of residual ||Ax-b|| = 8.47455e-16 Norm of error ||x1-x2|| = 4.1616e-16 R83S_DIF2_TEST R83S_DIF2 sets up the second difference matrix in R83S format. We check three cases, MN. The R83S matrix: Col: 1 2 3 4 5 Row --- 1 2.000000 -1.000000 2 -1.000000 2.000000 -1.000000 3 -1.000000 2.000000 -1.000000 The R83S matrix: Col: 1 2 3 4 5 Row --- 1 2.000000 -1.000000 2 -1.000000 2.000000 -1.000000 3 -1.000000 2.000000 -1.000000 4 -1.000000 2.000000 -1.000000 5 -1.000000 2.000000 The R83S matrix: Col: 1 2 3 Row --- 1 2.000000 -1.000000 2 -1.000000 2.000000 -1.000000 3 -1.000000 2.000000 4 -1.000000 R83S_GS_SL_TEST R83S_GS_SL uses Gauss-Seidel iteration on an R83S matrix. Matrix order N = 10 Iterations per call = 25 Current estimated solution: 0 0.634606 1 1.326504 2 2.095665 3 2.954509 4 3.907410 5 4.950882 6 6.074424 7 7.261901 8 8.493295 9 9.746648 Current estimated solution: 0 0.953567 1 1.914504 2 2.885328 3 3.867570 4 4.861732 5 5.867333 6 6.883019 7 7.906746 8 8.935991 9 9.967996 Current estimated solution: 0 0.994126 1 1.989185 2 2.985494 3 3.983248 4 4.982509 5 5.983218 6 6.985202 7 7.988203 8 8.991903 9 9.995951 R83S_INDICATOR_TEST R83S_INDICATOR sets up an R83S indicator matrix. The R83S indicator matrix: Col: 1 2 3 4 Row --- 1 2.000000 3.000000 2 1.000000 2.000000 3.000000 3 1.000000 2.000000 3.000000 4 1.000000 2.000000 5 1.000000 R83S_JAC_SL_TEST R83S_JAC_SL solves a linear system using Jacobi iteration, for an R83S matrix. Matrix order N = 10 Iterations per call = 25 Current estimated solution: 0 0.315171 1 0.727797 2 1.140423 3 1.827580 4 2.514737 5 3.590467 6 4.666197 7 6.128201 8 7.590205 9 9.295102 Current estimated solution: 0 0.757545 1 1.515091 2 2.349361 3 3.183631 4 4.147287 5 5.110943 6 6.215812 7 7.320681 8 8.533664 9 9.746648 Current estimated solution: 0 0.910021 1 1.834323 2 2.758626 3 3.721242 4 4.683859 5 5.696660 6 6.709462 7 7.768388 8 8.827314 9 9.913657 R83S_MTV_TEST R83S_MTV computes b=A'*x, where A is an R83S matrix. We check three cases, MN. Product comparison: 0: 2.61534 2.61534 1: 4.61958 4.61958 2: 3.30579 3.30579 3: 0.655255 0.655255 4: 0 0 Product comparison: 0: 2.61534 2.61534 1: 4.61958 4.61958 2: 6.62383 6.62383 3: 8.62807 8.62807 4: 5.65526 5.65526 Product comparison: 0: 2.61534 2.61534 1: 4.61958 4.61958 2: 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: 0: 1.39315 1.39315 1: 3.3974 3.3974 2: 5.40164 5.40164 Product comparison: 0: 1.39315 1.39315 1: 3.3974 3.3974 2: 5.40164 5.40164 3: 7.40589 7.40589 4: 8.09962 8.09962 Product comparison: 0: 1.39315 1.39315 1: 3.3974 3.3974 2: 4.52797 4.52797 3: 2.48853 2.48853 4: 0 0 R83S_PRINT_TEST R83S_PRINT prints an R83S matrix. The R83S matrix: Col: 1 2 3 4 Row --- 1 2.000000 3.000000 2 1.000000 2.000000 3.000000 3 1.000000 2.000000 3.000000 4 1.000000 2.000000 5 1.000000 R83S_PRINT_SOME_TEST R83S_PRINT_SOME prints some of an R83S matrix. Rows 1-4, Cols 1-3: Col: 1 2 3 Row --- 1 2.000000 3.000000 2 1.000000 2.000000 3.000000 3 1.000000 2.000000 4 1.000000 R83S_RANDOM_TEST R83S_RANDOM randomizes an R83S matrix. The R83 matrix: Col: 1 2 3 4 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 5 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: 0 0.000000 1 0.000000 2 0.000000 Residual A*x-b: 0 0.000000 1 0.000000 2 0.000000 3 0.000000 4 0.000000 Residual A*x-b: 0 0.000000 1 0.000000 2 0.000000 3 0.000000 4 0.000000 R83S_TO_R8GE_TEST R83S_TO_R8GE converts an R83S matrix to R8GE format. We check three cases, MN. The 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 The 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 The 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 The 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 The 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 The 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. The R83S matrix: Col: 1 2 3 4 Row --- 1 0.000000 0.000000 2 0.000000 0.000000 0.000000 3 0.000000 0.000000 0.000000 4 0.000000 0.000000 5 0.000000 r83s_test(): Normal end of execution. 17 August 2022 07:26:01 AM