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

Re: [ALPS-users] Comp-phys-alps-users Digest, Vol 112, Issue 12
by Tadeusz Wasiutyński 27 Jul '15

27 Jul '15

I had about the same and the solution was given by Synge in mail from May 1st: Message: 1 Date: Thu, 30 Apr 2015 19:53:12 +0900 From: Synge Todo <wistaria(a)comp-phys.org> To: comp-phys-alps-users(a)lists.phys.ethz.ch Subject: Re: [ALPS-users] building alps from source Message-ID: <B8D338DA-8CC7-4C6A-9E25-C19CACE3C427(a)comp-phys.org> Content-Type: text/plain; charset=utf-8 This looks similar as the problem reported at https://svn.boost.org/trac/boost/ticket/8973 . You can solve the problem by using a newer version of Boost, or apply the patch provided at https://svn.boost.org/trac/boost/changeset/84950 to the Boost source code. Best, Synge 2015-07-27 12:00 GMT+02:00 <comp-phys-alps-users-request(a)lists.phys.ethz.ch> : > Send Comp-phys-alps-users mailing list submissions to > comp-phys-alps-users(a)lists.phys.ethz.ch > > To subscribe or unsubscribe via the World Wide Web, visit > https://lists.phys.ethz.ch/listinfo/comp-phys-alps-users > or, via email, send a message with subject or body 'help' to > comp-phys-alps-users-request(a)lists.phys.ethz.ch > > You can reach the person managing the list at > comp-phys-alps-users-owner(a)lists.phys.ethz.ch > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of Comp-phys-alps-users digest..." > > > Today's Topics: > > 1. make failed (Inzo Babaria) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Mon, 27 Jul 2015 11:31:42 +0530 > From: Inzo Babaria <inzo.babaria(a)gmail.com> > To: comp-phys-alps-users(a)lists.phys.ethz.ch > Subject: [ALPS-users] make failed > Message-ID: > < > CACz9evC+dK46UysXmVSo1VLGwkxt1BFFa2p-CwQuCaETWN3dSA(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Dear all users and developers, > > I was building and installing ALPS in my loptop. I have done each process > step by step. In the last step I created a build Folder and I unziped > alps-2.2.b3-r7462-src-with-boost. therefore I went to Build folder > > cd build > > and wrote the bellow: > > cmake -D > > Boost_ROOT_DIR:PATH=/home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost > /home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/alps > > After it finished, I typed make > > acerv3-571g@acerv3571g-Aspire-V3-571G:~/Downloads/build$ make > > it was going ok but I faced to these error message: > > [ 21%] Building CXX object > > src/boost/CMakeFiles/boost.dir/home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/libs/mpi/src/request.cpp.o > [ 21%] Building CXX object > > src/boost/CMakeFiles/boost.dir/home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/libs/mpi/src/text_skeleton_oarchive.cpp.o > [ 21%] Building CXX object > > src/boost/CMakeFiles/boost.dir/home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/libs/mpi/src/timer.cpp.o > [ 21%] Building CXX object > > src/boost/CMakeFiles/boost.dir/home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/libs/thread/src/future.cpp.o > [ 21%] Building CXX object > > src/boost/CMakeFiles/boost.dir/home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/libs/thread/src/pthread/once.cpp.o > In file included from > > /home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/boost/atomic.hpp:12:0, > from > > /home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/boost/thread/pthread/once_atomic.hpp:20, > from > > /home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/boost/thread/once.hpp:20, > from > > /home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/libs/thread/src/pthread/./once_atomic.cpp:9, > from > > /home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/libs/thread/src/pthread/once.cpp:8: > > /home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/boost/atomic/atomic.hpp:202:16: > error: ?uintptr_t? was not declared in this scope > typedef atomic<uintptr_t> atomic_uintptr_t; > ^ > > /home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/boost/atomic/atomic.hpp:202:25: > error: template argument 1 is invalid > typedef atomic<uintptr_t> atomic_uintptr_t; > ^ > > /home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/boost/atomic/atomic.hpp:202:43: > error: invalid type in declaration before ?;? token > typedef atomic<uintptr_t> atomic_uintptr_t; > ^ > src/boost/CMakeFiles/boost.dir/build.make:3159: recipe for target > > 'src/boost/CMakeFiles/boost.dir/home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/libs/thread/src/pthread/once.cpp.o' > failed > make[2]: *** > > [src/boost/CMakeFiles/boost.dir/home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/libs/thread/src/pthread/once.cpp.o] > Error 1 > CMakeFiles/Makefile2:966: recipe for target > 'src/boost/CMakeFiles/boost.dir/all' failed > make[1]: *** [src/boost/CMakeFiles/boost.dir/all] Error 2 > Makefile:147: recipe for target 'all' failed > make: *** [all] Error 2 > acerv3-571g@acerv3571g-Aspire-V3-571G:~/Downloads/build$ > > I don't know what is the problem and how I can remove that?! Thanks so much > to any body for getting help and guidance. > > Inzo >

