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

transverse Ising model with DMRG
by shuai cui 26 May '11

26 May '11

Dear Prof. Troyer I am trying to calculate the ground state energy of the transverse Ising-Model ( -J*4*Sx(i)*Sx(j)-h*2*Sz(i) ) on a period 1d-chain-lattice with the DMRG. Though I had use both the default model and a modified model, the results are not satisfactory. So I ask for why the DMRG doesn't give a good result in transverse field Ising model as in other 1D model. the modified model is: <MODELS> <SITEBASIS name="spin"> <PARAMETER name="local_spin" default="local_S"/> <PARAMETER name="local_S" default="1/2"/> <QUANTUMNUMBER name="S" min="local_spin" max="local_spin"/> <QUANTUMNUMBER name="Sz" min="-S" max="S"/> <OPERATOR name="Splus" matrixelement="sqrt(S*(S+1)-Sz*(Sz+1))"> <CHANGE quantumnumber="Sz" change="1"/> </OPERATOR> <OPERATOR name="Sminus" matrixelement="sqrt(S*(S+1)-Sz*(Sz-1))"> <CHANGE quantumnumber="Sz" change="-1"/> </OPERATOR> <OPERATOR name="Sz" matrixelement="Sz"/> </SITEBASIS> <BASIS name="spin"> <SITEBASIS ref="spin"> <PARAMETER name="local_spin" value="local_S#"/> <PARAMETER name="local_S#" value="local_S"/> <PARAMETER name="local_S" value="1/2"/> </SITEBASIS> <CONSTRAINT quantumnumber="Sz" value="Sz_total"/> </BASIS> <SITEOPERATOR name="Sx" site="x"> 1/2*(Splus(x)+Sminus(x)) </SITEOPERATOR> <HAMILTONIAN name="trans_Ising"> <PARAMETER name="J" default="0"/> <PARAMETER name="h" default="0"/> <BASIS ref="spin"/> <SITETERM site="i"> h*2*Sz(i) </SITETERM> <BONDTERM source="i" target="j"> J*4*Sx(i)*Sx(j) </BONDTERM> </HAMILTONIAN> </MODELS> Separately, the parameter files are: LATTICE="chain lattice" MODEL= "spin" Jz=-4 Gamma=2 local_S=1/2 SWEEPS=4 NUMBER_EIGENVALUES=1 L=8 {MAXSTATES=20} LATTICE="chain lattice" MODEL= "trans_Ising" J=-1 h=-1 local_S=1/2 SWEEPS=4 NUMBER_EIGENVALUES=1 L=8 {MAXSTATES=20} Correspondingly, the results are: FINAL SWEEP - LAST ITERATION NUMBER OF STATES: 8 2 2 8 256 E0 = -282.636028178 ITER = 3 ENTROPY = 1.78546641552 ------------------------------------------- Truncation error = -4.98732999343e-17 ------------------------------------------- Checking if it is finished: not yet, next check in 60 seconds ( 0% done). Halted Simulation 1 This task took 46 seconds. Checkpointing Simulation 1 Finished with everything. FINAL SWEEP - LAST ITERATION NUMBER OF STATES: 8 2 2 8 256 E0 = -9.11126807853 ITER = 3 ENTROPY = 0.931162271124 ------------------------------------------- Truncation error = 2.22044604925e-16 ------------------------------------------- Checking if it is finished: not yet, next check in 60 seconds ( 0% done). Halted Simulation 1 This task took 13 seconds. Checkpointing Simulation 1 Finished with everything. but the correct data should be -10.251661790966022 which calculated by sparsediag. I don't understand why DMRG give a low precision result, even I had set h=0. Remarkably, the default model gives a wrong data, which far away from the correct one. I had exercise the XX model, which gives a very well data. I wonder if the transverse field Ising model has some special feature, or the DMRG parameter I used is not right? best regards, Shuai Cui

1 0

tebd will not be built
by Cristian Degli Esposti Boschi 25 May '11

25 May '11

