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

Re: [ALPS-users] Alps 1.3.1 spinmc compile problem
by Synge Todo 31 May '08

31 May '08

Does anyone have an experience on Ubuntsu Linux? Synge From: "Miles Stoudenmire" <miles.stoudenmire(a)gmail.com> Date: Fri, 30 May 2008 08:55:53 -0700 > Hi Synge - I'm using Ubuntu Linux on a machine with an AMD processor. > I tried installing the lapack 3.0 package through apt-get and saying > --with-lapack-dir=/usr/lib --with-lapack=-llapack because the lapack > library files got installed in /usr/lib (I have no idea if those > arguments are sensible, they were just a guess). When that didn't > work, I tried installing the ACML library from AMD and using a similar > approach but the configure script still said that it couldn't find any > LAPACK support. > > On Fri, May 30, 2008 at 1:08 AM, Synge Todo <wistaria(a)comp-phys.org> wrote: > > Dear Miles Stoudenmire, > > > > From: "Miles Stoudenmire" <miles(a)physics.ucsb.edu> > > Date: Thu, 29 May 2008 22:38:58 -0700 > >> It looks like it's referring to some funny lapack function from what I > >> can tell, so I think I should mention that I did not build the > >> libraries with lapack support (though the configure script claimed > >> that I'd be ok to compile classical monte carlo). > > > > Thank you for your report. It seems that the classical Monte Carlo > > also uses a SVD routine from lapack. So, this is a bug of configure > > script. I will fixed it in the next version. > > > >> If it does turn out > >> to be because I need to add lapack support, I will have to ask some > >> questions about how to do that. > > > > What kind of platform are you using? > > > > Best, > > Synge > > > > > > > > -- > -=Miles Stoudenmire=- > miles(a)physics.ucsb.edu > miles.stoudenmire(a)gmail.com > http://www.crazytheory.com/

3 2

ALPS 1.3.1,dirloop_sse: symmetry breaking
by Beni Thielemann 30 May '08

30 May '08

Hi all, I'm trying to calculate the Neel temperature in an Heisenberg antiferromagnet and wonder how the symmetry breaking is considered in dirloop_sse: Do I have to define the direction of symmetry breaking by including a small (small enough to not affect the results) magnetic field? If so, is there a way to do simulations with transverse staggered magnetic fields in ALPS/diroop_sse? (the <HAMILTONIAN name="spin"> node in models.xml looks like it only allows for a uniform transverse field Gamma) How would one specify this field in case that the unit cell is not primitive? What does dirloop_sse calculate below the Neel temperature in case no direction of symmetry breaking is specified (do one get an average over many different directions?) ? I'd be glad about any advice, Beni Thielemann -- Benedikt Thielemann, Ph. D. Student Laboratory for Neutron Scattering ETH Zurich and Paul Scherrer Institute 5232 Villigen, Switzerland Office: +41 56 310 3698 Mobile: +41 78 611 8074 Fax: +41 56 310 2939

1 0

Alps 1.3.1 spinmc compile problem
by Miles Stoudenmire 30 May '08

30 May '08

I'm trying to build the latest version of ALPS, and I'm primarily interested in the classical monte carlo app. I built the new version of the libraries ok, but then when I tried to compile the spinmc app, I got the error below. It looks like it's referring to some funny lapack function from what I can tell, so I think I should mention that I did not build the libraries with lapack support (though the configure script claimed that I'd be ok to compile classical monte carlo). If it does turn out to be because I need to add lapack support, I will have to ask some questions about how to do that. make[1]: Entering directory `/miles/alps-applications-1.3.1/mc/spins' /bin/bash ../../libtool --tag=CXX --mode=link g++ -pthread -ftemplate-depth-150 -O3 -DBOOST_DISABLE_ASSERTS -fexceptions -w -L/miles/ALPS/lib -o spinmc spinmc.o abstract_fitter.o base_incr_fitter.o dummy_fitter.o est_grad_fitter.o fitting_scheduler.o spinmc_factory.o -lalps -lcomm-sgl -lpthread -lm g++ -pthread -ftemplate-depth-150 -O3 -DBOOST_DISABLE_ASSERTS -fexceptions -w -o spinmc spinmc.o abstract_fitter.o base_incr_fitter.o dummy_fitter.o est_grad_fitter.o fitting_scheduler.o spinmc_factory.o -L/miles/ALPS/lib /miles/ALPS/lib/libalps.so /miles/ALPS/lib/libcomm-sgl.so -lpthread -lm -Wl,--rpath -Wl,/miles/ALPS/lib -Wl,--rpath -Wl,/miles/ALPS/lib est_grad_fitter.o: In function `double solve_llsp<double, boost::numeric::ublas::basic_column_major<unsigned int, int>, boost::numeric::ublas::unbounded_array<double, std::allocator<double> > >(boost::numeric::ublas::matrix<double, boost::numeric::ublas::basic_column_major<unsigned int, int>, boost::numeric::ublas::unbounded_array<double, std::allocator<double> > > const&, boost::numeric::ublas::vector<double, boost::numeric::ublas::unbounded_array<double, std::allocator<double> > > const&, boost::numeric::ublas::vector<double, boost::numeric::ublas::unbounded_array<double, std::allocator<double> > >&, double)': est_grad_fitter.C:(.text._Z10solve_llspIdN5boost7numeric5ublas18basic_column_majorIjiEENS2_15unbounded_arrayIdSaIdEEEET_RKNS2_6matrixIS8_T0_T1_EERKNS2_6vectorIS8_NS5_IS8_SaIS8_EEEEERSI_d[double solve_llsp<double, boost::numeric::ublas::basic_column_major<unsigned int, int>, boost::numeric::ublas::unbounded_array<double, std::allocator<double> > >(boost::numeric::ublas::matrix<double, boost::numeric::ublas::basic_column_major<unsigned int, int>, boost::numeric::ublas::unbounded_array<double, std::allocator<double> > > const&, boost::numeric::ublas::vector<double, boost::numeric::ublas::unbounded_array<double, std::allocator<double> > > const&, boost::numeric::ublas::vector<double, boost::numeric::ublas::unbounded_array<double, std::allocator<double> > >&, double)]+0x14a): undefined reference to `dgesvd_' est_grad_fitter.C:(.text._Z10solve_llspIdN5boost7numeric5ublas18basic_column_majorIjiEENS2_15unbounded_arrayIdSaIdEEEET_RKNS2_6matrixIS8_T0_T1_EERKNS2_6vectorIS8_NS5_IS8_SaIS8_EEEEERSI_d[double solve_llsp<double, boost::numeric::ublas::basic_column_major<unsigned int, int>, boost::numeric::ublas::unbounded_array<double, std::allocator<double> > >(boost::numeric::ublas::matrix<double, boost::numeric::ublas::basic_column_major<unsigned int, int>, boost::numeric::ublas::unbounded_array<double, std::allocator<double> > > const&, boost::numeric::ublas::vector<double, boost::numeric::ublas::unbounded_array<double, std::allocator<double> > > const&, boost::numeric::ublas::vector<double, boost::numeric::ublas::unbounded_array<double, std::allocator<double> > >&, double)]+0x5ac): undefined reference to `dgesvd_' collect2: ld returned 1 exit status make[1]: *** [spinmc] Error 1 make[1]: Leaving directory `/miles/alps-applications-1.3.1/mc/spins' make: *** [all-recursive] Error 1 -- -=Miles Stoudenmire=- miles(a)physics.ucsb.edu miles.stoudenmire(a)gmail.com http://www.crazytheory.com/