1 0

make failed
by Inzo Babaria 27 Jul '15

27 Jul '15

Dear all users and developers, I was building and installing ALPS in my loptop. I have done each process step by step. In the last step I created a build Folder and I unziped alps-2.2.b3-r7462-src-with-boost. therefore I went to Build folder cd build and wrote the bellow: cmake -D Boost_ROOT_DIR:PATH=/home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost /home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/alps After it finished, I typed make acerv3-571g@acerv3571g-Aspire-V3-571G:~/Downloads/build$ make it was going ok but I faced to these error message: [ 21%] Building CXX object src/boost/CMakeFiles/boost.dir/home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/libs/mpi/src/request.cpp.o [ 21%] Building CXX object src/boost/CMakeFiles/boost.dir/home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/libs/mpi/src/text_skeleton_oarchive.cpp.o [ 21%] Building CXX object src/boost/CMakeFiles/boost.dir/home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/libs/mpi/src/timer.cpp.o [ 21%] Building CXX object src/boost/CMakeFiles/boost.dir/home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/libs/thread/src/future.cpp.o [ 21%] Building CXX object src/boost/CMakeFiles/boost.dir/home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/libs/thread/src/pthread/once.cpp.o In file included from /home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/boost/atomic.hpp:12:0, from /home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/boost/thread/pthread/once_atomic.hpp:20, from /home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/boost/thread/once.hpp:20, from /home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/libs/thread/src/pthread/./once_atomic.cpp:9, from /home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/libs/thread/src/pthread/once.cpp:8: /home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/boost/atomic/atomic.hpp:202:16: error: ‘uintptr_t’ was not declared in this scope typedef atomic<uintptr_t> atomic_uintptr_t; ^ /home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/boost/atomic/atomic.hpp:202:25: error: template argument 1 is invalid typedef atomic<uintptr_t> atomic_uintptr_t; ^ /home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/boost/atomic/atomic.hpp:202:43: error: invalid type in declaration before ‘;’ token typedef atomic<uintptr_t> atomic_uintptr_t; ^ src/boost/CMakeFiles/boost.dir/build.make:3159: recipe for target 'src/boost/CMakeFiles/boost.dir/home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/libs/thread/src/pthread/once.cpp.o' failed make[2]: *** [src/boost/CMakeFiles/boost.dir/home/acerv3-571g/Downloads/alps-2.2.b3-r7462-src-with-boost/boost/libs/thread/src/pthread/once.cpp.o] Error 1 CMakeFiles/Makefile2:966: recipe for target 'src/boost/CMakeFiles/boost.dir/all' failed make[1]: *** [src/boost/CMakeFiles/boost.dir/all] Error 2 Makefile:147: recipe for target 'all' failed make: *** [all] Error 2 acerv3-571g@acerv3571g-Aspire-V3-571G:~/Downloads/build$ I don't know what is the problem and how I can remove that?! Thanks so much to any body for getting help and guidance. Inzo

1 0

mps_optim for inhomogeneous Bose-Hubbard
by Dominique Delande 23 Jul '15

23 Jul '15

