- Comp-phys-alps-users - lists.phys.ethz.ch

TEBD application: How to access the Schmidt eigenvalues ?
by Grégoire Misguich 24 Feb '11

24 Feb '11

Dear alps users, I would like to know if one can simply have access the to the Schmidt eigenvalues which describe the wave-function in (or at the end of) the TEBD application. It there a "measurement" for that ? Thank you ! Gregoire Misguich

1 0

continue or extend a run using directloop_sse
by Kubila 18 Feb '11

18 Feb '11

Hi ALPS 2.0 users I was wondering did anyone succeed in extending or continuing a run using directloop_sse in ALPS 2.0? Using job.out.xml seems to work with loop (./loop job.out.xml works), but not with dirctloop_sse (./directloop_sse job.out.xml fails).... Any help will appreciate very much! Sincerely, Kubila

2 1

quantum XY model with an uniform magnetic field in the x-direction
by Kubila 11 Feb '11

11 Feb '11

Hi ALPS users I was wondering can one use looper to simulate quantum XY model with an uniform magnetic field in the x-direction. Thanks.

3 6

TEBD tutorial "2a" (kink & XX chain)
by Grégoire Misguich 09 Feb '11

09 Feb '11

Dear ALPS community, As my first attempt to use the TEBD code provided in alps-2.0, I ran a tutorial concerning the simulation of the XX spin chain (tutorial_2a.py), which simulates the dynamics starting from a domain wall spin configuration: http://alps.comp-phys.org/mediawiki/index.php/ALPS_2_Tutorials:TEBD-02_kink I got the following figures (attached .png files). The Wiki says "We see that the magnetization agrees very well to visual accuracy", but in the plot I obtained, it is not quite the case (very "noisy"). I played a bit with some parameters (increasing the number of time steps and the max Schmidt dimension), but it still get some ugly curves. So, I suspect that something might we wrong... Any help or suggestion would be very much appreciated ! Gregoire Misguich.

1 1

RESTRICT_MEASUREMENTS[Sz] dirloop_sse alps 1.35
by Marcelo Guimarães 08 Feb '11

08 Feb '11

Hello, I would like to know if the parameter RESTRICT_MEASUREMENTS[Sz] works for dirloop_sse as indicated at http://alps.comp-phys.org/mediawiki/index.php/Documentation:dirloop_sse. I'm trying to obtain the gap for (3d), H=J( S+(i)S(i+1)- + S-(i)S+(i+1) ) +D(Sz)^2 at D > Dc by calculating energy difference from E(Mz=1) , E(Mz=0) at low temperature. Since <Sz(i)Sz(i+1)> ~ 0 , and dirloop does not measure off-diagonal correlations , I was wandering if I could use this procedure in dirloop to estimate the gap. Thanks, Marcelo Guimaraes.

2 1

Re: [ALPS-users] Extracting dates
by Xiao Quan Yu 07 Feb '11

07 Feb '11

