Hi,all
I am trying to sparsediag the Kondo lattice in ALPS, but the command line return "Called get_quantumnumber with illegal index". I feel that there may be a misunderstanding of the constraint in my code, but I just can not find it out now. So, dear users and developers please help me!!!

The parameter file is,

LATTICE="chain lattice"
MODEL="Kondo lattice"
LATTICE_LIBRARY="mylattices.xml"
MODEL_LIBRARY="mymodels.xml"
CONSERVED_QUANTUMNUMBERS="Sz"
Sz_total=0
t=1
J=2
{L=4}

the "Kondo lattice" code in mymodels.xml is almostly copied from the "Kondo lattice" code in ALPS models.xml, the "only" change I made (no changes of sitebasis, operators or hamiltonian) is adding a constraint of the total Sz (Sz_spin+Sz_fermion) to the BASIS as shown below:

<BASIS name="Kondo lattice">
  <SITEBASIS type="0" ref="fermion"/>
  <SITEBASIS type="1" ref="spin"/>
  <CONSTRAINT quantumnumber="Sz" value="Sz_total"/>
</BASIS>

The lattice in mylattices.xml is defined as:

<LATTICES>

<LATTICE name="chain lattice" dimension="1">
  <PARAMETER name="a" default="1"/>
  <BASIS><VECTOR>a</VECTOR></BASIS>
  <RECIPROCALBASIS><VECTOR>2*pi/a</VECTOR></RECIPROCALBASIS>
</LATTICE>
<UNITCELL name="hf" dimension="1">
  <VERTEX type="0"/>
  <VERTEX type="1"/>
  <EDGE type="0"><SOURCE vertex="1" offset="0"/><TARGET vertex="1" offset="1"/></EDGE>
  <EDGE type="1"><SOURCE vertex="1" offset="0"/><TARGET vertex="2" offset="0"/></EDGE>
</UNITCELL>
<LATTICEGRAPH name = "chain lattice">
  <FINITELATTICE>
    <LATTICE ref="chain lattice"/>
    <PARAMETER name="L" default="4"/>
    <EXTENT size="L"/>
    <BOUNDARY type="periodic"/>
  </FINITELATTICE>
  <UNITCELL ref="hf"/>
</LATTICEGRAPH>

</LATTICES>

where VERTEX type=0 stands for fermion and VERTEX type=1 stands for spin.

Any advice is welcome!

Bo-Nan
--
Stay foolish,Stay hungry.