CUDA programs for solving the time-dependent dipolar Gross-Pitaevskii equation in an anisotropic trap

• Published in: arXiv.org
• Main Authors: Loncar, Vladimir, Balaz, Antun, Bogojevic, Aleksandar, Skrbic, Srdjan, Muruganandam, Paulsamy, Adhikari, Sadhan K
• Format: Paper Journal Article
• Language: English
• Published: Ithaca Cornell University Library, arXiv.org 01.08.2022
• Subjects: Algorithms; Fast Fourier transformations; Fourier transforms; Physics - Computational Physics; Physics - Pattern Formation and Solitons; Physics - Quantum Gases; Physics - Quantum Physics; Time dependence
• ISSN: 2331-8422
• DOI: 10.48550/arxiv.1601.04640

In this paper we present new versions of previously published numerical programs for solving the dipolar Gross-Pitaevskii (GP) equation including the contact interaction in two and three spatial dimensions in imaginary and in real time, yielding both stationary and non-stationary solutions. New versions of programs were developed using CUDA toolkit and can make use of Nvidia GPU devices. The algorithm used is the same split-step semi-implicit Crank-Nicolson method as in the previous version (R. Kishor Kumar et al., Comput. Phys. Commun. 195, 117 (2015)), which is here implemented as a series of CUDA kernels that compute the solution on the GPU. In addition, the Fast Fourier Transform (FFT) library used in the previous version is replaced by cuFFT library, which works on CUDA-enabled GPUs. We present speedup test results obtained using new versions of programs and demonstrate an average speedup of 12 to 25, depending on the program and input size.