We present an effective adaptive procedure for the numerical approximation of
the steady-state Gross-Pitaevskii equation. Our approach is solely based on
energy minimization, and consists of a combination of gradient flow iterations
and adaptive finite element mesh refinements. Numerical tests show that this
strategy is able to provide highly accurate results, with optimal convergence
rates with respect to the number of freedom
