Solving the ghostgluon system of YangMills theory on GPUs
We solve the ghostgluon system of YangMills theory using graphics processing units (GPUs). Working in the Landau gauge, we use the DysonSchwinger formalism for the mathematical description, as this approach is well suited to directly benefit from the computing power of the GPUs. With the help of a Chebyshev expansion for the dressing functions and a subsequent appliance of a NewtonRaphson method, the non-linear system of coupled integral equations is linearized. The resulting Newton matrix is generated in parallel using OpenMPI and CUDA. Our results show that it is possible to cut down the run time by two orders of magnitude as compared to a sequential version of the code. This makes the proposed techniques well suited for DysonSchwinger calculations on more complicated systems where the YangMills sector of QCD serves as a starting point. In addition, the computation of Schwinger functions using GPU devices is studied.
ghostgluon system yangmills graphics processing units dysonschwinger equations parallel computing cudaprogramming model