Hybrid Monte Carlo with Wilson Dirac operator on the Fermi GPU
Electra Design Automation, SDF Building, SaltLake Sec-V, Kolkata – 700091
arXiv:1207.2223v1 [hep-lat] (10 Jul 2012)
@article{2012arXiv1207.2223C,
author={Chakrabarty}, A. and {Majumdar}, P.},
title={"{Hybrid Monte Carlo with Wilson Dirac operator on the Fermi GPU}"},
journal={ArXiv e-prints},
archivePrefix={"arXiv"},
eprint={1207.2223},
primaryClass={"hep-lat"},
keywords={High Energy Physics – Lattice},
year={2012},
month={jul},
adsurl={http://adsabs.harvard.edu/abs/2012arXiv1207.2223C},
adsnote={Provided by the SAO/NASA Astrophysics Data System}
}
In this article we present our implementation of a Hybrid Monte Carlo algorithm for Lattice Gauge Theory using two degenerate flavours of Wilson-Dirac fermions on a Fermi GPU. We find that using registers instead of global memory speeds up the code by almost an order of magnitude. To map the array variables to scalars, so that the compiler puts them in the registers, we use code generators. Our final program is more than 10 times faster than a generic single CPU.
July 11, 2012 by hgpu
