<div dir="ltr"><div class="gmail_default" style="font-size:small">I tried to follow your instructions, but I am doing something wrong and cannot find out what it is. In order to test creating a graph with different coupling types, I created a file called 'my_lattices.xml' implementing a 9-site dimerized chain with the following contents:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><div class="gmail_default"><font face="monospace, monospace"><LATTICES></font></div><div class="gmail_default"><font face="monospace, monospace"> <GRAPH name = "open chain dimerized lattice" dimension="1" vertices="9" edges="8"></font></div><div class="gmail_default"><font face="monospace, monospace">  <VERTEX id="1"/></font></div><div class="gmail_default"><font face="monospace, monospace">  <VERTEX id="2"/></font></div><div class="gmail_default"><font face="monospace, monospace">  <VERTEX id="3"/></font></div><div class="gmail_default"><font face="monospace, monospace">  <VERTEX id="4"/></font></div><div class="gmail_default"><font face="monospace, monospace">  <VERTEX id="5"/></font></div><div class="gmail_default"><font face="monospace, monospace">  <VERTEX id="6"/></font></div><div class="gmail_default"><font face="monospace, monospace">  <VERTEX id="7"/></font></div><div class="gmail_default"><font face="monospace, monospace">  <VERTEX id="8"/></font></div><div class="gmail_default"><font face="monospace, monospace">  <VERTEX id="9"/></font></div><div class="gmail_default"><font face="monospace, monospace">  <EDGE source="1" target="2" id="1" type="0"/></font></div><div class="gmail_default"><font face="monospace, monospace">  <EDGE source="2" target="3" id="2" type="1"/></font></div><div class="gmail_default"><font face="monospace, monospace">  <EDGE source="3" target="4" id="3" type="0"/></font></div><div class="gmail_default"><font face="monospace, monospace">  <EDGE source="4" target="5" id="4" type="1"/></font></div><div class="gmail_default"><font face="monospace, monospace">  <EDGE source="5" target="6" id="5" type="0"/></font></div><div class="gmail_default"><font face="monospace, monospace">  <EDGE source="6" target="7" id="6" type="1"/></font></div><div class="gmail_default"><font face="monospace, monospace">  <EDGE source="7" target="8" id="7" type="0"/></font></div><div class="gmail_default"><font face="monospace, monospace">  <EDGE source="8" target="9" id="8" type="1"/></font></div><div class="gmail_default"><font face="monospace, monospace"> </GRAPH></font></div><div class="gmail_default"><font face="monospace, monospace"></LATTICES></font></div></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I then adapted a Python file from the tutorials in order to calculate the susceptibility of an antiferromagnetic Heisenberg model sitting on the sites of the dimerized lattice, using 'J0' and 'J1' to refer to the couplings associated with edge types 0 and 1, which I set equal to each other for testing: </div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><div class="gmail_default"><font face="monospace, monospace">import pyalps</font></div><div class="gmail_default"><font face="monospace, monospace">import matplotlib.pyplot as plt</font></div><div class="gmail_default"><font face="monospace, monospace">import pyalps.plot</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">#prepare the input parameters</font></div><div class="gmail_default"><font face="monospace, monospace">#skip this part if you already ran the simulation from the command line</font></div><div class="gmail_default"><font face="monospace, monospace">parms = []</font></div><div class="gmail_default"><font face="monospace, monospace">for t in [0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.25, 1.5, 1.75, 2.0]:</font></div><div class="gmail_default"><font face="monospace, monospace">    parms.append(</font></div><div class="gmail_default"><font face="monospace, monospace">        { </font></div><div class="gmail_default"><font face="monospace, monospace">          'LATTICE_LIBRARY': 'my_lattice.xml',</font></div><div class="gmail_default"><font face="monospace, monospace">          'LATTICE'        : 'open chain dimerized lattice', </font></div><div class="gmail_default"><font face="monospace, monospace">          'T'              : t,</font></div><div class="gmail_default"><font face="monospace, monospace">          'local_S'        : 0.5,</font></div><div class="gmail_default"><font face="monospace, monospace">          'J0'             : 1.0,</font></div><div class="gmail_default"><font face="monospace, monospace">          'J1'             : 1.0,</font></div><div class="gmail_default"><font face="monospace, monospace">          'THERMALIZATION' : 10000,</font></div><div class="gmail_default"><font face="monospace, monospace">          'SWEEPS'         : 500000,</font></div><div class="gmail_default"><font face="monospace, monospace">          'ALGORITHM'      : 'loop',</font></div><div class="gmail_default"><font face="monospace, monospace">          'MODEL'          : 'spin'</font></div><div class="gmail_default"><font face="monospace, monospace">        }</font></div><div class="gmail_default"><font face="monospace, monospace">    )</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">#write the input file and run the simulation</font></div><div class="gmail_default"><font face="monospace, monospace">input_file = pyalps.writeInputFiles('parm2a',parms)</font></div><div class="gmail_default"><font face="monospace, monospace">pyalps.runApplication('loop',input_file,Tmin=5)</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">#load the susceptibility and collect it as function of temperature T</font></div><div class="gmail_default"><font face="monospace, monospace">data = pyalps.loadMeasurements(pyalps.getResultFiles(prefix='parm2a'),'Susceptibility')</font></div><div class="gmail_default"><font face="monospace, monospace">susceptibility = pyalps.collectXY(data,x='T',y='Susceptibility')</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">#make plot</font></div><div class="gmail_default"><font face="monospace, monospace">plt.figure()</font></div><div class="gmail_default"><font face="monospace, monospace">plt.title("Antiferromagnetic Heisenberg chain")</font></div><div class="gmail_default"><font face="monospace, monospace">pyalps.plot.plot(susceptibility)</font></div><div class="gmail_default"><font face="monospace, monospace">plt.show()</font></div><div class="gmail_default"><br></div><div class="gmail_default">When I run the code with alpspython, I get a segmentation fault, but if I use the standard lattice library and an 'open chain lattice' with 9 sites, everything works as it should. </div><div class="gmail_default"><br></div><div class="gmail_default">Is there anything obviously silly in what I am doing?</div></div></div><div class="gmail_extra"><br><div class="gmail_quote">2017-08-02 19:14 GMT-03:00 Matthias Troyer <span dir="ltr"><<a href="mailto:troyer@phys.ethz.ch" target="_blank">troyer@phys.ethz.ch</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">You need to create a graph with different coupling types and specify all couplings individually.<div><br><div><blockquote type="cite"><div><div class="h5"><div>On 2 Aug 2017, at 13:45, Andre Vieira <<a href="mailto:apvieira@if.usp.br" target="_blank">apvieira@if.usp.br</a>> wrote:</div><br class="m_910827451660558412Apple-interchange-newline"></div></div><div><div><div class="h5"><div dir="ltr"><div class="gmail_default" style="font-size:small">Hi,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">is there a way of using the 'inhomogeneous chain lattice' definition to implement an arbitrary nonperiodic but deterministic lattice for QMC applications? To be specific, I would like to implement a spin model in which, for instance, J(1), J(2), J(3), etc. are all different, but not defined by a simple function, nor randomly chosen.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Thanks in advance for your help,</div><div class="gmail_default" style="font-size:small">Andre.</div></div>
<br><br></div></div>----<br>Comp-phys-alps-users Mailing List for the ALPS Project<br><a href="http://alps.comp-phys.org/" target="_blank">http://alps.comp-phys.org/</a><br><br>List info: <a href="https://lists.phys.ethz.ch//listinfo/comp-phys-alps-users" target="_blank">https://lists.phys.ethz.ch//<wbr>listinfo/comp-phys-alps-users</a><br>Archive: <a href="https://lists.phys.ethz.ch//pipermail/comp-phys-alps-users" target="_blank">https://lists.phys.ethz.ch//<wbr>pipermail/comp-phys-alps-users</a><br><br>Unsubscribe by writing a mail to <a href="mailto:comp-phys-alps-users-leave@lists.phys.ethz.ch" target="_blank">comp-phys-alps-users-leave@<wbr>lists.phys.ethz.ch</a>.</div></blockquote></div><br></div></div><br><br>
<br>
----<br>
Comp-phys-alps-users Mailing List for the ALPS Project<br>
<a href="http://alps.comp-phys.org/" rel="noreferrer" target="_blank">http://alps.comp-phys.org/</a><br>
<br>
List info: <a href="https://lists.phys.ethz.ch//listinfo/comp-phys-alps-users" rel="noreferrer" target="_blank">https://lists.phys.ethz.ch//<wbr>listinfo/comp-phys-alps-users</a><br>
Archive: <a href="https://lists.phys.ethz.ch//pipermail/comp-phys-alps-users" rel="noreferrer" target="_blank">https://lists.phys.ethz.ch//<wbr>pipermail/comp-phys-alps-users</a><br>
<br>
Unsubscribe by writing a mail to <a href="mailto:comp-phys-alps-users-leave@lists.phys.ethz.ch">comp-phys-alps-users-leave@<wbr>lists.phys.ethz.ch</a>.<br></blockquote></div><br></div>