Optimization of the Oktay-Kronfeld Action Conjugate Gradient Inverter
Lattice Gauge Theory Research Center, CTP, and FPRD, Department of Physics and Astronomy, Seoul National University, Seoul, 151-747, South Korea
arXiv:1311.5029 [hep-lat], (20 Nov 2013)
@article{2013arXiv1311.5029J,
author={Jang}, Y.-C. and {Bailey}, J.~A. and {Lee}, W. and {DeTar}, C. and {Oktay}, M.~B. and {Kronfeld}, A.~S.},
title={"{Optimization of the Oktay-Kronfeld Action Conjugate Gradient Inverter}"},
journal={ArXiv e-prints},
archivePrefix={"arXiv"},
eprint={1311.5029},
primaryClass={"hep-lat"},
keywords={High Energy Physics – Lattice},
year={2013},
month={nov},
adsurl={http://adsabs.harvard.edu/abs/2013arXiv1311.5029J},
adsnote={Provided by the SAO/NASA Astrophysics Data System}
}
Improving the Fermilab action to third order in heavy quark effective theory yields the Oktay-Kronfeld action, a promising candidate for precise calculations of the spectra of heavy quark systems and weak matrix elements relevant to searches for new physics. We have optimized the bi-stabilized conjugate gradient inverter in the SciDAC QOPQDP library and are developing a GPU code. The action is rewritten and the needed gauge-link combinations are precalculated. In tests with a MILC coarse lattice, this procedure accelerates the inverter by a factor of four. The remaining floating-point operations are mostly simple matrix multiplications between gauge links and fermion vectors, which we accelerate by more than an order of magnitude by using CUDA. Further gains could be achieved by using QUDA.
November 22, 2013 by hgpu
