Implementation of the conjugate gradient algorithm for heterogeneous systems

INSPIRE: 1980671
arXiv: 2111.14958
DOI: 10.22323/1.396.0507

Authors: Cali, Salvatore, Detmold, William, Korcyl, Grzegorz, Korcyl, Piotr, Shanahan, Phiala

Submitted: 29 November 2021

Subjects:

• hep-lat
• Lattice

Journal reference: PoS LATTICE2021 507 (2022)

Citations: 1

Abstract

Lattice QCD calculations require significant computational effort, with the dominant fraction of resources typically spent in the numerical inversion of the Dirac operator. One of the simplest methods to solve such large and sparse linear systems is the conjugate gradient (CG) approach. In this work we present an implementation of CG that can be executed on different devices, including CPUs, GPUs, and FPGAs. This is achieved by using the SYCL/DPC++ framework, which allows the execution of the same source code on heterogeneous systems.