GSGPEs-v1.1: A MATLAB code for computing the ground state of systems of Gross–Pitaevskii equations

It can compute the ground state in the defocusing case, for any number of equations with harmonic or quasi-harmonic trapping potentials, in spatial dimension one, two or three. The computation is based on a spectral decomposition of the solution into Hermite functions and direct minimization of the energy functional through a Newton-like method with an approximate line-search strategy.

This new version is due to a change in the function eig of Matlab® which requires a new way to compute Gauss–Hermite quadrature nodes and weights. The previous version of this program (AENT_v1_0) may be found at dx.doi.org/10.1016/j.cpc.2012.10.007.