Efficient single-precision simulations of nematohydrodynamics

Simulations of nematohydrodynamics on graphics processing units (GPUs) are typically performed using double precision, which ensures accuracy but significantly increases computational cost. However, consumer-grade GPUs are optimized for single-precision calculations, making double-precision simulations inefficient on widely available hardware. In this work, we demonstrate that single-precision simulations can achieve the same accuracy as double-precision methods while delivering a 27-fold increase in the computational speed. To achieve this, we introduce two key improvements: (i) the shifted distribution function in the lattice Boltzmann method, which mitigates precision loss at low velocities, and (ii) the use of larger time steps in the finite-difference solver, which reduces numerical errors and improves the overall accuracy. We find that, unlike in double precision, accuracy in single-precision simulations follows a non-monotonic trend with respect to the finite-difference time step, revealing an optimal regime for precise computations. To illustrate the effectiveness of our approach, we simulate the dynamics of single and multiple skyrmionic tubes in Poiseuille flow. Our results confirm that optimized single-precision simulations enable fast and accurate modeling of complex nematohydrodynamic systems, making large-scale simulations feasible on standard gaming GPUs.