2 1

Monte Carlo - Extending Sweeps
by Miles Stoudenmire 29 May '08

29 May '08

Hi, I've been working with the ALPS classical monte carlo program, and according to the ALPS website/wiki, it should be possible to take a finished simulation, edit number of SWEEPS in the simulation_name.task#.out.xml files, set the status of the task from "finished" to "running" in the simulation_name.out.xml file and then do "spinmc simulation_name.out.xml". However, the program immediately throws an exception saying "Cannot change parameter SWEEPS". I determined that the exception is being thrown from line 300 of worker.C in alps-1.3.1/src/alps/scheduler/. Is this a feature that just never got implemented or am I going about the procedure incorrectly? Also, could ALPS (or does ALPS already) include a feature that allows it to run in a mode where it automatically advances to the next task if the current task has converged in all measurements, even if the full range of sweeps has not been performed? I think this would be an incredibly userful feature! -- -=Miles Stoudenmire=- miles(a)physics.ucsb.edu miles.stoudenmire(a)gmail.com http://www.crazytheory.com/

2 4

DMRG Truncation Error
by Niall Moran 26 May '08

26 May '08

Hi, I am using the ALPS DMRG code to look at a 1D system. I am comparing the results for small system sizes with the results from exact diagonalisation to verify them. The values match to 5 - 6 decimal places in most cases. However I always get a truncation error or 1. Is this normal behavior? I am using ALPS 1.3.1. Many Thanks, Niall.

1 0

ALPS 1.3.1 released
by Matthias Troyer 21 May '08

21 May '08

Dear colleagues, In preparation for the Sherbrooke summer school we releases ALPS 1.3.1 including the following bug fixes: * the hexagonal lattice is now correctly defined a a 3D lattice * the default random number generator was switched to a Mersenne Twister * a bug in calculating total energies in a classical Heisenberg model with single ion anisotropy was fixed It can be downloaded from the ALPS server at http://alps.comp-phys.org/ Matthias Troyer

1 0

What's ALPS/ISUM?
by Zhenhui Liu 15 May '08

15 May '08

hi, does anyone knows the ALPS/ISUM code? best regards

2 1

Can not download files from "User Workshops"To whom it may concern:
by dulei 15 May '08

15 May '08

To whom it may concern: Almost every file can not be downloaded form http://www.comp-phys.org/lugano04/Talks.html except Wang-Landau sampling(QWL.pdf) and Path Integrals and Quantum Monte Carlo (qmc.pdf). Another files are "Not Found. Requested URL /lugano04/Talks/IntroToCpp.pdf was not found on this server. Additionally, a 404 Not Found error was encountered while trying to use an ErrorDocument to handle the request." Please tell me where I can download these files. Thank you very much. With best regards. Yours Sincerely. 2008-05-15 dulei

2 1