Dear Prof. Troyer, Thank you for your replying. One of my out put file is the following: <?xml version="1.0" encoding="UTF-8"?> <?xml-stylesheet type="text/xsl" href="file:///u/he/warter/ALPS/lib/xml/ALPS.xsl"?> <SIMULATION xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://xml.comp-phys.org/2002/10/ALPS.xsd"> <PARAMETERS> <PARAMETER name="LATTICE_LIBRARY">lattices.xml</PARAMETER> <PARAMETER name="LATTICE">inhomogeneous simple cubic lattice</PARAMETER> <PARAMETER name="L">5</PARAMETER> <PARAMETER name="MODEL_LIBRARY">models.xml</PARAMETER> <PARAMETER name="MODEL">hardcore boson</PARAMETER> <PARAMETER name="t">1.0</PARAMETER> <PARAMETER name="T">0.1</PARAMETER> <PARAMETER name="delta">10.0</PARAMETER> <PARAMETER name="SWEEPS">500000</PARAMETER> <PARAMETER name="THERMALIZATION">10000</PARAMETER> <PARAMETER name="EPSILON">0.25</PARAMETER> <PARAMETER name="DISORDERSEED">65897</PARAMETER> <PARAMETER name="mu">delta*2*(random()-0.5)</PARAMETER> <PARAMETER name="SEED">668368703</PARAMETER> </PARAMETERS> <AVERAGES> <SCALAR_AVERAGE name="Density"> <COUNT>497664</COUNT> <MEAN method="jackknife">0.52646836</MEAN> <ERROR converged="yes" method="jackknife">2.09e-05</ERROR> <VARIANCE method="simple">5.45e-05</VARIANCE> <AUTOCORR method="jackknife">1.5</AUTOCORR> </SCALAR_AVERAGE> <SCALAR_AVERAGE name="Density^2"> <COUNT>497664</COUNT> <MEAN method="jackknife">34.652926</MEAN> <ERROR converged="yes" method="jackknife">0.00276</ERROR> <VARIANCE method="simple">0.944</VARIANCE> <AUTOCORR method="jackknife">1.5</AUTOCORR> </SCALAR_AVERAGE> <SCALAR_AVERAGE name="Energy"> <COUNT>497664</COUNT> <MEAN method="jackknife">-251.9492</MEAN> <ERROR converged="yes" method="jackknife">0.0221</ERROR> <VARIANCE method="simple">64.1</VARIANCE> <AUTOCORR method="jackknife">1.4</AUTOCORR> </SCALAR_AVERAGE> <SCALAR_AVERAGE name="Energy Density"> <COUNT>497664</COUNT> <MEAN method="jackknife">-2.0155936</MEAN> <ERROR converged="yes" method="jackknife">0.000177</ERROR> <VARIANCE method="simple">0.0041</VARIANCE> <AUTOCORR method="jackknife">1.4</AUTOCORR> </SCALAR_AVERAGE> <SCALAR_AVERAGE name="Stiffness"> <COUNT>497664</COUNT> <MEAN method="jackknife">2.17262</MEAN> <ERROR converged="yes" method="jackknife">0.00305</ERROR> <VARIANCE method="simple">3.16</VARIANCE> 1,1 Top ..................... In the mean time, I will try to use ALPS 2.0. Thank you. Best, Xiaoquan > Today's Topics: > > 1. Re: Extracting dates (Matthias Troyer) > 2. Re: Dzyaloshinskii-Moriya in spinmc (Matthias Troyer) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sun, 6 Feb 2011 23:03:46 +0200 > From: Matthias Troyer <troyer(a)phys.ethz.ch> > Subject: Re: [ALPS-users] Extracting dates > To: comp-phys-alps-users(a)lists.phys.ethz.ch > Message-ID: <ACE5D63D-7EEA-43D9-BDAC-862222FFEBD5(a)phys.ethz.ch> > Content-Type: text/plain; charset=us-ascii > > Can you send your output files? The problem is probably that mu is not > just a number. You might also look into using ALPS 2 instead which has > much better plotting and evaluation tools. > > Matthias > > On Feb 4, 2011, at 12:33 PM, Xiao Quan Yu wrote: > >> Dear ALPS community, >> >> I met a trouble for extracting the calculated results from the >> XML output files. My setting is the following: >> >> plot2.xml file: >> >> <?xml version="1.0" encoding="UTF-8"?> >> >> <plot name="Stiffness versus disorder for three-dimesnional >> HardcoreBoson >> models"> >> >> <legend show="true"/> >> <xaxis label="Disorder" type="PARAMETER" name="mu"/> >> <yaxis label="Stiffness" type="SCALAR_AVERAGE"/> >> >> <set label="Inhomogeneous simple cubic lattice"/> >> >> </plot> >> >> and my "input" file is >> >> LATTICE_LIBRARY="lattices.xml"; >> LATTICE="inhomogeneous simple cubic lattice"; >> L=5; >> MODEL_LIBRARY="models.xml"; >> MODEL="hardcore boson"; >> t = 1.0; >> T = 0.1; >> delta=10.0; >> SWEEPS=500000; >> THERMALIZATION=10000; >> EPSILON=0.25; >> { DISORDERSEED = 34275; mu=delta*2*(random()-0.5); } >> { DISORDERSEED = 49802; mu=delta*2*(random()-0.5); } >> { DISORDERSEED = 82529; mu=delta*2*(random()-0.5); } >> { DISORDERSEED = 13564; mu=delta*2*(random()-0.5); } >> { DISORDERSEED = 65897; mu=delta*2*(random()-0.5); }. >> >> When I typed the command "extracttext plot2.xml input.task*.out.xml", >> the >> following thing happend: >> >> >> plot2.xml:4: parser error : XML declaration allowed only at the start of >> the document >> <?xml version="1.0" encoding="UTF-8"?> >> ^ >> unable to parse plot2.xml >> /tmp/archive2plot3652.xsl:1: parser error : Document is empty >> >> ^ >> /tmp/archive2plot3652.xsl:1: parser error : Start tag expected, '<' not >> found >> >> ^ >> cannot parse /tmp/archive2plot3652.xsl >> /tmp/plot3652:1: parser error : Document is empty >> >> ^ >> /tmp/plot3652:1: parser error : Start tag expected, '<' not found >> >> ^ >> unable to parse /tmp/plot3652 >> >> >> I don't know what is going on here, could you help me figure it out? >> Thank >> you in advance! >> >> >> Best, >> Xiaoquan >> >> >> >> ---------------------------------------------------------------- >> SISSA Webmail https://webmail.sissa.it/ >> Powered by SquirrelMail http://www.squirrelmail.org/ >> > > > > ------------------------------ > > Message: 2 > Date: Sun, 6 Feb 2011 23:11:48 +0200 > From: Matthias Troyer <troyer(a)phys.ethz.ch> > Subject: Re: [ALPS-users] Dzyaloshinskii-Moriya in spinmc > To: comp-phys-alps-users(a)lists.phys.ethz.ch > Message-ID: <574253B4-92C0-4A0F-97AB-0ADF062B8A7A(a)phys.ethz.ch> > Content-Type: text/plain; charset=iso-8859-1 > > Dear Francesco, > > The source code needs to be modified, but it is very easy to write such a > code. Anyone in my group at ETH can teach you how to do that. I'm on > vacation and without good internet this week but we can meet next week, or > you can ask others on my group whether they can help you. Writing a local > update code for your model should not take more than half a day at the > most. > > Matthias > > > > > On Feb 2, 2011, at 6:32 PM, Francesco Casola wrote: > >> Dear ALPS community, >> >> could you please let me know about the possibility of including >> Dzyaloshinskii-Moriya in spinmc, possibly through the presence of an >> off-diagonal exchange coupling constant? >> >> Is it necessary to modify the source code for such a case, or can I >> manage it in some other way not explicitly mentioned in the ALPS >> tutorial? >> >> Best wishes >> Francesco >> >> -- >> Francesco Casola, Ph. D. Student >> Laboratory for Solid State Physics >> ETH Z?rich H?nggerberg >> HPF E15.1 Schafmattstrasse 16 >> CH-8093 Switzerland >> >> Tel. Office: +41 (0)44 633 22 96 / 03 >> Tel. Lab: +41 (0)44 633 23 01 >> E-mail: fcasola(a)phys.ethz.ch >> >> Group website: https://protos.ethz.ch/dscm >> > > > > End of Comp-phys-alps-users Digest, Vol 59, Issue 4 > *************************************************** > ---------------------------------------------------------------- SISSA Webmail https://webmail.sissa.it/ Powered by SquirrelMail http://www.squirrelmail.org/

