March 2012 - 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 2.0.1 released
by Matthias Troyer 19 Sep '12

19 Sep '12

Dear ALPS users, ALPS 2.0.1 has been released. This is the version included with the ALPS paper: B. Bauer et al., J. Stat. Mech. (2011) P05001, which is available at : http://www.iop.org/EJ/abstract/1742-5468/2011/05/P05001 .Besides references to this paper, this release contains many minor bug fixes. The MacOS X 10.5 binary release will be uploaded shortly. The others are ready for download. Best regards Matthias Troyer

4 33

help
by #ABDUL MALEK BIN RAZALI# 10 Apr '12

10 Apr '12

Hi, I am using the honeycomb lattice and boson hubbard model from the alps library on vispython. I wish to measure the staggered structure factor given by the equation S = 1/N * Summation [ e(i)e(j) <n(i)n(j)> ] where e(i) = 1 or -1 for i on sublattice A(B) in terms of density-density correlations. I know that alps measure the density structure factor using wave vectors. How can i make alps measure the structure factor given above ? Regards, Mike

2 7

tests failed
by Rajeev Singh 22 Mar '12

22 Mar '12

Hi, I have installed Alps on a computer running SUSE Linux Enterprise Server 10 (x86_64) using Intel Compilers. The following is the result of 'make test' The following tests FAILED: 8 - hello_fortran (Failed) 9 - ising_fortran (Failed) 42 - lattice_constant_test (Failed) 80 - h5d_cpx_flt_ (Failed) 81 - h5dc_cpx_flt_ (Failed) 82 - h5a_cpx_flt_ (Failed) 83 - h5ac_cpx_flt_ (Failed) 88 - h5d_cpx_lng_dbl_ (Failed) 89 - h5dc_cpx_lng_dbl_ (Failed) 90 - h5a_cpx_lng_dbl_ (Failed) 91 - h5ac_cpx_lng_dbl_ (Failed) 536 - h5d_pr_dbl_cpx_dbl_ (Failed) 537 - h5dc_pr_dbl_cpx_dbl_ (Failed) 538 - h5a_pr_dbl_cpx_dbl_ (Failed) 539 - h5ac_pr_dbl_cpx_dbl_ (Failed) 612 - python_hlist_test (Failed) 613 - python_mcdata (Failed) 614 - python_pyhdf5 (Failed) Errors while running CTest Which of these should I worry about? All the programs in the tutorial seem to run fine. Regards, Rajeev

3 5

Re: [ALPS-users] Transverse magnetisation measurement
by Matthias Troyer 19 Mar '12

19 Mar '12

Hi Simon, with SU(2) symmetry the mean transverse magnetization is always 0 by definition. Matthias On 19 Mar 2012, at 09:14, Simon Ward wrote: > Hi, > > I'm using the spin model and SU(2) symmetry is currently preserved. > > Regards, > > Simon > > > On Mon, 2012-03-19 at 07:31 -0700, Matthias Troyer wrote: >> Hi Simon >> >> What is your model? Is it SU(2)symmetric? >> >> Matthias >> >> >> >> >> On Mar 19, 2012, at 2:55, "Ward Simon" <Simon.Ward(a)psi.ch> wrote: >> >>> Hello all, >>> >>> I have defined a set of S=1/2 3D lattices and would like to measure the transverse magnetisation (mx) in a longitudinal field using the dirloop_sse application. I presume a modification to the source/model files are needed to achieve this but I currently do not know where to start. >>> >>> Many thanks for your help and time in advance, >>> >>> S Ward >>> >

1 0

Re: [ALPS-users] Transverse magnetisation measurement
by Matthias Troyer 19 Mar '12

19 Mar '12

Hi Simon What is your model? Is it SU(2)symmetric? Matthias On Mar 19, 2012, at 2:55, "Ward Simon" <Simon.Ward(a)psi.ch> wrote: > Hello all, > > I have defined a set of S=1/2 3D lattices and would like to measure the transverse magnetisation (mx) in a longitudinal field using the dirloop_sse application. I presume a modification to the source/model files are needed to achieve this but I currently do not know where to start. > > Many thanks for your help and time in advance, > > S Ward >

1 0

spimc questions
by daniel.aravena＠antares.qi.ub.edu 14 Mar '12

14 Mar '12

Dear all, I am performing some tests using spinmc (alps 2.0, Heisenberg model) on uncoupled spins changing the number and magnitude of spins at different fields. There are some things in the results that I dont understand: The energy is always zero, independently of the magnitude of the field for the uncoupled case, if I put a coupling constant, then I get as expected a non-zero value for the energy, but it is still non-dependent on field as it should be because of the Zeeman term. In a previous post, it is said that the susceptibility is calculated from magnetization^2, so it should work while M is linear on H, is there a way to get the dependence of susceptibility against field as the system approaches to saturation? And regarding the direction of field, setting H=1 means: Hz=1 Hx, Hy and Hz=1 Hx^2+Hy^2+Hz^2=1^2 , Hx=Hy=Hz or something else? I also tried to set H="0 0 1" and the program runs, but i dont know what I am actually getting. Any comments are welcome Kind Regards Daniel

3 14

