GPU-Based Conjugate Gradient Solver for Lattice QCD with Domain-Wall Fermions
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)
@article{2011arXiv1101.0423C,
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},
archivePrefix={“arXiv”},
eprint={1101.0423},
primaryClass={“hep-lat”},
keywords={High Energy Physics – Lattice},
year={2011},
month={jan},
adsurl={http://adsabs.harvard.edu/abs/2011arXiv1101.0423C},
adsnote={Provided by the SAO/NASA Astrophysics Data System}
}
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).
January 4, 2011 by hgpu
