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

HarmonicTrap2D[];
sizes={64,64};
timeStep=0.001;
k=-0.3;
contactInteractionFactor=4 Pi k;
initialize[];
gridReport

calcAllTable
p1=plotProjections

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

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

reportEnergies

reportMoments