We investigate the use of half-precision floating-point numbers (FP16) in mixed-precision linear solvers for lattice QCD simulations. Since the emergence of GPUs for general-purpose, mixed-precision algorithms that combine single-precision (FP32) with double-precision (FP64) arithmetics have become widely used in this field and others. While FP32-based methods are now well established, we examine the practicality of using FP16. In this work, we introduce rescaling steps in both the outer iterative refinement step and the inner BiCGStab solver to avoid numerical instability. In our experiments with a simple Wilson kernel, the solver shows improved stability, and the additional iteration count compared to the FP64 version remains within 20\%, indicating that the FP16 version is practical for use. We believe that the proposed rescaling methods can also benefit other mixed precision preconditioners in avoiding underflows.
@article{arxiv.2602.14450,
title = {Mixed precision solvers with half-precision floating point numbers for Lattice QCD on A64FX processor},
author = {Issaku Kanamori and Hideo Matsufuru and Tatsumi Aoyama and Kazuyuki Kanaya and Yusuke Namekawa and Hidekatsu Nemura and Keigo Nitadori},
journal= {arXiv preprint arXiv:2602.14450},
year = {2026}
}
Comments
11 pages, 9 figures, contribution to the International Workshop on Arm-based HPC: Practice and Experience 2026 (IWAHPCE2026) at HPC-Asia 2026