Dear all, I try to use mps_optim for the 1D Bose-Hubbard model in the presence of a trapping potential. I use the latest SVN version (from last week), the latticel model "inhomogeneous chain lattice" and a homemade model "harm boson Hubbard" with the site term: -mu#*n(i)+U#*n(i)*(n(i)-1)/2+K#*n(i)*(x-0.5*(L-1))^2 which differs simply from the usual Bose-Hubbard model by the addition of the last term K#*n(i)*(x-0.5*(L-1))^2. Such a model works with the worm QMC algorithm, but fails with mps_optim with the error message: Exception thrown: Cannot evaluate x I tried to replace "x" by "i" without success. How should I proceed to have mps_optim work for an inhomogeneous system? More generally, I find it always **VERY** difficult to obtain the information on what should work for a specific algorithm/model in ALPS and what should not work. The only solution I found is trial and error, which is a bit frustrating. Thanks in advance for your answer(s). Dominique -- Dominique Delande (Dominique.Delande(a)lkb.upmc.fr) Laboratoire Kastler-Brossel - Case 74 - Universite P. et M. Curie 4, place Jussieu, F-75252 Paris Cedex 05, FRANCE Phone : +33 1 44 27 27 97 - Fax : +33 1 44 27 38 45 Acces : Tour 13, Couloir 12-13, 3eme etage - Bureau 316

2 1

follow up to problems with dwa
by Jakub Zakrzewski 23 Jul '15

23 Jul '15

Dear all, As another example of using dwa-algorithm I tried to run 3D case in trap. The relevant part of "model.xml" is: <SITETERM site="i"> <PARAMETER name="mu#" default="mu"/> <PARAMETER name="U#" default="U"/> <PARAMETER name="K#" default="K"/> -mu#*n(i)+U#*n(i)*(n(i)-1)/2+K#*n(i)*((x-(L-1)/2)^2+(y-(L-1)/2)^2+(z-(L-1)/2)^2) </SITETERM> which works with both worm and dwa. But suppose I want a QUARTIC trap and not harmonic so I modify to <SITETERM site="i"> <PARAMETER name="mu#" default="mu"/> <PARAMETER name="U#" default="U"/> <PARAMETER name="K#" default="K"/> -mu#*n(i)+U#*n(i)*(n(i)-1)/2+K#*n(i)*((x-(L-1)/2)^2+(y-(L-1)/2)^2+(z-(L-1)/2)^2)^2 </SITETERM> Note that the only difference is a power 2 on squared radius. Worm manages fine. Dwa gives error message "Cannot evaluate y" What I am doing wrong? Is the failure of dwa due to an additional power 2 (since "dwa" in polish means "two")? I know the answer - dwa fails for any power different than 1.... Help appreciated best kuba Jakub Zakrzewski Zakład Optyki Atomowej, Instytut Fizyki im. Mariana Smoluchowskiego, UJ ul. Lojasiewicza 11, 30-348 Kraków, Poland +48-12-6644555, fax +48126644905 http://chaos.if.uj.edu.pl/~kuba

2 1

worm versus dwa
by Jakub Zakrzewski 23 Jul '15

23 Jul '15

Dear All, I use worm with a harmonic trap defined in a standard version as in ALPS2.0 paper with the following addition in <SITETERM> (for 2D vesion) K#*n(i)*((x-(L-1)/2)^2+(y-(L-1)/2)^2) works perfectly. I tried dwa getting the error message: cannot evaluate y Trying the one dimensional version with just K#*n(i)*(x-(L-1)/2)^2 works for both worm and dwa. any idea why 2D harmonic trap cannot be used with dwa? All the best kuba Jakub Zakrzewski Zakład Optyki Atomowej, Instytut Fizyki im. Mariana Smoluchowskiego, UJ ul. Lojasiewicza 11, 30-348 Kraków, Poland +48-12-6644555, fax +48126644905 http://chaos.if.uj.edu.pl/~kuba

1 0

Error when using 'singlesite' optimization and 'storagedir' together for MPS based DMRG
by Wei.Su 22 Jul '15

22 Jul '15

I occured an error "Intel MKL ERROR: Parameter 10 was incorrect on entry to DGEMM " when using 'singlesite' optimization and 'storagedir' together. I'm not sure if this is an error of MKL or one of the MPS code in ALPS However it runs ok if i turn off the "storagedir" option or use the "twosite" optimization schedule. Here comes the python script #--------------------------------------------------------------------------------------------# import pyalps import numpy as np import matplotlib.pyplot as plt import pyalps.plot #prepare the input parameters parms = [ { 'optimization' : 'singlesite', 'LATTICE' : 'open chain lattice', 'L' : 20, 'MODEL' : 'spin', 'local_S0' : '0.5', 'local_S1' : '1', 'CONSERVED_QUANTUMNUMBERS' : 'N,Sz', 'Sz_total' : 0, 'J' : 1, 'SWEEPS' : 4, 'NUMBER_EIGENVALUES' : 1, 'MAXSTATES' : 50, 'storagedir' : 'storage' } ] #write the input file and run the simulation input_file = pyalps.writeInputFiles('parm_spin_one',parms) res = pyalps.runApplication('mps_optim',input_file,writexml=True) #----------------------------------------------------------------------------------------# -------Su Wei-------