2 1

Dzyaloshinskii-Moriya in spinmc
by Francesco Casola 06 Feb '11

06 Feb '11

Dear ALPS community, could you please let me know about the possibility of including Dzyaloshinskii-Moriya in spinmc, possibly through the presence of an off-diagonal exchange coupling constant? Is it necessary to modify the source code for such a case, or can I manage it in some other way not explicitly mentioned in the ALPS tutorial? Best wishes Francesco -- Francesco Casola, Ph. D. Student Laboratory for Solid State Physics ETH Zürich Hönggerberg HPF E15.1 Schafmattstrasse 16 CH-8093 Switzerland Tel. Office: +41 (0)44 633 22 96 / 03 Tel. Lab: +41 (0)44 633 23 01 E-mail: fcasola(a)phys.ethz.ch Group website: https://protos.ethz.ch/dscm

2 1

Extracting dates
by Xiao Quan Yu 06 Feb '11

06 Feb '11

Dear ALPS community, I met a trouble for extracting the calculated results from the XML output files. My setting is the following: plot2.xml file: <?xml version="1.0" encoding="UTF-8"?> <plot name="Stiffness versus disorder for three-dimesnional HardcoreBoson models"> <legend show="true"/> <xaxis label="Disorder" type="PARAMETER" name="mu"/> <yaxis label="Stiffness" type="SCALAR_AVERAGE"/> <set label="Inhomogeneous simple cubic lattice"/> </plot> and my "input" file is LATTICE_LIBRARY="lattices.xml"; LATTICE="inhomogeneous simple cubic lattice"; L=5; MODEL_LIBRARY="models.xml"; MODEL="hardcore boson"; t = 1.0; T = 0.1; delta=10.0; SWEEPS=500000; THERMALIZATION=10000; EPSILON=0.25; { DISORDERSEED = 34275; mu=delta*2*(random()-0.5); } { DISORDERSEED = 49802; mu=delta*2*(random()-0.5); } { DISORDERSEED = 82529; mu=delta*2*(random()-0.5); } { DISORDERSEED = 13564; mu=delta*2*(random()-0.5); } { DISORDERSEED = 65897; mu=delta*2*(random()-0.5); }. When I typed the command "extracttext plot2.xml input.task*.out.xml", the following thing happend: plot2.xml:4: parser error : XML declaration allowed only at the start of the document <?xml version="1.0" encoding="UTF-8"?> ^ unable to parse plot2.xml /tmp/archive2plot3652.xsl:1: parser error : Document is empty ^ /tmp/archive2plot3652.xsl:1: parser error : Start tag expected, '<' not found ^ cannot parse /tmp/archive2plot3652.xsl /tmp/plot3652:1: parser error : Document is empty ^ /tmp/plot3652:1: parser error : Start tag expected, '<' not found ^ unable to parse /tmp/plot3652 I don't know what is going on here, could you help me figure it out? Thank you in advance! Best, Xiaoquan ---------------------------------------------------------------- SISSA Webmail https://webmail.sissa.it/ Powered by SquirrelMail http://www.squirrelmail.org/

