Implementation of the twisted mass fermion operator in the QUDA library

TL;DR

This paper extends the QUDA library to efficiently implement the Wilson twisted mass fermion operator, demonstrating high-performance results on NVIDIA GPUs for both degenerate and non-degenerate flavor doublets, supporting hadron structure research.

Contribution

The paper introduces a new implementation of the twisted mass fermion operator in QUDA, optimizing performance for GPU computing in lattice QCD simulations.

Findings

Achieves up to 856 Gflops for degenerate doublets at half precision.

Reaches 879 Gflops for non-degenerate doublets at half precision.

Supports production-level calculations for hadron structure studies.

Abstract

We discuss an extension of the QUDA library for the Wilson twisted mass operator. A performance analysis is presented for both degenerate and non-degenerate flavor doublets. The degenerate twisted mass fermion operator runs at up to 190, 487 and 856 Gflops, for double, single and half precisions respectively on recent NVIDIA Kepler GPUs, while our implementation for the non-degenerate flavor doublet allows to reach 163, 516 and 879 GFlops, respectively. The code is currently in production for the hadron structure study.