Victor N. Datsyuk High Performance Computers Department Rostov State University Rostov-on-Don, 28-54-66 email: root@rsusu1.rnd.runnet.ru Internet: http://rsusu1.rnd.runnet.ru Aztec example driver This program create an MSR matrix corresponding to a 7pt discrete approximation to the 3D Poisson operator on an n x n x n square and solve it various methods. n must be <= 50. input number grid point on each direction n = ********************** Method solve is =cg ********************** Dimension matrices = 27000 NPROCS = 32 iter: 190 residual = 8.417737e-13 Solution time: 5.327474 (sec.) total iterations: 190 scale_flops: 0.000000e+00 iter_flops: 1.333556e+08 precond_flops: 0.000000e+00 total_flops: 1.333556e+08 Solver MFlop rate: 25.031681 MFlops/sec. Solver processor MFlop rate: 0.782240 MFlops/sec. i= 0 x(i)= 1.00000000 i= 26999 x(i)= 27000.00000000 ********************** Method solve is =gmres ********************** Dimension matrices = 27000 NPROCS = 32 iter: 553 residual = 9.542542e-13 Solution time: 26.872881 (sec.) total iterations: 553 scale_flops: 0.000000e+00 iter_flops: 1.611694e+09 precond_flops: 0.000000e+00 total_flops: 1.611694e+09 Solver MFlop rate: 59.974725 MFlops/sec. Solver processor MFlop rate: 1.874210 MFlops/sec. i= 0 x(i)= 0.99999997 i= 26999 x(i)= 26999.99999998 ********************** Method solve is =cgs ********************** Dimension matrices = 27000 NPROCS = 32 iter: 119 residual = 7.539097e-13 Solution time: 6.317644 (sec.) total iterations: 119 scale_flops: 0.000000e+00 iter_flops: 1.447851e+08 precond_flops: 0.000000e+00 total_flops: 1.447851e+08 Solver MFlop rate: 22.917582 MFlops/sec. Solver processor MFlop rate: 0.716174 MFlops/sec. i= 0 x(i)= 1.00000000 i= 26999 x(i)= 27000.00000000 ********************** Method solve is =tfqmr ********************** Dimension matrices = 27000 NPROCS = 32 iter: 119 residual = 2.636488e-13 Solution time: 6.803927 (sec.) total iterations: 119 scale_flops: 0.000000e+00 iter_flops: 1.709193e+08 precond_flops: 0.000000e+00 total_flops: 1.709193e+08 Solver MFlop rate: 25.120677 MFlops/sec. Solver processor MFlop rate: 0.785021 MFlops/sec. i= 0 x(i)= 1.00000000 i= 26999 x(i)= 27000.00000000 ********************** Method solve is =bicgstab ********************** Dimension matrices = 27000 NPROCS = 32 iter: 126 residual = 8.453235e-13 Solution time: 6.886447 (sec.) total iterations: 126 scale_flops: 0.000000e+00 iter_flops: 1.566413e+08 precond_flops: 0.000000e+00 total_flops: 1.566413e+08 Solver MFlop rate: 22.746314 MFlops/sec. Solver processor MFlop rate: 0.710822 MFlops/sec. i= 0 x(i)= 0.99999999 i= 26999 x(i)= 26999.99999999