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

HarmonicTrap1D[];
sizes={1024};
timeStep=0.001;
k=-0.4;
contactInteractionFactor=4 Pi k;
initialize[];
gridReport

calcAllTable
p1=plotProjections

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

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

reportMoments

Turn off harmonic potential

newPotential[(0&)];
plot:=AppendTo[lplots,{time,plotProjections}];
reset[];
AbsoluteTiming[evolve["ite",5,20,20]]

ListPlot[lmoments[[All,{1,3}]],Joined->True,PlotRange->All]

Show[lplots[[All,2]],PlotRange->All]