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

HarmonicTrap3D[0.5,0.7,1.0];
sizes={32,32,32};
timeStep=0.0005;
k=-0.5;
contactInteractionFactor=4 Pi k;
initialize[];
gridReport

calcAllTable
p1=plotProjections

AbsoluteTiming[evolve["ite",15,10]]

calcAllTable
p2=plotProjections;
Show[p1,p2,PlotRange->All]

reportEnergies

reportMoments