NotebookEvaluate[NotebookDirectory[]<>"/split_step.nb"];

HarmonicTrap3D[];
sizes={16,16,16};
timeStep=0.001;
k=0;
contactInteractionFactor=4 Pi k;
initialize[];
gridReport

calcAllTable
plotProjections

AbsoluteTiming[evolve["ite",2.5,5]]

calcAllTable
plotProjections

reportEnergies

reportMoments