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

HarmonicTrap1D[];
sizes={128};
timeStep=0.001;
k=0;
contactInteractionFactor=4 Pi k;
initialize[];
gridReport

calcAllTable
p1=plotProjections

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

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

reportEnergies

reportMoments

reportMisc