QCD simulations with staggered fermions on GPUs

Abstract

We report on our implementation of the RHMC algorithm for the simulation of lattice QCD with two staggered flavors on Graphics Processing Units, using the NVIDIA CUDA programming language. The main feature of our code is that the GPU is not used just as an accelerator, but instead the whole Molecular Dynamics trajectory is performed on it. After pointing out the main bottlenecks and how to circumvent them, we discuss the obtained performances. We present some preliminary results regarding OpenCL and multiGPU extensions of our code and discuss future perspectives.

22 pages, 14 eps figures, final version to be published in Computer Physics Communications

References in corpus (8)

• Accelerating Dynamical Fermion Computations using the Rational Hybrid Monte Carlo (RHMC) Algorithm with Multiple Pseudofermion Fields
• Lattice QCD as a video game
• Accelerating Staggered Fermion Dynamics with the Rational Hybrid Monte Carlo (RHMC) Algorithm
• The Roberge-Weiss endpoint in N_f = 2 QCD
• Multi-mass solvers for lattice QCD on GPUs
• Topological Susceptibility in Two Flavors Lattice QCD with the Optimal Domain-Wall Fermion
• Parallelizing the QUDA Library for Multi-GPU Calculations in Lattice Quantum Chromodynamics
• Blasting through lattice calculations using CUDA

Cited by in corpus (20)

• Curvature of the pseudocritical line in QCD: Taylor expansion matches analytic continuation
• The chiral phase transition in two-flavor QCD from imaginary chemical potential
• Massively parallel lattice-Boltzmann codes on large GPU clusters
• OpenFPM: A scalable open framework for particle and particle-mesh codes on parallel computers
• Magnetic Susceptibility of Strongly Interacting Matter across the Deconfinement Transition
• A Framework for Lattice QCD Calculations on GPUs
• Susceptibility of the QCD vacuum to CP-odd electromagnetic background fields
• Lattice QCD based on OpenCL
• Coulomb, Landau and Maximally Abelian Gauge Fixing in Lattice QCD with Multi-GPUs
• Roberge-Weiss endpoint and chiral symmetry restoration in $N_f=2+1$ QCD
• Computational Physics on Graphics Processing Units
• Design and optimization of a portable LQCD Monte Carlo code using OpenACC
• Generating SU(Nc) pure gauge lattice QCD configurations on GPUs with CUDA
• Gauge-invariant field-strength correlators for QCD in a magnetic background
• Portable multi-node LQCD Monte Carlo simulations using OpenACC
• Gauge-invariant screening masses and static quark free energies in $N_f=2+1$ QCD at non-zero baryon density
• Dependence of the static quark free energy on ${\mu }_B$ and the crossover temperature of $N_f=2+1$ QCD
• QCDGPU: open-source package for Monte Carlo lattice simulations on OpenCL-compatible multi-GPU systems
• Phase structure of compactified $SU(N)$ gauge theories in magnetic backgrounds
• LatticeQCD using OpenCL