Block-Iterative Methods for 3D Constant-Coefficient Stencils on GPUs and Multicore CPUs

Technical Report · Sun Jun 01 00:00:00 EDT 2014

DOI:https://doi.org/10.2172/1134156· OSTI ID:1134156

Philip, Bobby ^[1]; Wang, Zhen ^[1]; Berrill, Mark A ^[1]

---

1. ORNL

Block iterative methods are extremely important as smoothers for multigrid methods, as preconditioners for Krylov methods, and as solvers for diagonally dominant linear systems. Developing robust and efficient smoother algorithms suitable for current and evolving GPU and multicore CPU systems is a significant challenge. We address this issue in the case of constant-coefficient stencils arising in the solution of elliptic partial differential equations on structured 3D uniform and adaptively refined block structured grids. Robust, highly parallel implementations of block Jacobi and chaotic block Gauss-Seidel algorithms with exact inversion of the blocks are developed using different parallelization techniques. Experimental results for NVIDIA Fermi/Kepler GPUs and AMD multicore systems are presented.

• Research Organization: Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States). Oak Ridge Leadership Computing Facility (OLCF)
• Sponsoring Organization: USDOE Office of Science (SC)
• DOE Contract Number: DE-AC05-00OR22725
• OSTI ID: 1134156
• Report Number(s): ORNL/TM-2014/225; KJ0403000; ERKJR04
• Country of Publication: United States
• Language: English