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

dim=1;
NoPotential[];
sizes={256};
timeStep=0.001;
k=-0.5;
contactInteractionFactor=4 Pi k;
initialize[];
gridReport

k0=0;
L=max[[1]]-min[[1]]+delta[[1]];
Δk=(2Pi/L)2;
a=0.01;
(* Initial WF: very small modulation with wavenumber Δk *)
psi0[{x_}]:=Exp[I k0 x](1+a Exp[I Δk x]+a Exp[-I Δk x]);
Plot[Abs[psi0[{x}]],{x,min[[1]],max[[1]]},PlotRange->{0,1.1}]
wf=Map[psi0,mesh[],{dim}];
wf/=Sqrt[norm];

calcAllTable
p1=plotProjections;

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

calcAllTable
p2=plotProjections;
Show[p1,p2,PlotRange->All,AxesOrigin->{-5,0}]

reportEnergies

reportMoments

reportMisc