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

Quantum Monte Carlo for Ising in transverse field
by zhian asadzadeh 25 Jul '11

25 Jul '11

Dear Matthias, by \Sigma_{i,j}, I mean \sum_{i,j}, for the exact result, we know the model is integrable, so by fermionization and finally diagonalizing the Hamiltonian we will end up to an Hamiltonian with sum over independent modes, H=\sum_{k} E(k) (2 a_{k}^{dagger} a_{k}-1) for the T=0 ground state energy, I have compared with DMRG results which are exactly same, so I am sure up to this point the form of E(k) and my summation over k space is correct. so finally for E(T) simply we have: E(T)=\sum_{k} E(k)(2*n(k,T)-1) in which n(k,T)=1/(exp(E(k)/T)+1) ::with Boltzmann constant=1 Regards, Zhian. On Sun, Jul 24, 2011 at 12:00 PM, < comp-phys-alps-users-request(a)lists.phys.ethz.ch> wrote: > 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. Quantum Monte Carlo for Ising in transverse field > (zhian asadzadeh) > 2. Re: Quantum Monte Carlo for Ising in transverse field > (Matthias Troyer) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sat, 23 Jul 2011 16:29:39 +0200 > From: zhian asadzadeh <zhian.asadzadeh(a)gmail.com> > Subject: [ALPS-users] Quantum Monte Carlo for Ising in transverse > field > To: comp-phys-alps-users(a)lists.phys.ethz.ch > Message-ID: > <CAEw6G6sPavSO0RMreyy701xo+LN4M4LGrBzOLJkonEKs9S=Wkw(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Dear Matthias, > > my intention is doing finite T QMC for 2D ising in transverse field with > the following Hamiltonian: > H=-J \Sigma_{i,j} S_{i}^{z} S_{j}^{z}-h \Sigma_{i} S_{i}^{x} > but for checking the accuracy of the numeric I did some checks for 1D > case which we can have an exact results, > probably I am expecting too much accuracy, for example: > > L=100, T=0.32859 > Exact:: E(T)=-105.29789708270407, Mx(T)=0.27030367598168420 > looper:: E(T)=-105.61795 (0.00554), M(T)=0.4093185 (0.000138) ( I > think it is in z direction) > with SWEEPS:6000000, THERMALIZATION=2000000 > > I have checked also for other sizes and very large SWEEPS. > I don't expect to get exact result but I need much more accuracy in the > calculated energy. > > and also other question: that magnetization which looper calculates is > in the z direction? > can I also measure M_{x} just with declaring it in parm file? > > Thanks a lot, > Zhian. > > > On Sat, Jul 23, 2011 at 12:00 PM, < > comp-phys-alps-users-request(a)lists.phys.ethz.ch> wrote: > > > 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. restarting Lanczos iteration (jessica.alfonsi(a)unipd.it) > > 2. Quantum Monte Carlo for Ising in transvers field (zhian asadzadeh) > > 3. Re: Quantum Monte Carlo for Ising in transvers field > > (Matthias Troyer) > > 4. Re: restarting Lanczos iteration (Matthias Troyer) > > > > > > ---------------------------------------------------------------------- > > > > Message: 1 > > Date: Fri, 22 Jul 2011 12:22:47 +0200 > > From: jessica.alfonsi(a)unipd.it > > Subject: [ALPS-users] restarting Lanczos iteration > > To: comp-phys-alps-users(a)lists.phys.ethz.ch > > Message-ID: > > <5b0ae17b0813c63a97720e0fbc6cc1c6.squirrel(a)webmail.unipd.it> > > Content-Type: text/plain;charset=utf-8 > > > > Hi all, > > I'd like to modify the Lanczos *.cpp examples in IETL example folder to > be > > able to restart the iteration from any interruption point during the > > Lanczos iteration. I suppose first one has to dump some information, for > > instance the cofficients of T-matrix and the starting vector, however the > > more difficult part would be the resuming of the iteration. Can somebody > > give me some hints on how accomplishing this task? Thanks in advance... > > > > Best regards, > > > > Jessica Alfonsi > > > > > > > > > > > > ------------------------------ > > > > Message: 2 > > Date: Fri, 22 Jul 2011 13:03:23 +0200 > > From: zhian asadzadeh <zhian.asadzadeh(a)gmail.com> > > Subject: [ALPS-users] Quantum Monte Carlo for Ising in transvers field > > To: comp-phys-alps-users(a)lists.phys.ethz.ch > > Message-ID: > > < > CAEw6G6vLtP2G8kviZ97cxgZy9Dv1x_MqyUzFWs9ejnM8WvwnmA(a)mail.gmail.com > > > > > Content-Type: text/plain; charset="iso-8859-1" > > > > Dear All, > > > > I am trying to do QMC at finite temperature for calculating the > > magnetization of *Ising model in transverse field*, but I realized just > > looper code is working for this model, ( the others codes doesn't work > with > > off diagonal matrix elements) but the result are not so good. > > I am interested to modify dirloop_sse for this model, is it possible? > > I appreciate your help, > > Thanks, > > zhian. > >

3 2

Quantum Monte Carlo for Ising in transverse field
by zhian asadzadeh 25 Jul '11

25 Jul '11

Dear Matthias, my intention is doing finite T QMC for 2D ising in transverse field with the following Hamiltonian: H=-J \Sigma_{i,j} S_{i}^{z} S_{j}^{z}-h \Sigma_{i} S_{i}^{x} but for checking the accuracy of the numeric I did some checks for 1D case which we can have an exact results, probably I am expecting too much accuracy, for example: L=100, T=0.32859 Exact:: E(T)=-105.29789708270407, Mx(T)=0.27030367598168420 looper:: E(T)=-105.61795 (0.00554), M(T)=0.4093185 (0.000138) ( I think it is in z direction) with SWEEPS:6000000, THERMALIZATION=2000000 I have checked also for other sizes and very large SWEEPS. I don't expect to get exact result but I need much more accuracy in the calculated energy. and also other question: that magnetization which looper calculates is in the z direction? can I also measure M_{x} just with declaring it in parm file? Thanks a lot, Zhian. On Sat, Jul 23, 2011 at 12:00 PM, < comp-phys-alps-users-request(a)lists.phys.ethz.ch> wrote: > 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. restarting Lanczos iteration (jessica.alfonsi(a)unipd.it) > 2. Quantum Monte Carlo for Ising in transvers field (zhian asadzadeh) > 3. Re: Quantum Monte Carlo for Ising in transvers field > (Matthias Troyer) > 4. Re: restarting Lanczos iteration (Matthias Troyer) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Fri, 22 Jul 2011 12:22:47 +0200 > From: jessica.alfonsi(a)unipd.it > Subject: [ALPS-users] restarting Lanczos iteration > To: comp-phys-alps-users(a)lists.phys.ethz.ch > Message-ID: > <5b0ae17b0813c63a97720e0fbc6cc1c6.squirrel(a)webmail.unipd.it> > Content-Type: text/plain;charset=utf-8 > > Hi all, > I'd like to modify the Lanczos *.cpp examples in IETL example folder to be > able to restart the iteration from any interruption point during the > Lanczos iteration. I suppose first one has to dump some information, for > instance the cofficients of T-matrix and the starting vector, however the > more difficult part would be the resuming of the iteration. Can somebody > give me some hints on how accomplishing this task? Thanks in advance... > > Best regards, > > Jessica Alfonsi > > > > > > ------------------------------ > > Message: 2 > Date: Fri, 22 Jul 2011 13:03:23 +0200 > From: zhian asadzadeh <zhian.asadzadeh(a)gmail.com> > Subject: [ALPS-users] Quantum Monte Carlo for Ising in transvers field > To: comp-phys-alps-users(a)lists.phys.ethz.ch > Message-ID: > <CAEw6G6vLtP2G8kviZ97cxgZy9Dv1x_MqyUzFWs9ejnM8WvwnmA(a)mail.gmail.com > > > Content-Type: text/plain; charset="iso-8859-1" > > Dear All, > > I am trying to do QMC at finite temperature for calculating the > magnetization of *Ising model in transverse field*, but I realized just > looper code is working for this model, ( the others codes doesn't work with > off diagonal matrix elements) but the result are not so good. > I am interested to modify dirloop_sse for this model, is it possible? > I appreciate your help, > Thanks, > zhian. >

3 2

restarting Lanczos iteration
by jessica.alfonsi＠unipd.it 22 Jul '11

22 Jul '11

Hi all, I'd like to modify the Lanczos *.cpp examples in IETL example folder to be able to restart the iteration from any interruption point during the Lanczos iteration. I suppose first one has to dump some information, for instance the cofficients of T-matrix and the starting vector, however the more difficult part would be the resuming of the iteration. Can somebody give me some hints on how accomplishing this task? Thanks in advance... Best regards, Jessica Alfonsi

2 1

Quantum Monte Carlo for Ising in transvers field
by zhian asadzadeh 22 Jul '11

22 Jul '11

Dear All, I am trying to do QMC at finite temperature for calculating the magnetization of *Ising model in transverse field*, but I realized just looper code is working for this model, ( the others codes doesn't work with off diagonal matrix elements) but the result are not so good. I am interested to modify dirloop_sse for this model, is it possible? I appreciate your help, Thanks, zhian.

2 1

installation problem
by zhian asadzadeh 23 Jun '11

23 Jun '11

Dear All, I am trying to install *ALPS.2.01 *on cluster, for the installation I have followed the provided instruction on the main page of the project, those three step for installation has been down without error, but surprisingly finally I have nothing on the specified path for the installed code. these are the out put of each command: $ *cmake* -D Boost_ROOT_DIR=/u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/boost/ -D Boost_INCLUDE_DIR=/u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/boost/ -D FFTW_LIBRARIES=/opt/fftw/3.2.2/intel/lib -D FFTW_INCLUDE_DIR=/opt/fftw/3.2.2/intel/include -D MKL_INCLUDE_DIR=/opt/intel/mkl/10.0.5.025/include/ -D LAPACK_LIBRARY=/opt/intel/Compiler/11.1/069/mkl/lib/em64t/libmkl_core.so -D LAPACK_64_BIT=ON -D HDF5_INCLUDE_DIR=/opt/hdf5/1.8.5p1/gnu/include/ -D HDF5_LIBRARIES=/opt/hdf5/1.8.5p1/gnu/lib/libhdf5.so.6 -D BLAS_LIBRARY=/opt/rinave/blas-lapack/libblas.so.3.0.3 -D MPI_LIBRARY=/opt/openmpi/1.3.4/intel-11.0/lib -DCMAKE_INSTALL_PREFIX=/u/cm/asadzade/ALPS2/ /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/* -- Build type: Release -- Python interpreter /opt/python/2.6.4/bin/python -- Python interpreter ok : version 2.6.4 -- PYTHON_INCLUDE_DIRS = /opt/python/2.6.4/include/python2.6 -- PYTHON_NUMPY_INCLUDE_DIR = /opt/numpy/1.5.1/lib/python2.6/site-packages/numpy/core/include -- PYTHON_SITE_PKG = /opt/python/2.6.4/lib/python2.6/site-packages -- PYTHON_LIBRARY = /opt/python/2.6.4/lib/python2.6/config/libpython2.6.a -- PYTHON_EXTRA_LIBS =-lpthread -ldl -lutil -- PYTHON_LINK_FOR_SHARED = -Xlinker -export-dynamic -- ALPS version: 2.0.1 -- Looking for Boost Source -- Found Boost Source: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/boost/ -- Boost Version: 1_46_1 -- LPSolve Library: not found -- Using ALPS-installed ZLIB -- Using ALPS-installed SZIP -- Could NOT find SZIP (missing: SZIP_LIBRARIES SZIP_INCLUDE_DIRS) -- Using ALPS-installed HDF5 in . -- HDF5 version found is 1.8.5-patch1 in /opt/hdf5/1.8.5p1/gnu/include -- Python interpreter /opt/python/2.6.4/bin/python -- Python interpreter ok : version 2.6.4 -- PYTHON_INCLUDE_DIRS = /opt/python/2.6.4/include/python2.6 -- PYTHON_NUMPY_INCLUDE_DIR = /opt/numpy/1.5.1/lib/python2.6/site-packages/numpy/core/include -- PYTHON_SITE_PKG = /opt/python/2.6.4/lib/python2.6/site-packages -- PYTHON_LIBRARY = /opt/python/2.6.4/lib/python2.6/config/libpython2.6.a -- PYTHON_EXTRA_LIBS =-lpthread -ldl -lutil -- PYTHON_LINK_FOR_SHARED = -Xlinker -export-dynamic -- Numpy include in /opt/numpy/1.5.1/lib/python2.6/site-packages/numpy/core/include -- ALPS XML dir is ./u/cm/asadzade/ALPS2/xml -- Adding Boost dir: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/boost/ -- Configuring done -- Generating done -- Build files have been written to: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps *$* make[ 73%] Built target boost [ 75%] Built target boost_prg_exec_monitor [ 87%] Built target boost_test_exec_monitor [100%] Built target boost_unit_test_framework $ make install (just some part of it) [ 73%] Built target boost [ 75%] Built target boost_prg_exec_monitor [ 87%] Built target boost_test_exec_monitor [100%] Built target boost_unit_test_framework Install the project... -- Install configuration: "Release" -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/Welcome.txt -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/CITATIONS.txt -- Installing: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/CMakeCache.txt -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/LICENSE-package.txt -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/LICENSE.txt -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/LICENSE_1_0.txt -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/README-light.txt -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/LICENSE-applications.txt -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/LICENSE-light.txt -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/README-package.txt -- Installing: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/install_manifest.txt -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/CMakeLists.txt -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/README.txt -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/script/hdf5.sh -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565 . . .-- Installing: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/tutorials/./dmft-02-hybridization/beta_10_U3_tsqrt2 -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/tutorials/./dmft-02-hybridization/beta_10_U3_tsqrt2/G0_omega_input -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/tutorials/./dmft-02-hybridization/beta_10_U3_tsqrt2/hybrid.parm -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/tutorials/./dmft-02-hybridization/G0_tau_input_beta_12.0 -- Installing: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/tutorials/./dmft-02-hybridization/beta_12_U3_tsqrt2 -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/tutorials/./dmft-02-hybridization/beta_12_U3_tsqrt2/G0_omega_input -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/tutorials/./dmft-02-hybridization/beta_12_U3_tsqrt2/hybrid.parm -- Installing: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/tutorials/./dmft-02-hybridization/beta_6_U3_tsqrt2 -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/tutorials/./dmft-02-hybridization/beta_6_U3_tsqrt2/G0_omega_input -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/tutorials/./dmft-02-hybridization/beta_6_U3_tsqrt2/hybrid.parm -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/tutorials/./dmft-02-hybridization/dmft-02-inputs.zip -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/tutorials/./dmft-02-hybridization/G0_tau_input_beta_14.0 -- Installing: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/tutorials/./dmft-02-hybridization/beta_14_U3_tsqrt2 -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/tutorials/./dmft-02-hybridization/beta_14_U3_tsqrt2/G0_omega_input -- Up-to-date: /u/cm/asadzade/alps-2.0.1-r5565-src-with-boost/alps/./alps/tutorials/./dmft-02-hybridization/beta_14_U3_tsqrt2/hybrid.parm *could you please help me to know which part is not complete, Thanks, zhian*. *

2 1

Re: [ALPS-users] tebd will not be built (Matthias Troyer)
by Cristian Degli Esposti Boschi 03 Jun '11

03 Jun '11

> Did you specify a Fortran-90 compiler? I guess so, at least I installed gfortran (4.4) via the Debian procedure. However, now I am worried about these two things: 1) (from the installation guidelines on the ALPS website) The tebd code requires a fortran compiler and that HDF5 be built with the --enable-fortran configure option. Will this be the case with Debian's gfortran? How to check a posteriori in HDF5 in my machine has been built with the --enable-fortran option? 2) what is the role of the two scripts $ $HOME/src/alps2/script/cmake.sh $HOME/opt $HOME/tmp $ $HOME/src/alps2/script/hdf5.sh $HOME/opt $HOME/tmp if all the dependent libraries are already ok, are these scripts necessary? Thanks again for your help. Cristian On 05/26/2011 12:00 PM, comp-phys-alps-users-request(a)lists.phys.ethz.ch wrote: > 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. tebd will not be built (Cristian Degli Esposti Boschi) > 2. Re: tebd will not be built (Matthias Troyer) > 3. transverse Ising model with DMRG (shuai cui) -- 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

Define two-site spin correlation and controll the output in DMRG simulation
by 潘丽华 02 Jun '11

02 Jun '11

Dear ALPS users and developers, I am now studying the properties of the frustated Heisenberg spin chain. And the correlation functions is considered. The calculation result give (N+1)^2 terms of the two-sites' correlation value. Maybe we are interested in the result about C(l)=<Sz(N+2-l/2)Sz(N/2+l/2)>. I'd like to ask for your help. Can we modify the file /ALPS/alps-application/dmrg/dmrg/dmrg.c to get just the <Sz(N+2-l/2)Sz(N/2+l/2)> result? Or if there is another way to do so? The main reason we do not like the available one is that we find hat the output data of a single program is huge. So I'd like to consult you that if there's some other method to reduce the output data of the program? Thank you in advance for your time and effort to response. Best Regards! Your Sincerely, Pan Lihua

1 0

Modify the source code(dmrg.C)
by muzi352 30 May '11

30 May '11

Dear editor: i always use ALPS (DMRG) to calculate correlation function, my topics: when i calculate correlation chain long L=300, the block_ALPS_.dat is very large. how to chang the dmrg.C in order to calculate what i need corrletion,for example <S0Sr> <SiSi+1> thank you

1 0

Re: [ALPS-users] Problem of installing ALPS on Ubuntu: Falling back to CMake provied LAPACK/BLAS detection
by Brains Chen 27 May '11

27 May '11

> To whom it may concern, > > I need to install the ALPS 2.0 in a Linux Ubuntu-10.04 system. So far I get > an error message when using > > cmake -D > Boost_ROOT_DIR=/opt/local/alps-2.0/alps-2.0.1-r5565-src-with-boost/boost/ > -DCMAKE_INSTALL_PREFIX=/opt/local/alps-2.0/ > /opt/local/alps-2.0/alps-2.0.1-r5565-src-with-boost/alps > > -- Build type: Release > -- Python interpreter /usr/bin/python > -- Python interpreter ok : version 2.6.5 > -- PYTHON_INCLUDE_DIRS = /usr/include/python2.6 > -- PYTHON_NUMPY_INCLUDE_DIR = > /usr/lib/python2.6/dist-packages/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 > -- ALPS version: 2.0.1 > -- Looking for Boost Source > -- Found Boost Source: > /opt/local/alps-2.0/alps-2.0.1-r5565-src-with-boost/boost/ > -- Boost Version: 1_46_1 > 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. > -- LPSolve Library: not found > -- Could NOT find SZIP (missing: SZIP_LIBRARIES SZIP_INCLUDE_DIRS) > -- HDF5 version found is 1.8.4 in /usr/include > -- Python interpreter /usr/bin/python > -- Python interpreter ok : version 2.6.5 > -- PYTHON_INCLUDE_DIRS = /usr/include/python2.6 > -- PYTHON_NUMPY_INCLUDE_DIR = > /usr/lib/python2.6/dist-packages/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/python2.6/dist-packages/numpy/core/include > -- ALPS XML dir is /opt/local/alps-2.0/lib/xml > -- Adding Boost dir: > /opt/local/alps-2.0/alps-2.0.1-r5565-src-with-boost/boost/ > -- HDF5 version found is 1.8.4 in /usr/include > -- HDF5 version found is 1.8.4 in /usr/include > -- tebd will not be built > -- Configuring done > -- Generating done > -- Build files have been written to: > /opt/local/alps-2.0/alps-2.0.1-r5565-src-with-boost > > Since Dr. Matthias Troyer told me that it in not required to have the > lp_solve library for ALPS 2.0, I still have one error for Lapack (Falling > back to CMake provied LAPACK/BLAS detection). If I don't fix it and try to > install it, I will get an error message as > > ...... > SSE.Directed.cpp:(.text+0x2745): undefined reference to `add_constraint' > CMakeFiles/dirloop_sse_v1.dir/SSE.Directed.cpp.o:SSE.Directed.cpp:(.text+0x27a1): > more undefined references to `add_constraint' follow > CMakeFiles/dirloop_sse_v1.dir/SSE.Directed.cpp.o: In function > `SSE::find_minbounce(double*, double*, double*, double*)': > SSE.Directed.cpp:(.text+0x29cd): undefined reference to `set_obj_fn' > SSE.Directed.cpp:(.text+0x29d5): undefined reference to `solve' > SSE.Directed.cpp:(.text+0x2ae0): undefined reference to `delete_lp' > SSE.Directed.cpp:(.text+0x2b57): undefined reference to `add_constraint' > SSE.Directed.cpp:(.text+0x2bb6): undefined reference to `add_constraint' > SSE.Directed.cpp:(.text+0x2c74): undefined reference to `add_constraint' > SSE.Directed.cpp:(.text+0x2cd3): undefined reference to `add_constraint' > SSE.Directed.cpp:(.text+0x2d64): undefined reference to `add_constraint' > CMakeFiles/dirloop_sse_v1.dir/SSE.Directed.cpp.o:SSE.Directed.cpp:(.text+0x2dc0): > more undefined references to `add_constraint' follow > CMakeFiles/dirloop_sse_v1.dir/SSE.Directed.cpp.o: In function > `SSE::find_minbounce(double*, double*, double*, double*)': > SSE.Directed.cpp:(.text+0x2fe9): undefined reference to `set_upbo' > SSE.Directed.cpp:(.text+0x3008): undefined reference to `set_upbo' > SSE.Directed.cpp:(.text+0x3027): undefined reference to `set_upbo' > SSE.Directed.cpp:(.text+0x3040): undefined reference to `set_upbo' > SSE.Directed.cpp:(.text+0x305a): undefined reference to `set_upbo' > CMakeFiles/dirloop_sse_v1.dir/SSE.Directed.cpp.o:SSE.Directed.cpp:(.text+0x3073): > more undefined references to `set_upbo' follow > CMakeFiles/dirloop_sse_v1.dir/SSE.Directed.cpp.o: In function > `SSE::find_minbounce(double*, double*, double*, double*)': > SSE.Directed.cpp:(.text+0x312f): undefined reference to `add_constraint' > SSE.Directed.cpp:(.text+0x3189): undefined reference to `add_constraint' > SSE.Directed.cpp:(.text+0x31ff): undefined reference to `add_constraint' > SSE.Directed.cpp:(.text+0x325d): undefined reference to `add_constraint' > SSE.Directed.cpp:(.text+0x32de): undefined reference to `add_constraint' > CMakeFiles/dirloop_sse_v1.dir/SSE.Directed.cpp.o:SSE.Directed.cpp:(.text+0x333c): > more undefined references to `add_constraint' follow > collect2: ld returned 1 exit status > make[2]: *** [applications/qmc/sse/dirloop_sse_v1] Error 1 > make[1]: *** [applications/qmc/sse/CMakeFiles/dirloop_sse_v1.dir/all] Error > 2 > make[1]: *** Waiting for unfinished jobs.... > > It sounds like some dynamical libraries are not correct. Do you know how to > fix it? > > By the way, do "tebd will not be built" and "Could NOT find SZIP" matter? > > Thank you, > -- > Kuang-Shing Chen > > -- Kuang-Shing Chen

5 11

Re: [ALPS-users] Comp-phys-alps-users Digest, Vol 62, Issue 13
by zhian asadzadeh 26 May '11

26 May '11

Hi shuai cui, your model is correct, and the result for the second run is fair but not so near to exact result, I suggest you increase the number of *maxstates (=300) *, and also your *lattice size*, you will get a better result, just for check point I give you the exact energy (per site) at h=1 (which is critical h in this model) J=-1, h=1 L=100, *E0=-1.2732919* ( This is exact energy not DMRG result) Zhian. On Thu, May 26, 2011 at 12:00 PM, < comp-phys-alps-users-request(a)lists.phys.ethz.ch> wrote: > 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. tebd will not be built (Cristian Degli Esposti Boschi) > 2. Re: tebd will not be built (Matthias Troyer) > 3. transverse Ising model with DMRG (shuai cui) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 25 May 2011 13:16:52 +0200 > From: Cristian Degli Esposti Boschi <degliesposti(a)bo.imm.cnr.it> > Subject: [ALPS-users] tebd will not be built > To: comp-phys-alps-users(a)lists.phys.ethz.ch > Message-ID: <4DDCE524.6070605(a)bo.imm.cnr.it> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > 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 > > > ------------------------------ > > Message: 2 > Date: Wed, 25 May 2011 06:51:39 -1000 > From: Matthias Troyer <troyer(a)phys.ethz.ch> > Subject: Re: [ALPS-users] tebd will not be built > To: comp-phys-alps-users(a)lists.phys.ethz.ch > Message-ID: <291CA21C-BF47-48EC-881E-27DD21C9F49D(a)phys.ethz.ch> > Content-Type: text/plain; charset=us-ascii > > Did you specify a Fortran-90 compiler? > > On May 25, 2011, at 1:16 AM, Cristian Degli Esposti Boschi wrote: > > > 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 > > > > > > ------------------------------ > > Message: 3 > Date: Thu, 26 May 2011 11:28:42 +0800 (CST) > From: shuai cui <cuishuai_757(a)yahoo.com.cn> > Subject: [ALPS-users] transverse Ising model with DMRG > To: comp-phys-alps-users(a)lists.phys.ethz.ch > Message-ID: <606168.86645.qm(a)web15904.mail.cnb.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > 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 > ? >

2 1