Dear users, I'm installing ALPS 2 on a Debian 6 (squeeze) machine. After issuing cmake -DCMAKE_INSTALL_PREFIX... I get the following (I extracted the last relevant part): ... Falling back to CMake provied LAPACK/BLAS detection. -- A library with BLAS API found. -- A library with BLAS API found. -- A library with LAPACK API found. -- SQLite Library: not found -- LPSolve Library: not found -- Could NOT find SZIP (missing: SZIP_LIBRARIES SZIP_INCLUDE_DIRS) -- HDF5 version found is 1.8.4-patch1 in /usr/include -- Python interpreter /usr/bin/python -- Python interpreter ok : version 2.6.6 -- PYTHON_INCLUDE_DIRS = /usr/include/python2.6 -- PYTHON_NUMPY_INCLUDE_DIR = /usr/lib/pymodules/python2.6/numpy/core/include -- PYTHON_SITE_PKG = /usr/lib/python2.6/dist-packages -- PYTHON_LIBRARY = /usr/lib/python2.6/config/libpython2.6.so -- PYTHON_EXTRA_LIBS =-lssl -lcrypto -lssl -lcrypto -L/usr/lib -lz -lpthread -ldl -lutil -- PYTHON_LINK_FOR_SHARED = -Xlinker -export-dynamic -Wl,-O1 -Wl,-Bsymbolic-functions -- Numpy include in /usr/lib/pymodules/python2.6/numpy/core/include -- ALPS XML dir is /usr/local/ALPS/alps-2.0.1-r5565-src-with-boost/alps/lib/xml -- Adding Boost dir: /usr/local/ALPS/alps-2.0.1-r5565-src-with-boost/boost -- HDF5 version found is 1.8.4-patch1 in /usr/include -- HDF5 version found is 1.8.4-patch1 in /usr/include -- tebd will not be built -- Configuring done -- Generating done -- Build files have been written to: /usr/local/ALPS/alps-2.0.1-r5565-src-with-boost Now, my problem is that "tebd will not be built" and I wonder why. Is it related to the problem of "Falling back to CMake provied LAPACK/BLAS detection" or maybe to what reported in the installation guidelines that "The tebd code requires a fortran compiler and that HDF5 be built with the --enable-fortran configure option." In the latter case how can I fix the things in Debian-provided HDF5? Thanks for your time, Cristian Degli Esposti Boschi CNR-IMM, Sezione di Bologna, via Gobetti, 101, 40129, Bologna, Italia tel. ++39 051 6399152, fax ++39 051 6399216 degliesposti -AT- bo.imm.cnr.it

2 1

parameters in spinmc
by rosales＠fisica.unlp.edu.ar 23 May '11

23 May '11

Dear all, I have a few questions about parameters in the spinmc application. In the Tutorias: Parameters web page, I found THERMALIZATION: Number of Monte Carlo sweeps for thermalization SWEEPS: number of Monte Carlo steps (after thermalization) -I don't understand what is the number of "copys" used to compute the averages. -Between copy and copy to compute the averages, how many MC steps are there? - Suppose I run with different magnetic fields. Is it possible to use the previous state as initial state for the next run? Best regards, > Dear all, > > I have a few questions about the observables that can be measured by the > spinmc application. > - It is possible to compute the static magnetic structure factor? > S(q)=Sum_{ij} <S_i.S_j> Exp(i*q(ri-rj)) > or > > - It is possible to compute the average of correlations <S_i.S_j> ? if > not, there is possibility to include now such terms? > > Best regards, > > Hector > > > -- > Dr. Hector Diego Rosales > IFLP - Dto. de Física UNLP, > Calle 49 y 115 s/n, > 1900 La Plata, > Argentina > Tel: +54 221 424 6062; > Fax: +54 221 423 0122 > e-Mail: rosales(a)fisica.unlp.edu.ar > -- Dr. Hector Diego Rosales IFLP - Dto. de Física UNLP, Calle 49 y 115 s/n, 1900 La Plata, Argentina Tel: +54 221 424 6062; Fax: +54 221 423 0122 e-Mail: rosales(a)fisica.unlp.edu.ar

3 5

tools (extracttext) missing?
by Martin Hohenadler 22 May '11

22 May '11

Dear all, I just installed the ALPS 2.0.1 Ubuntu package alps2all on Natty. Things work out of the box, but I can not find the xml tools such as extracttext anywhere. I don't remember doing anything extra the last time I installed ALPS. Thanks, Martin

3 4

tutorials
by Tadeusz Wasiutyński 12 May '11

12 May '11

Hello, In the tutorial on the alps web page line: import pyalps.pyplot should be rather: import pyalps.plot as it is in the files stored in distribution. Regards Tadeusz Wasiutynski

