December 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

the entanglement Hamiltonian
by hamid mosadegh 05 Jan '15

05 Jan '15

Dear ALPS user I want to calculate the left reduced density matrix \rho_L=Tr_R |\Phi><\Phi| using DMRG and then evaluate the entanglement spectrum consists of the energies of the entanglement Hamiltonian H_E=-Ln(\rho_L) But we don't have direct access to the many body wave function in the DMRG calculations using ALPS. Is there any way to evaluate these parameters in the DMRG calculations using ALPS? Thanks -- H. Mosadeq Shahr-e-Kord University (SKU), Shahr-e-Kord , Iran

2 1

problems regarding the MPI runs in the ALPS
by D´Souza, Sunil Wilfred 12 Dec '14

12 Dec '14

Dear Alps users, I have installed ALPS2.2 on my local machine containing 12 processing units. I have tried running a simple tutorial given in the alps examples files. when i am using the following command "pyalps.runApplication('spinmc',input_file,Tmin=5)" , everything works fine and the ALPS program runs (Serial mode) and gives the proper output as expected. On the other hand when i use the command "pyalps.runApplication('spinmc',input_file,Tmin=5,MPI=7)" , i am getting a lot of error messages related to segmentation fault. The error messages are as follows: mpirun -np 7 spinmc --mpi parm2a.in.xml --Tmin 5 [LIN21:04682] *** Process received signal *** [LIN21:04682] Signal: Segmentation fault (11) [LIN21:04682] Signal code: Address not mapped (1) [LIN21:04682] Failing at address: (nil) [LIN21:04682] [ 0] /lib64/libpthread.so.0(+0xf9f0) [0x7f356ffcf9f0] [LIN21:04682] [ 1] /opt/openmpi/1.4.5/gcc/lib64/openmpi/mca_pml_v.so(+0x1b62) [0x7f3569f92b62] [LIN21:04682] [ 2] /opt/openmpi/1.4.5/gcc/lib64/libopen-pal.so.0(mca_base_components_close+0x76) [0x7f356f2616c6] [LIN21:04682] [ 3] /opt/openmpi/1.4.5/gcc/lib64/libmpi.so.0(mca_pml_base_select+0x301) [0x7f356f765d11] [LIN21:04682] [ 4] /opt/openmpi/1.4.5/gcc/lib64/libmpi.so.0(+0x39acc) [0x7f356f725acc] [LIN21:04682] [ 5] /opt/openmpi/1.4.5/gcc/lib64/libmpi.so.0(MPI_Init+0x16b) [0x7f356f745a8b] [LIN21:04682] [ 6] /home/sunil/ALPS2.2/lib/libalps.so.2(_ZN4alps9comm_initERiRPPcb+0x58) [0x7f3570f60098] [LIN21:04682] [ 7] /home/sunil/ALPS2.2/lib/libalps.so.2(_ZN4alps9scheduler5startEiPPcRKNS0_7FactoryE+0x37) [0x7f35711a08e7] [LIN21:04682] [ 8] spinmc(main+0x15) [0x449a95] [LIN21:04682] [ 9] /lib64/libc.so.6(__libc_start_main+0xf5) [0x7f356d184be5] [LIN21:04682] [10] spinmc() [0x449de5] [LIN21:04682] *** End of error message *** -------------------------------------------------------------------------- mpirun noticed that process rank 1 with PID 4682 on node LIN21 exited on signal 11 (Segmentation fault). -------------------------------------------------------------------------- I am very kindly requesting the ALPs users to help me out in solving the above mentioned issue. I am looking forward for your suggestions to overcome this problem. with best regards Sunil

2 1

Exact Diagonalization Fermion Hubbard
by Sriluckshmy Viswanathan 12 Dec '14

12 Dec '14

Dear all, I am trying to solve the fermion hubbard model on a square lattice with a small cluster size of 2x2 using exact diagonalization. I want to look at the energy levels at half filling. I used to N_total parameter to fix the filling, but the energy levels results dont change if I change this parameter. Could you please help me with this problem? I have attached the param file that I used. I am using fulldiag of ALPS. Regards Sriluckshmy

3 5

Re: [ALPS-users] Comp-phys-alps-users Digest, Vol 105, Issue 1
by Mateusz Łącki 01 Dec '14

01 Dec '14

ok, to weź poproś - chciałem Ci skopiować oprócz skryptów wyniki minimalizacji energii dla tych danych dla których om chciał wyciągnąć „Sz“. Nie miałem póki co czasu tego zrobić (bo Kuba chciał inne rzeczy „na wczoraj“, a także tutaj mam sporo czasochłonnych obowiązków). > On 01 Dec 2014, at 12:00, 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. Exact Diagonalization Fermion Hubbard (Sriluckshmy Viswanathan) > 2. Re: Exact Diagonalization Fermion Hubbard (Michele Dolfi) > 3. Re: Exact Diagonalization Fermion Hubbard > (Sriluckshmy Viswanathan) > 4. Re: bug in mps_optim (Michele Dolfi) > 5. Re: Exact Diagonalization Fermion Hubbard (Matthias Troyer) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sun, 30 Nov 2014 16:54:43 +0530 > From: Sriluckshmy Viswanathan <sriluckus(a)gmail.com> > To: comp-phys-alps-users <comp-phys-alps-users(a)lists.phys.ethz.ch> > Subject: [ALPS-users] Exact Diagonalization Fermion Hubbard > Message-ID: > <CAJTRBgcmbU1Fz=zr+fPxRm9id6JAD66RPBYgn_JEouVkbbjuXw(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Dear all, > > I am trying to solve the fermion hubbard model on a square lattice with a > small cluster size of 2x2 using exact diagonalization. I want to look at > the energy levels at half filling. I used to N_total parameter to fix the > filling, but the energy levels results dont change if I change this > parameter. Could you please help me with this problem? I have attached the > param file that I used. I am using fulldiag of ALPS. > > Regards > Sriluckshmy >

1 0