2 1

ising scheduler example
by Evgeny Burovskiy 31 Jan '11

31 Jan '11

Ladies and Gentlemen, I'm having a bit of a difficulty with something probably very very simple: I'm writing a simple MC code, which I just derive from the Ising scheduler example of ALPS 1.3.5. The problem appears once I try using the evaluate.C code from the same example. Even if I just take the code from an example, it throws a bad lexical cast exception. Looks like I can nail it down to a line, but the line itself is a bit of alea magic and I simply have no clue how to fix it. Hoping that somebody would spend a minute of their time and have a look at it, I'm attaching a minimal example which shows this behaviour. Here I just copy the sources from the scheduler example directory and only patch the source of evaluate.C. Below is the listing of what I'm doing. Any suggestions would be kindly appreciated, Zhenya ===================================================================== br@ymir:~/sweethome/adatoms/ising/from_scheduler$ make ising g++ -DNDEBUG -I/usr/local/boost -I/usr/local/ALPS-1.3.5/include -g -O2 -o ising ising.C main.C -L/usr/local/ALPS-1.3.5/lib -lalps -lcomm-sgl -llapack -lf77blas -latlas -lpthread -lm <.... a lot of complaints on deprecated headers...> br@ymir:~/sweethome/adatoms/ising/from_scheduler$ make evaluate g++ -DNDEBUG -I/usr/local/boost -I/usr/local/ALPS-1.3.5/include -g -O2 -o evaluate evaluate.C -L/usr/local/ALPS-1.3.5/lib -lalps -lcomm-sgl -llapack -lf77blas -latlas -lpthread -lm <.... a lot of complaints on deprecated headers...> br@ymir:~/sweethome/adatoms/ising/from_scheduler$ parameter2xml par1 Converting parameter file par1 to par1.in.xml br@ymir:~/sweethome/adatoms/ising/from_scheduler$ ./ising par1.in.xml Ising simulation example program using the ALPS Monte Carlo library copyright(c) 1994-2006 by Matthias Troyer <troyer(a)comp-phys.org> using the ALPS parallelizing scheduler copyright (c) 1994-2006 by Matthias Troyer <troyer(a)comp-phys.org>. see Lecture Notes in Computer Science, Vol. 1505, p. 191 (1998). based on the ALPS libraries version 1.3.5 available from http://alps.comp-phys.org/ copyright (c) 1994-2009 by the ALPS collaboration. Consult the web page for license details. For details see the publication: A.F. Albuquerque et al., J. of Magn. and Magn. Materials 310, 1187 (2007). parsing task files ... Created run 1 locally Starting task 1. Checking if it is finished: not yet, next check in 60 seconds ( 0% done). Halted Simulation 1 This task took 1 seconds. Checkpointing Simulation 1 Finished with everything. br@ymir:~/sweethome/adatoms/ising/from_scheduler$ ./evaluate par1.task1.out.xml parsing task files ... line 1 of evaluate line 2 of evaluate Caught exception: bad lexical cast: source type value could not be interpreted as target br@ymir:~/sweethome/adatoms/ising/from_scheduler$ ======================================================================