2 1

Re: [ALPS-users] another installation error
by Kubila 11 May '11

11 May '11

Dear Prof. Troyer Thank you very much for your email. In the attachment please find a copy of the file you request. I appreciate your help a lot. Sincerely, kubila --- On Tue, 5/10/11, Matthias Troyer <troyer(a)phys.ethz.ch> wrote: From: Matthias Troyer <troyer(a)phys.ethz.ch> Subject: Re: [ALPS-users] another installation error To: comp-phys-alps-users(a)lists.phys.ethz.ch Date: Tuesday, May 10, 2011, 5:53 PM It seems that the pthreads library is not available on your machine or not linked. Can you send a copy of your CMakeCache.txt? Matthias On May 10, 2011, at 2:29 AM, Kubila wrote: Dear Prof. Troyer Using the official release gcc 4.1.2, this issue is gone. On the other hand, when installing the latest alps 2.0.1, I got another installation error as follows [ 52%] Building CXX object tool/CMakeFiles/compactrun.dir/compactrun.C.o Linking CXX executable compactrun ../src/boost/libboost.so: undefined reference to `pthread_create' ../src/boost/libboost.so: undefined reference to `pthread_key_create' ../src/boost/libboost.so: undefined reference to `pthread_getspecific' ../src/boost/libboost.so: undefined reference to `pthread_once' ../src/boost/libboost.so: undefined reference to `pthread_join' ../src/boost/libboost.so: undefined reference to `pthread_setspecific' ../src/boost/libboost.so: undefined reference to `pthread_detach' collect2: ld returned 1 exit status make[2]: *** [tool/compactrun] Error 1 make[1]: *** [tool/CMakeFiles/compactrun.dir/all] Error 2 make: *** [all] Error 2 Any help will be greatly appreciated. Sincerely, kubila --- On Mon, 5/2/11, Matthias Troyer <troyer(a)phys.ethz.ch> wrote: From: Matthias Troyer <troyer(a)phys.ethz.ch> Subject: Re: [ALPS-users] installation error To: comp-phys-alps-users(a)lists.phys.ethz.ch Date: Monday, May 2, 2011, 5:58 AM Since the linker crashes it would be either a bug in the compiler or the linker. Can you try a more recent compiler than a prerelease of gcc 4.1? On May 2, 2011, at 12:12 AM, Kubila wrote: Dear Prof. Troyer Thank you very much for your email. The operating system is Debian GNU/Linux 4.0 (etch). The compiler is gcc (gcc version 4.1.2 20061115 (prerelease) (Debian 4.1.1-21)). Will the operating system be too " old " so that it is impossible to install ALPS 2.0.1? Thanks again for your reply. Sincerely, kubila From: Matthias Troyer <troyer(a)phys.ethz.ch> To: comp-phys-alps-users(a)lists.phys.ethz.ch Sent: Mon, May 2, 2011 12:11:17 AM Subject: Re: [ALPS-users] installation error Hi, this is the linker of your system crashing. Which operating systems is this and which compiler? Matthias On May 1, 2011, at 8:43 AM, Kubila wrote: Hi all I tried to install the nightly built ALPS 2.0.1 in a PC cluster and ended up with the following error message Linking CXX shared library libalps.so collect2: ld terminated with signal 11 [Segmentation fault] /usr/bin/ld: make[2]: *** [src/alps/libalps.so.2.0.1] Error 1 make[1]: *** [src/alps/CMakeFiles/alps.dir/all] Error 2 make: *** [all] Error 2 I am not the system administration so I tried to install ALPS into my home directory. The installation process moved on to about 51% and I encountered the error quoted above..... any help will appreciate much. Sincerely, kubila

1 0

installation error
by Kubila 11 May '11

11 May '11

Hi all I tried to install the nightly built ALPS 2.0.1 in a PC cluster and ended up with the following error message Linking CXX shared library libalps.so collect2: ld terminated with signal 11 [Segmentation fault] /usr/bin/ld: make[2]: *** [src/alps/libalps.so.2.0.1] Error 1 make[1]: *** [src/alps/CMakeFiles/alps.dir/all] Error 2 make: *** [all] Error 2 I am not the system administration so I tried to install ALPS into my home directory. The installation process moved on to about 51% and I encountered the error quoted above..... any help will appreciate much. Sincerely, kubila

2 5

multi-site terms support in hamiltonian
by Oles Zaburannyi 02 May '11

02 May '11

Hello, I have question about multi-site terms support in hamiltonian definition. I found some mention about that in ALPS workshop (2006) In recent wiki announced that `Extensions to more complex terms, such as 3- and 4-site terms are in preparation and will be included here as soon as the ALPS libraries support such terms.` There is possibility to include now such terms? Will arise huge problems in DMRG methods after including such milti-site terms? Thanks, Oles

2 3

Simulation of O(N) models
by Rachele Nerattini 28 Apr '11

28 Apr '11

Hi ALPS community, first of all sorry for my English and sorry if the questions that I'm going to make are trivial but I'm not a real expert in programming and using web libraries. I'm interested in doing classical MC simulations of O(N) models. I saw on the ALPS page Documentation:ClassicalMCSimulations and in particular ALPS application: spinmc that is possible to extend the application to additional O(N) models "...by editing the file mc/spins/spinmc_factory.C " I need to run a simulation for a O(4) model so I tried to edit the file spinmc_factory.C to this purpose in the following way /***************************************************************************** * * ALPS Project Applications * * Copyright (C) 2009 by Matthias Troyer <troyer(a)comp-phys.org>, * Synge Todo <wistaria(a)comp-phys.org> * * This software is part of the ALPS Applications, published under the ALPS * Application License; you can use, redistribute it and/or modify it under * the terms of the license, either version 1 or (at your option) any later * version. * * You should have received a copy of the ALPS Application License along with * the ALPS Applications; see the file LICENSE.txt. If not, the license is also * available from http://alps.comp-phys.org/. * * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, * FITNESS FOR A PARTICULAR PURPOSE, TITLE AND NON-INFRINGEMENT. IN NO EVENT * SHALL THE COPYRIGHT HOLDERS OR ANYONE DISTRIBUTING THE SOFTWARE BE LIABLE * FOR ANY DAMAGES OR OTHER LIABILITY, WHETHER IN CONTRACT, TORT OR OTHERWISE, * ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER * DEALINGS IN THE SOFTWARE. * *****************************************************************************/ /* $Id: spinmc_factory.C 3124 2009-02-26 09:40:13Z wistaria $ */ #include "factory.h" #include "abstractspinsim.h" #include "spinsim.h" #include "ising.h" #include "on.h" #include "potts.h" void SpinFactory::print_copyright(std::ostream& out) const { out << "Generic classical Monte Carlo program using local or cluster updates\n" << " available from http://alps.comp-phys.org/\n" << " copyright(c) 1999-2007 by Matthias Troyer <troyer(a)comp-phys.org>\n" << " Mathias Koerner < mkoerner(a)comp-phys.org>\n" << " for details see the publication:\n" << " A.F. Albuquerque et al., J. of Magn. and Magn. Materials 310, 1187 (2007).\n\n"; } alps::scheduler::MCSimulation* SpinFactory::make_task( const alps::ProcessList& w, const boost::filesystem::path& fn) const { return new alps::scheduler::MCSimulation(w,fn); } /** * Returns the number of white-space separated elements in the string. * * \param str the string to be analyzed. */ int SpinFactory::countElements(const std::string& str) const { int counter = 0; int pos = 0; if (str[0] == '\0') return 0; while (true) { while (str[pos] == ' ') pos++; if (str[pos] == '\0') return counter; counter++; while ((str[pos] !='\0') && (str[pos] != ' ')) pos++; } } /** * Returns the number of elements in the string with the largest number of elements. */ int SpinFactory::findDominantMatrixString(const alps::Parameters& parms) const { int count = 0; /* should be done more elegantly ... */ if (parms.defined("J")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["J"])); if (parms.defined("J'")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["J'"])); if (parms.defined("J0")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["J0"])); if (parms.defined("J1")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["J1"])); if (parms.defined("J2")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["J2"])); if (parms.defined("J3")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["J3"])); if (parms.defined("D")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["D"])); if (parms.defined("D'")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["D'"])); if (parms.defined("D0")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["D0"])); if (parms.defined("D1")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["D1"])); if (parms.defined("D2")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["D2"])); if (parms.defined("D3")) count = std::max BOOST_PREVENT_MACRO_SUBSTITUTION (count, countElements(parms["D3"])); return count; } /** * Outputs an error message and aborts the computation if an invalid number of * elements is provided in the input for a matrix parameter. */ void SpinFactory::produceError(const alps::Parameters& parms) const { std::cerr << "Invalid parameters given:\n" << "Cannot produce model " << parms["MODEL"] << " with strings " << parms["D"] << " for D and " << parms["J"] << " for J.\nAborting\n"; boost::throw_exception(std::runtime_error("invalid parameters")); std::exit(-5); } /** * Creates a worker, according to the parameters. The type of the produced * worker depends on the model given in the parameter file and on the number of * input elements for the matrix parameters. The chosen matrix type is the most * general one (over all matrix parameters) of all most specialized ones that * are compatible with the user input. */ alps::scheduler::Worker* SpinFactory::make_worker( const alps::ProcessList& where, const alps::Parameters& parms,int node) const { int maxElemCount = findDominantMatrixString(parms); if (parms["MODEL"]=="Ising") if (maxElemCount == 1) return new SpinSim<IsingMoment,MIdMatrix<double,1> >(where,parms,node); else { produceError(parms); return 0; } else if (parms["MODEL"]=="O(4)") if (maxElemCount == 1) return new SpinSim< ONMoment<4>, MIdMatrix<double,4> >(where, parms, node); else { produceError(parms); return 0; } else if (parms["MODEL"]=="Heisenberg") switch (maxElemCount) { case 1: return new SpinSim< ONMoment<3>, MIdMatrix<double,3> >(where, parms, node); break; case 3: return new SpinSim< ONMoment<3>, DiagMatrix<double,3> >(where, parms, node); break; case 6: return new SpinSim< ONMoment<3>, SquareMatrix<double,3> >(where, parms, node); break; case 9: return new SpinSim< ONMoment<3>, SquareMatrix<double,3> >(where, parms, node); break; default: produceError(parms); return 0; } else if (parms["MODEL"]=="XY") switch (maxElemCount) { case 1: return new SpinSim< XYMoment, MIdMatrix<double,2> >(where, parms, node); break; case 2: return new SpinSim< XYMoment, DiagMatrix<double,2> >(where, parms, node); break; case 3: return new SpinSim< XYMoment, SquareMatrix<double,2> >(where, parms, node); break; case 4: return new SpinSim< XYMoment, SquareMatrix<double,2> >(where,parms,node); break; default: produceError(parms); return 0; } else if (parms["MODEL"]=="Potts") switch (int(parms["q"])) { case 3: return new SpinSim< PottsMoment<3>, MIdMatrix<double,2> >(where,parms,node); case 4: return new SpinSim< PottsMoment<4>, MIdMatrix<double,2> >(where,parms,node); case 10: return new SpinSim< PottsMoment<10>, MIdMatrix<double,2> >(where,parms,node); default: boost::throw_exception(std::runtime_error(std::string(parms["q"]) + "-state Potts model not implemented in file factory.C")); } else boost::throw_exception(std::runtime_error(std::string(parms["MODEL"]) + " model not implemented in file factory.C")); return 0; } unfortunately I don't know this programming language so I'm not sure to have edited it in the right way. Is it correct? And if it is so, what do I have to do after having edited the file to run the simulation? Thank you for the help. Best regards, Rachele Nerattini

2 1

evaluating simulation results in python in ALPS 2
by Evgeny Burovskiy 26 Apr '11

26 Apr '11

Dear all, In an attempt to get to grips with ALPS 2, I've ran the example from example/parapack/exchange directory. When compiled and run from the command line, it produces both xml and hdf output files: br@ymir:~/sweethome/adatoms/par_temp/parapack_example$ ll total 10064 drwxr-xr-x 2 br br 4096 2011-04-25 18:57 ./ drwxr-xr-x 3 br br 4096 2011-04-25 16:06 ../ -rw-r--r-- 1 br br 13276 2011-04-25 16:17 ALPS.xsl -rw-r--r-- 1 br br 1862 2011-04-25 16:11 ising.C -rw-r--r-- 1 br br 6696 2011-04-25 16:06 ising.h -rw-r--r-- 1 br br 1353 2011-04-25 16:05 main.C -rw-r--r-- 1 br br 483 2011-04-25 16:09 Makefile -rw-r--r-- 1 br br 10068 2011-04-25 16:20 mpir.txt -rwxr-xr-x 1 br br 752957 2011-04-25 16:11 pa_ex* -rw-r--r-- 1 br br 217 2011-04-25 16:06 params-ising -rw-r--r-- 1 br br 463 2011-04-25 16:12 params-ising.in.xml -rw-r--r-- 1 br br 729 2011-04-25 16:13 params-ising.out.xml -rw-r--r-- 1 br br 46708 2011-04-25 16:13 params-ising.task1.clone1 -rw-r--r-- 1 br br 2831368 2011-04-25 16:13 params-ising.task1.clone1.h5 -rw-r--r-- 1 br br 780 2011-04-25 16:12 params-ising.task1.in.xml -rw-r--r-- 1 br br 168452 2011-04-25 18:57 params-ising.task1.out.xml -rw-r--r-- 1 br br 398984 2011-04-25 18:57 params-ising.task1.replica10.h5 -rw-r--r-- 1 br br 398984 2011-04-25 18:57 params-ising.task1.replica11.h5 -rw-r--r-- 1 br br 398984 2011-04-25 18:57 params-ising.task1.replica12.h5 -rw-r--r-- 1 br br 398984 2011-04-25 18:57 params-ising.task1.replica13.h5 -rw-r--r-- 1 br br 398984 2011-04-25 18:57 params-ising.task1.replica14.h5 -rw-r--r-- 1 br br 398984 2011-04-25 18:57 params-ising.task1.replica15.h5 -rw-r--r-- 1 br br 388560 2011-04-25 18:57 params-ising.task1.replica16.h5 -rw-r--r-- 1 br br 407752 2011-04-25 18:57 params-ising.task1.replica1.h5 -rw-r--r-- 1 br br 398984 2011-04-25 18:57 params-ising.task1.replica2.h5 -rw-r--r-- 1 br br 398984 2011-04-25 18:57 params-ising.task1.replica3.h5 -rw-r--r-- 1 br br 398984 2011-04-25 18:57 params-ising.task1.replica4.h5 -rw-r--r-- 1 br br 398984 2011-04-25 18:57 params-ising.task1.replica5.h5 -rw-r--r-- 1 br br 398984 2011-04-25 18:57 params-ising.task1.replica6.h5 -rw-r--r-- 1 br br 398984 2011-04-25 18:57 params-ising.task1.replica7.h5 -rw-r--r-- 1 br br 398984 2011-04-25 18:57 params-ising.task1.replica8.h5 -rw-r--r-- 1 br br 398984 2011-04-25 18:57 params-ising.task1.replica9.h5 -rwxr--r-- 1 br br 227 2011-04-25 18:54 readin.py* Then I'm trying to load the simulation results using python tools, and that's where I seem to hit a snag. I follow the tutorial at <http://alps.comp-phys.org/mediawiki/index.php/ALPS_using_Python#Loading_the…>, here's the exact script: br@ymir:~/sweethome/adatoms/par_temp/parapack_example$ cat readin.py #!/usr/bin/python import pyalps import matplotlib.pyplot as plt import pyalps.plot result_files = pyalps.getResultFiles(prefix='params-ising') print result_files obslist=pyalps.loadObservableList(result_files) This fails with the following error message: br@ymir:~/sweethome/adatoms/par_temp/parapack_example$ ./readin.py ['./params-ising.task1.out.xml'] Traceback (most recent call last): File "./readin.py", line 8, in <module> obslist=pyalps.loadObservableList(result_files) File "/home/br/ALPS2/lib/pyalps/load.py", line 611, in loadObservableList res = ll.GetProperties(files,proppath,respath,verbose=verbose) File "/home/br/ALPS2/lib/pyalps/load.py", line 126, in GetProperties fs = self.GetFileNames(flist) File "/home/br/ALPS2/lib/pyalps/load.py", line 108, in GetFileNames log( "FILE ", f, "DOES NOT EXIST!") TypeError: log() takes exactly 1 argument (3 given) Running the same sequence of commands in alpspython interpreter produces the same error message. From what I can see, there are two issues here: 1) GetFileNames() function of the Hdf5Loader class does not handle the replica#.h5 file endings. 2) there's a problem with log() function which the interpreter reports, Or am I doing something wrong here? Or, more to the point, how do I load the simulation results into python? Thanks in advance, Zhenya

2 1