3 2

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

Acquiring data to susceptibility for Ising`s Model
by Felipe da Cruz Rodrigues 17 Jul '15

17 Jul '15

Hello all. Someone know how can I to acquire the data for susceptibility in Ising's Model (ex: tutorial 7)? p.s.: I don't want the generated graphic for the susceptibility in vistrails, but only the data. Thanks a lot! -- Felipe da Cruz Rodrigues Mestrando em Física pela Universidade Federal de Lavras Graduado em Física (Licenciatura) - UFLA

2 3

mps_optim example 1: itinerant fermions in the Hubbard ladder
by Dominique Delande 09 Jul '15

09 Jul '15

Dear all, I am trying to use mps_optim and have serious problems. I try to run Example 1 (itinerant fermions in the Hubbard ladder) presented in M. Dolfi et al, CPC 185, 3430 (2014). I downloaded the mmcs1.zip file which should contain all scripts and results. When I run the python script for Example 1, my results are quite different from what they should be: * No "EnergyVariance" variable in the result files sim.taskX.out.h5 * All "pair field" results are zero, or vanishingly small (of the order of 1.e-300). * In the result files extracted from mmcs1.zip (e.g. sim.task1.out.h5), there are additional results such as "Energy^2", "Local density down", "dens corr down-down", etc. * In the result files extracted from mmcs1.zip (e.g. sim.task1.out.h5), in the "parameters" section, there are things like MEASURE_HALF_NN_CORRELATIONS which are not defined in the simulation parameters (most probably it it the same as MEASURE_HALF_CORRELATIONS). I thus have the impression that the results files provided in mmcs1.zip have NOT been obtained using the provided python scripts. More importantly, how to obtain correct "EnergyVariance" and "pair field" results ? I am using alps-2.2.b3-r7462, compiled from source on a Fedora 21 x86_64 machine, using the GNU compilers and Intel MKL libraries. Thanks in advance for your help. Dominique -- Dominique Delande (Dominique.Delande(a)lkb.upmc.fr) Laboratoire Kastler-Brossel - Case 74 - Universite P. et M. Curie 4, place Jussieu, F-75252 Paris Cedex 05, FRANCE Phone : +33 1 44 27 27 97 - Fax : +33 1 44 27 38 45 Acces : Tour 13, Couloir 12-13, 3eme etage - Bureau 316

2 1

Error when searching for excited states with MPS code
by Wei.Su 08 Jul '15

08 Jul '15

I was using MPS code in ALPS to calculate exited states of the fermion Hubbard model when there throwed an error -------------------------------------------------------------------------------------------- Sweep 0, optimizing sites 19 and 20 Input <MPS|O[0]> : 0 Output <MPS|O[0]> : -nan JCD used 8 iterations. Time elapsed in JCD: 0.19799719 MPS overlap: 0 Energy 1 2.341151156e-310 Exception thrown: Error in SVD! -------------------------------------------------------------------------------------------------- This error vanished when the "total number" was less than 10 or great than 20 (or around) or the 'MAXSTATES' is greater than about 80. In the attachment I deliver my python script. -------------------------------------------------------------------------------------------------------------- import pyalps import numpy as np parms=[] for L in [40]: for P in [12]: # for P in [55,56,57,58,59,60,61,62]: parms.append({ 'LATTICE' : "open chain lattice", 'L' : L, 'MODEL' : "fermion Hubbard", 'CONSERVED_QUANTUMNUMBERS' : 'Nup,Ndown', 'Nup_total' : P, 'Ndown_total' : P, 'mu' : 0.0, 't' : 1.0, 'U' : 0.0, 'SWEEPS' : 4, 'MAXSTATES' : 50, 'NUMBER_EIGENVALUES' : 2, 'storagedir' : 'storage', }) #write the input file and run the simulation input_file = pyalps.writeInputFiles('parm_tsc',parms) res = pyalps.runApplication('mps_optim',input_file,writexml=True) -------------------------------------------------------------------------------------------------------------------- -- -------苏威-------

2 1