GPU-Based Conjugate Gradient Solver for Lattice QCD with Domain-Wall Fermions

Ting-Wai Chiu, Tung-Han Hsieh, Yao-Yuan Mao, Kenji Ogawa
TWQCD Collaboration, Department of Physics, and Center for Theoretical Sciences, National Taiwan University, Taipei 10617, Taiwan
The XXVIII International Symposium on Lattice Field Theory, Lattice 2010, June 14-19, 2010, arXiv:1101.0423 [hep-lat] (2 Jan 2011)


   author={Chiu}, {T.-W.} and {Hsieh}, {T.-H.} and {Mao}, {Y.-Y.} and {Ogawa}, K.},

   title={“{GPU-Based Conjugate Gradient Solver for Lattice QCD with Domain-Wall Fermions}”},

   journal={ArXiv e-prints},




   keywords={High Energy Physics – Lattice},




   adsnote={Provided by the SAO/NASA Astrophysics Data System}


Download Download (PDF)   View View   Source Source   



We present the first GPU-based conjugate gradient (CG) solver for lattice QCD with domain-wall fermions (DWF). It is well-known that CG is the most time-consuming part in the Hybrid Monte Carlo simulation of unquenched lattice QCD, which becomes even more computational demanding for lattice QCD with exact chiral symmetry. We have designed a CG solver for the general 5-dimensional DWF operator on NVIDIA CUDA architecture with mixed-precision, using the defect correction as well as the reliable updates algorithms. We optimize our computation by even-odd preconditioning in the 4D space-time lattice, plus several innovative techniques for CUDA kernels. For NVIDIA GeForce GTX 285/480, our CG solver attains 180/233 Gflops (sustained).
No votes yet.
Please wait...

* * *

* * *

HGPU group © 2010-2021 hgpu.org

All rights belong to the respective authors

Contact us: