March 2014 - Comp-phys-alps-users

Inconsistent comment in generate_spline_matrix in the DMFT folder
by Kuangshing Chen 22 Jul '15

22 Jul '15

Hi, The function generate_spline_matrix in fouriertransform.C in the dmft folder shows the comments as follows: // A is the matrix whose inverse defines spline_matrix // // 6 6 // 1 4 1 // 1 4 1 // A = ... // // 1 4 1 // -2 0 0 2 However, the following codes, dense_matrix A = 4*dt/6.*boost::numeric::ublas::identity_matrix<double>(Np1); for (int i=1; i<Np1-1; i++) { A(i,i-1) = dt/6.; A(i,i+1) = dt/6.; } A(0,0) = 1.; A(0, Np1-1) = 1.; A(Np1-1, 0) = -2.*dt/6.; A(Np1-1, 1) = -1.*dt/6.; A(Np1-1, Np1-2) = 1*dt/6.; A(Np1-1, Np1-1) = 2*dt/6.; the red lines will make the matrix looks like // 6 6 // 1 4 1 // 1 4 1 // A = ... // // 1 4 1 // -2 *-1* *1* 2 Is the comment correct or the code correct? Any reference for this algorithm? Thank you, Kuang-Shing Chen

5 25

ALPS worm simulations hung after re-starting from checkpoints
by vvarma＠ictp.it 15 Apr '14

15 Apr '14

Dear ALPS users, I have a large list of tasks run from a given job file *parm.in.xml*, which got terminated after some tasks were completed; my *parm.out.xml* looks like this after the termination: ===========<SNIP-1>================= /......// //......// // <TASK status="finished">// // <INPUT file="parm.task62.out.xml"/>// // </TASK>// // <TASK status="finished">// // <INPUT file="parm.task63.out.xml"/>// // </TASK>// // <TASK status="new">// // <INPUT file="parm.task64.out.xml"/>// // </TASK>// // <TASK status="new">// //......// //....../ ============<SNIP-1>================ To restart the job from task64, I execute the command (without changing anything in the output files) *worm**--Tmin 10 parm.out.xml* Could you please confirm that this is correct? I ask because nothing seems to proceed after a while after executing the above: ============<SNIP-2>=============== /.........// //.........// //Loading information about run 1 from file parm.task62.out.run1// //Loading information about run 1 from file parm.task63.out.run1// //Task 1 finished.// //........// //........// / /Task 63 finished.// //Created run 1 locally// //Starting task 64.// //Checking if it is finished: not yet, next check in 10 seconds ( 0% done).// //Checking if it is finished: not yet, next check in 203 seconds ( 1% done).// //Checking if it is finished: not yet, next check in 174 seconds ( 23% done)./ =============<SNIP-2>============== >From thereon the simulation just hangs for many hours (I've tried this procedure repeatedly), which it should not because my SWEEPS and THERMALIZATION are not excessively large. Any ideas on why the re-started simulations don't seem to proceed are appreciated. With regards, Vipin

3 4

Complex number
by Tianyou YI 28 Mar '14

28 Mar '14

Hello, My model need to used complex number. Then I wrote the value of a parameter as "s=0.5*sqrt(-1)" and I got the error message "Cannot evaluate I". What is wrong with me? What can I do to input a complex number? Thank you for your reply. Best wishes Tianyou

1 0

Problems downloading libraries
by Niraj Ghimire 27 Mar '14

27 Mar '14

Dear alps users, I have been trying to build alps package from source on Windows 7 OS. I have been following the outdated instructions on the page https://alps.comp-phys.org/trac/wiki/AlpsOnWindows. I could not find the proper versions of Zlib, SZlib, Lpsolve and SQlite. There are links on this webpage that read "Check out the source code from <>" where <> is, for instance, https://alps.comp-phys.org/svn/alps1/trunk/wistaria/windows/lpsolve for "Lpsolve". I don't know how to download the source codes using these links. I tried to download the source codes from a trusted webpage like sourceforge.net, some of them work while others don't and I keep getting errors. Please help me. Regards, Niraj Ghimire Graduate Student Department of Physics University of Connecticut

3 24

complexe number
by Tianyou YI 27 Mar '14

27 Mar '14

Hello, In my model, some parameters need to be a complexe number. How can I give a complexe number in the model file. Thank you for your reply. Best wishes Tianyou

1 0

Nearest-neighbor interaction
by Tianyou YI 25 Mar '14

25 Mar '14

Hello, I need to create my own Hamiltonian. I have a term of nearest-neighbor interaction, and I write in the model file </BONDTERM> <BONDTERM source="i" target="j"> -t*bdag(i)*b(j) </BONDTERM> Here in the field of the target should I write target="i+1" or should I write simply target="j". What is the difference between the two cases ? In the second case, the source i will go through all the sites and the target j will go through alll the sites too, I have then interaction of next-nearest-neighbor terms and so on ? Thank you for your reply. Best wishes Tianyou

2 1

ALPS 2.2b1 'make' error on Linux Mint 13 Maya
by vvarma＠ictp.it 25 Mar '14

25 Mar '14

Dear ALPS users/developers, After cmake and ccmake, I am trying to 'make' the ALPS 2.2b1 files on my computer (running on Linux Mint 13 Maya, which is in turn based on Ubuntu 12.04); however I get the following error message: /*<SNIP>*/ Linking CXX static library libspinmc_impl.a [ 67%] Built target spinmc_impl Scanning dependencies of target spinmc [ 68%] Building CXX object applications/mc/spins/CMakeFiles/spinmc.dir/spinmc.C.o In file included from /afs/ictp.it/home/v/vvarma/Desktop/alps-2.2.b1-r7195-src-with-boost/boost/boost/python/detail/wrap_python.hpp:50:0, from /afs/ictp.it/home/v/vvarma/Desktop/alps-2.2.b1-r7195-src-with-boost/boost/boost/python/detail/prefix.hpp:13, from /afs/ictp.it/home/v/vvarma/Desktop/alps-2.2.b1-r7195-src-with-boost/boost/boost/python/args.hpp:8, from /afs/ictp.it/home/v/vvarma/Desktop/alps-2.2.b1-r7195-src-with-boost/boost/boost/python.hpp:11, from /afs/ictp.it/home/v/vvarma/Desktop/alps-2.2.b1-r7195-src-with-boost/alps/src/alps/alea/mcdata.hpp:79, from /afs/ictp.it/home/v/vvarma/Desktop/alps-2.2.b1-r7195-src-with-boost/alps/src/alps/alea.h:34, from /afs/ictp.it/home/v/vvarma/Desktop/alps-2.2.b1-r7195-src-with-boost/alps/src/alps/alea/observableset.h:49, from /afs/ictp.it/home/v/vvarma/Desktop/alps-2.2.b1-r7195-src-with-boost/alps/src/alps/scheduler/montecarlo.h:40, from /afs/ictp.it/home/v/vvarma/Desktop/alps-2.2.b1-r7195-src-with-boost/alps/src/alps/scheduler.h:31, from /afs/ictp.it/home/v/vvarma/Desktop/alps-2.2.b1-r7195-src-with-boost/alps/applications/mc/spins/spinmc.C:29: /usr/include/python2.7/pyconfig.h:1161:0: warning: "_POSIX_C_SOURCE" redefined [enabled by default] /usr/include/features.h:164:0: note: this is the location of the previous definition /usr/include/python2.7/pyconfig.h:1183:0: warning: "_XOPEN_SOURCE" redefined [enabled by default] /usr/include/features.h:166:0: note: this is the location of the previous definition make[2]: *** No rule to make target `/usr/lib/libhdf5.so', needed by `applications/mc/spins/spinmc'. Stop. make[1]: *** [applications/mc/spins/CMakeFiles/spinmc.dir/all] Error 2 make: *** [all] Error 2 /*<SNIP>*/ I have tried reinstalling hdf5 (libhdf5-serial-dev) and the ALPS multiple times on this computer but there is always some error related to the hdf5 files, either during the make process or the cmake configuration. If you could please help me debug this particular error, that would be greatly appreciated. A few technical notes are appended at the end. Thanks so much. Regards, Vipin ====================================== 1) Surprisingly I am able to build and install ALPS without any problems on a laptop running native Ubuntu 12.04. 2) *cmake command* cmake -D Boost_ROOT_DIR:PATH=~/Desktop/alps-2.2.b1-r7195-src-with-boost/boost -DALPS_INSTALL_VISTRAILS_PACKAGES=ON -D VISTRAILS_APP_DIR:PATH=~/Desktop/vistrails-src-nightly-09cb1cfff309/ -DCMAKE_INSTALL_PREFIX=~/Desktop/alps-2.2.b1-r7195-src-with-boost/opt_alps ~/Desktop/alps-2.2.b1-r7195-src-with-boost/alps 3) *Snippet of CMakeCache.txt* //HDF5 C++ Wrapper compiler. Used only to detect HDF5 compile // flags. HDF5_CXX_COMPILER_EXECUTABLE:FILEPATH=HDF5_CXX_COMPILER_EXECUTABLE-NOTFOUND //HDF5 Wrapper compiler. Used only to detect HDF5 compile flags. HDF5_C_COMPILER_EXECUTABLE:FILEPATH=HDF5_C_COMPILER_EXECUTABLE-NOTFOUND //Path to a file. HDF5_C_INCLUDE_DIR:PATH=/usr/include //HDF5 file differencing tool. HDF5_DIFF_EXECUTABLE:FILEPATH=HDF5_DIFF_EXECUTABLE-NOTFOUND //Path to a file. HDF5_DLL:FILEPATH=HDF5_DLL-NOTFOUND //Path to a file. HDF5_HL_DLL:FILEPATH=HDF5_HL_DLL-NOTFOUND //HDF5 library compiled with parallel IO support HDF5_IS_PARALLEL:BOOL=FALSE //The HDF5_hdf5 library HDF5_hdf5_LIBRARY:FILEPATH=/usr/include //Path to a library. HDF5_hdf5_LIBRARY_DEBUG:FILEPATH=HDF5_hdf5_LIBRARY_DEBUG-NOTFOUND //Path to a library. HDF5_hdf5_LIBRARY_RELEASE:FILEPATH=/usr/lib/libhdf5.so //The HDF5_hdf5_hl library HDF5_hdf5_hl_LIBRARY:FILEPATH=/usr/include //Path to a library. HDF5_hdf5_hl_LIBRARY_DEBUG:FILEPATH=HDF5_hdf5_hl_LIBRARY_DEBUG-NOTFOUND //Path to a library. HDF5_hdf5_hl_LIBRARY_RELEASE:FILEPATH=/usr/lib/libhdf5_hl.so

4 5

How to build package from source
by Niraj Ghimire 22 Mar '14

22 Mar '14

Dear alps users, I want to use alps package with my 32 bit version of Python 2.7 on Windows 7 OS. I have been told that I need to build the package from source but I have no idea how to do this. Could you please provide me with a clear step-by-step procedure ? I would really appreciate your help. Regards, Niraj Ghimire Graduate Student Department of Physics University of Connecticut

2 2

Using DISORDERSEED in Bose glass worm simulations: Python script
by vvarma＠ictp.it 14 Mar '14

14 Mar '14

Dear ALPS users, I am trying to run a very simple worm simulation with disorder (site-dependent 'mu') for the Bose-Hubbard model using a Python script. If I simply translate a working parameter file (obtained from the ALPS website) into the below Python code, 'random()' is unrecognized (which is not so for the parameter file where the simulation goes through). But if I include 'import random' at the beginning of the Python script then the simulation does go through but a different model gets simulated (i.e. a random 'mu' but equal for all sites). I was wondering, should not the below Python script translated directly from a similar working parameter file run as is for the desired model? /*<SNIP>*/ import pyalps import matplotlib.pyplot as plt import pyalps.plot #prepare the input parameters parms = [] for DISORDERSEED in [34275, 49802, 82529]: parms.append( { 'LATTICE' : "inhomogeneous square lattice periodic", 'MODEL' : "boson Hubbard", 'T' : 0.1, 'L' : 4 , 't' : 1.0 , 'DISORDERSEED' : DISORDERSEED, 'mu' : 5.0*2*(random()-0.5), 'U' : 1.0 , 'NONLOCAL' : 0 , 'Nmax' : 2 , 'THERMALIZATION' : 100, 'SWEEPS' : 5000 } ) /*<SNIP>*/ Thanks, Vipin

2 3

how to run alps in Win7?
by Tianyou YI 06 Mar '14

06 Mar '14

Hello everyone, I have installed ALPS in Win7 and I also installed Pythonxy2.7.6.0, but I have not installed Vistrails because I cannot download it. But I do not know how to run ALPS, can anyone help me ? Thank you very much! Tianyou Yi

2 1