We present an OpenCL-based Lattice QCD application using a heatbath algorithm for the pure gauge case and Wilson fermions in the twisted mass formulation. The implementation is platform independent and can be used on AMD or NVIDIA GPUs, as well as on classical CPUs. On the AMD Radeon HD 5870 our double precision dslash implementation performs at 60 GFLOPS over a wide range of lattice sizes. The hybrid Monte-Carlo presented reaches a speedup of four over the reference code running on a server CPU.
@article{arxiv.1209.5942,
title = {Lattice QCD based on OpenCL},
author = {Matthias Bach and Volker Lindenstruth and Owe Philipsen and Christopher Pinke},
journal= {arXiv preprint arXiv:1209.5942},
year = {2013}
}