(*

Eddy Ardonne and Nicolas Regnault

Squeeze Routines for polarized and spinful quantum Hall states,
bosonic as well as fermionic. Routines for the spinful case
are based on

The structure of spinful quantum Hall states: a squeezing perspective,
E. Ardonne, N. Regnault,
arXiv:1107.2232

See this paper for references for the polarized case, and for reference
about particle and orbital entanglement spectra.

Routines making use of the additional symmetries in case of L=0 and/or
S=0 states are not implemented in this version. A routine to calculate
only the number of solutions to the consistency condition via an
eigenvalue problem is not included either. Please contact E. Ardonne
about these (and possible other) issues.

*)

<<SqueezeRoutines`
<<ParticleEntanglement`
<<OrbitalEntanglement`

?SqueezeRoutines`*

?ParticleEntanglement`*

?OrbitalEntanglement`*

bosonsqueeze[{2,0,0,2,0,0,2,0,0,2}]
hsdim
bosonlplus
listlength
solveall

normalization="sphere";
bosonpes[4];

lzstates
lmultiplets[lzstates]
Plus@@lmultiplets[lzstates]

datal=datalsorted;

plot=ListPlot[data,Frame->True,Axes->False,PlotRange->{{Min[data[[All,1]]]-1,Max[data[[All,1]]]+1},{Min[data[[All,2]]]-1/2,Max[data[[All,2]]]+1/2}},PlotStyle->{Red}]

plot=ListPlot[datal,Frame->True,Axes->False,PlotRange->{{Min[datal[[All,1]]]-1,Max[datal[[All,1]]]+1},{Min[datal[[All,2]]]-1/2,Max[datal[[All,2]]]+1/2}},PlotStyle->{Red}]

normalization="sphere";
oes[5]
nvalues

range=nvalues;
data=Flatten[Table[ndata[i],{i,range}],1];
plot=ListPlot[data,Frame->True,Axes->False,PlotStyle->Red,PlotRange->{{Min[data[[All,1]]]-1,Max[data[[All,1]]]+1},{Floor[Min[data[[All,2]]]]-1/2,Ceiling[Max[data[[All,2]]]]+1/2}}]

normalization="cft";
bosonpes[4];

plot=ListPlot[data,Frame->True,Axes->False,PlotRange->{{Min[data[[All,1]]]-1,Max[data[[All,1]]]+1},{Min[data[[All,2]]]-1/2,Max[data[[All,2]]]+1/2}},PlotStyle->{Red}]

bosonsqueeze[{2,0,0,2,0,0,2}]
hsdim
bosonlplus
listlength
solveall

b[1]=4;
Table[b[i],{i,1,hsdim}]

pol=bosonpol;

Factor[pol/.{z[2]->z[1]}]

spinbosonsqueeze[{2,2,0,2,2},4,4]
hsdim
spinbosonlplus
spinbosonsplus
listlength
solveall

b[4]=8;
Table[b[i],{i,1,hsdim}]

normalization="sphere";
spinoes[2]

normalization="sphere";
spinbosonpes[4];

plus={Line[{{0,-1/2},{0,1/2}}],Line[{{-1/2,0},{1/2,0}}]};
cross={Line[{{-1/2,-1/2},{1/2,1/2}}],Line[{{1/2,-1/2},{-1/2,1/2}}]};
square=Line[{{-1/2,-1/2},{-1/2,1/2},{1/2,1/2},{1/2,-1/2},{-1/2,-1/2}}];
dot=Disk[{0,0},.05];

plot=ListPlot[data,Frame->True,Axes->False,PlotRange->{{Min[data[[All,All,1]]]-1,Max[data[[All,All,1]]]+1},{Min[data[[All,All,2]]]-1/2,Max[data[[All,All,2]]]+1/2}},PlotMarkers->{{Graphics[{Directive[Red],cross}],0.03},{Graphics[{Directive[Blue],plus}],0.03},{Graphics[{Directive[Green],square}],0.03}}]

datals=datalssorted;

plot=ListPlot[datals,Frame->True,Axes->False,PlotRange->{{Min[datals[[All,All,1]]]-1,Max[datals[[All,All,1]]]+1},{Min[datals[[All,All,2]]]-1/2,Max[datals[[All,All,2]]]+1/2}},PlotMarkers->{{Graphics[{Directive[Red],cross}],0.03},{Graphics[{Directive[Blue],plus}],0.03},{Graphics[{Directive[Green],square}],0.03}}]

lsmultiplets[lzszstates]//MatrixForm