C++ CUDA programs for solving the time-dependent dipolar Gross-Pitaevskii equation

This paper introduces an improved version of our previous CUDA programs [1] for solving the dipolar Gross-Pitaevskii equation in three spatial dimensions, now incorporating a quantum fluctuation term [2, 3, 4]. This enhancement is vital for accurately modeling quantum dipolar droplets [5], [6] in dipolar Bose-Einstein condensates. Originally developed in C, the code has been transitioned to C++ to take advantage of its features. Speedup tests were conducted on two types of GPU cards – one commercial and one HPC-optimized – and the results were compared to earlier versions on GPUs and a modern HPC cluster. Both configurations showed an increase in speed.