1 0

Fwd: FW: FW: error at installation for pyalea_c.so
by Krunoslav Prsa 23 Jan '11

23 Jan '11

Indeed, that was also the error in my case. Now it works, thanks! Kruno -------- Original Message -------- Subject: FW: [ALPS-users] FW: error at installation for pyalea_c.so Date: Sat, 22 Jan 2011 17:56:25 +0100 From: "Prsa Krunoslav" <Krunoslav.Prsa(a)psi.ch> To: <kprsa(a)phys.ethz.ch> ------------------------------------------- From: comp-phys-alps-users-bounces(a)lists.phys.ethz.ch on behalf of Munehisa Matsumoto[SMTP:MATSUMOTO@PHYSICS.UCDAVIS.EDU] Sent: Saturday, January 22, 2011 5:56:18 PM To: comp-phys-alps-users(a)lists.phys.ethz.ch Subject: Re: [ALPS-users] FW: error at installation for pyalea_c.so Auto forwarded by a Rule Dear Benoît, Thank you for sending the file "cmake_install.cmake" which was taken from lib/pyalps under your build directory as I understood. It looks like your build directory is the same as your installation directory. Is that the case? In line 41 of the file your install destination reads ${CMAKE_INSTALL_PREFIX}/lib/pyalps with ${CMAKE_INSTALL_PREFIX} being "/home/gremaud/alps/alps-2.0.0-r5363" as defined in line 5, while from the build directory "pyalea_c.so" is referenced as "/home/gremaud/alps/alps-2.0.0-r5363/lib/pyalps/pyalea_c.so" I guess your installation moves the file "pyalea.so" from your build directory to your installation directory and thus the file "pyalea.so" is deleted after this stage of installation is supposed to be passed. If the build directory and the installation directory is the same, it might happen that a target file is copied onto itself and then deleted. To avoid such complication having a separate installation directory might help. At least now I am reproducing the problem by setting the build directory and the installation directory to be the same. best regards, Munehisa -- ------------------------------------------------- Dr. Krunoslav Prsa, Postdoc Dynamics of strongly correlated materials Lab. for Solid State Physics Building/Office: HPF/E15.1 Schafmattstrasse 16, ETH Zürich Hönggerberg, CH-8093 Switzerland. Office phone: +41 44 633 22 96 Lab phone: +41 44 633 23 01 Mobile phone: +41 76 386 17 99 -------------------------------------------------

1 0