DMRG, different type of sites - Comp-phys-alps-users

7 May 2009


      Hello all,
This is my first post on this mailing list so I would like say hello to all
community members. Great thanks to developers for they work on such nice
thing like ALPS.
I'm currently using alps version 1.34 with boost version 1.34.1. I have
following problem.
I defined simple graph in lattices.xml file:
<GRAPH name="circle1d">
  <VERTEX id="1" type="0"/>
  <VERTEX id="2" type="1"/>
  <VERTEX id="3" type="2"/>
  <VERTEX id="4" type="3"/>
  <VERTEX id="5" type="4"/>
  <VERTEX id="6" type="5"/>
  <VERTEX id="7" type="6"/>
  <VERTEX id="8" type="7"/>
  <VERTEX id="9" type="8"/>
  <VERTEX id="10" type="9"/>
  <EDGE source="1" target="2" type="0"/>
  <EDGE source="2" target="3" type="0"/>
  <EDGE source="3" target="4" type="0"/>
  <EDGE source="4" target="5" type="0"/>
  <EDGE source="5" target="6" type="0"/>
  <EDGE source="6" target="7" type="0"/>
  <EDGE source="7" target="8" type="0"/>
  <EDGE source="8" target="9" type="0"/>
  <EDGE source="9" target="10" type="0"/>
  <EDGE source="10" target="1" type="0"/>
</GRAPH>
next I create parameter file in form:
LATTICE="circle1d";
SWEEPS=4;
MAXSTATES=100;
MODEL="spin";
NUMBER_EIGENVALUES=2;
J=1;
MEASURE_AVERAGE[Magnetization]=Sz;
MEASURE_AVERAGE[Exchange]=exchange;
MEASURE_LOCAL[Local magnetization]=Sz;
MEASURE_CORRELATIONS[Diagonal spin correlations]=Sz;
MEASURE_CORRELATIONS[Offdiagonal spin correlations]="Splus:Sminus";
{ h0=cos(0.0); h1=cos(2.0*pi/10); h2=cos(4.0*pi/10); h3=cos(6.0*pi/10);
h4=cos(8.0*pi/10); h5=cos(10.0*pi/10); h6=cos(12.0*pi/10);
h7=cos(14.0*pi/10); h8=cos(16.0*pi/10); h9=cos(18.0*pi/10); Gamma0=sin(0.0);
Gamma1=sin(2.0*pi/10); Gamma2=sin(4.0*pi/10); Gamma3=sin(6.0*pi/10);
Gamma4=sin(8.0*pi/10); Gamma5=sin(10.0*pi/10); Gamma6=sin(12.0*pi/10);
Gamma7=sin(14.0*pi/10); Gamma8=sin(16.0*pi/10); Gamma9=sin(18.0*pi/10); }
When I started this task I get "segmentation fault (core dumped)" error. I
notice that when using graph with less then 5 node types I don't have any
problem with this.
What I do wrong? Is there any restrictions for number of node types ?
Great Thanks,
Remigiusz Ziemkiewicz