Implementation of the Neuberger-Dirac operator on GPUs

INSPIRE: 874632
arXiv: 1010.5636
DOI: 10.22323/1.105.0044

Authors: Walk, Bjoern, Wittig, Hartmut, Dranischnikow, Egor, Schomer, Elmar

Submitted: 2010-10

Subjects:

• hep-lat
• Lattice

Journal reference: PoS LATTICE2010 044 (2010)

Citations: 5

Abstract

Recent developments have shown that a lot can be gained for QCD simulations from GPU hardware. This can be exploited especially in the case of Ginsparg-Wilson fermions when the com putational costs are particularly high. In this work, we use the Neuberger-Dirac operator as our realisation of Ginsparg-Wilson fermions, which greatly facilitate lattice investigations of decays like $K \to \pi\pi$. We report on the ongoing study of our GPU implementation of the Neuberger-Dirac operator including the exact treatment of the low lying eigenmodes of the Wilson-Dirac operator. Our benchmarks show that we achieve speed-up factors of around 23 and 16 in single and double precision, respectively.