OpenMP Fortran and C programs for solving the time-dependent Gross-Pitaevskii equation in an anisotropic trap

TL;DR

This paper introduces optimized OpenMP-parallelized Fortran and C programs for solving the Gross-Pitaevskii equation in multiple dimensions, with inputs and outputs aligned with experimental parameters, improving computational efficiency and usability.

Contribution

The paper presents new parallelized versions of existing programs with inputs based on experimental parameters and includes performance benchmarks.

Findings

Significant speedup in execution time on multicore processors.

Programs now accept experimentally relevant input parameters.

New output formats for density distributions.

Abstract

We present new version of previously published Fortran and C programs for solving the Gross-Pitaevskii equation for a Bose-Einstein condensate with contact interaction in one, two and three spatial dimensions in imaginary and real time, yielding both stationary and non-stationary solutions. To reduce the execution time on multicore processors, new versions of parallelized programs are developed using Open Multi-Processing (OpenMP) interface. The input in the previous versions of programs was the mathematical quantity nonlinearity for dimensionless form of Gross-Pitaevskii equation, whereas in the present programs the inputs are quantities of experimental interest, such as, number of atoms, scattering length, oscillator length for the trap, etc. New output files for some integrated one- and two-dimensional densities of experimental interest are given. We also present speedup test results for the new programs.