Split-step Fourier methods for the Gross-Pitaevskii equation

We perform a systematic study of the accuracy of split-step Fourier transform methods for the time dependent Gross-Pitaevskii equation using symbolic calculation. Provided the most recent approximation for the wave function is always used in the nonlinear atom-atom interaction potential energy, every split-step algorithm we have tried has the same-order time stepping error for the Gross-Pitaevskii equation and the Schroedinger equation.