help
by #ABDUL MALEK BIN RAZALI# 12 Mar '12

12 Mar '12

Hi, I would like to ask about the definition of the honeycomb lattice in the lib. Firstly, the honeycomb unitcell has 2 sites at positions (0,0) and (-0.28866...,0.57734...). I do not understand the second coordinate. Secondly, is it necessary to have different edge types although they represent the same bond connecting nearest neighbour? Regards, Mike

1 0

Re: [ALPS-users] Problems compiling alps-2.0.3-src-with-boost (Andreas Sindermann)
by Tadeusz Wasiutyński 08 Mar '12

08 Mar '12

Andreas, I have the same message after make on Scientific Linux 6.1, nevertheless make install goes smoothly and I did not notice any problems going through tutorial. Tadeusz 2012/2/16 <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. Problems compiling alps-2.0.3-src-with-boost (Andreas Sindermann) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 15 Feb 2012 13:19:02 +0100 (MET) > From: Andreas Sindermann <sinder(a)thp.Uni-Koeln.DE> > To: comp-phys-alps-users(a)lists.phys.ethz.ch > Subject: [ALPS-users] Problems compiling alps-2.0.3-src-with-boost > Message-ID: <201202151219.q1FCJ2p19838(a)sun80.thp.Uni-Koeln.DE> > Content-Type: text/plain; charset=US-ASCII > > > Dear all, > > I'm just about to compile and to install alps-2.0.3 from source with > libboost included. Compilation succeeds without problems. But > unfortunately 'make test' results in 4 failed tests: > > Test project /scratch/sinder/alps-2.0.3-src-with-boost/alps > Start 1: ising_single > 1/99 Test #1: ising_single .....................***Failed 0.32 sec > Start 2: loop_single > 2/99 Test #2: loop_single ......................***Failed 0.19 sec > Start 3: exchange > 3/99 Test #3: exchange .........................***Failed 0.25 sec > Start 40: lattice_constant_test > 40/99 Test #40: lattice_constant_test ............***Failed 0.08 sec > > > As I don't have experience with the Cmake building system I'd like to > ask for help where I should start looking. > > Thanks in advance > Andreas > > -- > Dr. Andreas Sindermann fon: +49 (221) 470-4201 > Institut fuer Theoretische Physik fax: +49 (221) 470-5159 > Universitaet zu Koeln > Zuelpicher Str. 77 mailto:sinder@thp.Uni-Koeln.DE > D-50937 Koeln, Germany http://www.thp.uni-koeln.de/~sinder > > > End of Comp-phys-alps-users Digest, Vol 71, Issue 9 > *************************************************** > -- Tadeusz Wasiutyński

4 17

use PBC condition
by muzi352 07 Mar '12

07 Mar '12

Dear Matthias: I begin with the Hamiltonian for a chain with spins S=1/2 , use periodic boundary conditions ,the total number of sites is N=20. calculate correlation : Offdiagonal spin correlations has problem, why 0-0 and 0-1 offdiagonal spin correlations value is 0 ?? | Offdiagonal spin correlations[( 0 ) -- ( 0 )] | 0 | | Offdiagonal spin correlations[( 0 ) -- ( 1 )] | 0 | | Offdiagonal spin correlations[( 1 ) -- ( 0 )] | 0.2499999999999999 | | Offdiagonal spin correlations[( 1 ) -- ( 1 )] | 0.01195256946918583 | | Offdiagonal spin correlations[( 0 ) -- ( 1 )] | -0.1410917216909388 | | Offdiagonal spin correlations[( 0 ) -- ( 2 )] | -0.02578705926914329 | | Offdiagonal spin correlations[( 1 ) -- ( 1 )] | 0.05437436640436828 | | Offdiagonal spin correlations[( 1 ) -- ( 2 )] | 0.01374443422199295 | | Offdiagonal spin correlations[( 0 ) -- ( 2 )] | -0.04722700778251831 | | Offdiagonal spin correlations[( 0 ) -- ( 3 )] | -0.01673969792911871 | | Offdiagonal spin correlations[( 1 ) -- ( 2 )] | 0.06699300813943075 | | Offdiagonal spin correlations[( 1 ) -- ( 3 )] | 0.002194337470049455 | | Offdiagonal spin correlations[( 0 ) -- ( 3 )] | -0.03664352593820266 | | Offdiagonal spin correlations[( 0 ) -- ( 4 )] | 0.002203905249778245 | | Offdiagonal spin correlations[( 1 ) -- ( 3 )] | 0.0331275686776057 | | Offdiagonal spin correlations[( 1 ) -- ( 4 )] | -0.007827950318578362 | | Offdiagonal spin correlations[( 0 ) -- ( 4 )] | -0.04668994676799755 | | Offdiagonal spin correlations[( 0 ) -- ( 5 )] | 0.01338526650557481 | | Offdiagonal spin correlations[( 1 ) -- ( 4 )] | 0.05418339350735944 | | Offdiagonal spin correlations[( 1 ) -- ( 5 )] | -0.02572536103232437 | | Offdiagonal spin correlations[( 0 ) -- ( 5 )] | 0 | | Offdiagonal spin correlations[( 0 ) -- ( 6 )] | 0 |

2 1