Comp-phys-alps-users

comp-phys-alps-users@lists.phys.ethz.ch

1244 discussions

dirloop_sse density vs green's function measurement
by Francisco Cordobés 30 Jan '13

30 Jan '13

Dear all, I have been doing some tests for the Bose Hubbard model in 1D with dirloop_sse and I find that the the output density and the value of the Green's function at (0,0) do not coincide, as I would expect. I understand that when measuring Green's functions for the Bose-Hubbard model the output is that of <b^\dagger_0 b_j> with j going over all lattice points. b^\dagger_i is a creation operator and b_j is a destruction at sites i and j respectively. If that is correct, then, the following should be correct too and I should get <b^\dagger_0 b_0>=<n_0>=\rho. Where I denote the density by \rho. When I do the calculation, using the following input: <?xml version="1.0" encoding="UTF-8"?> <?xml-stylesheet type="text/xsl" href="ALPS.xsl"?> <SIMULATION xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://xml.comp-phys.org/2003/8/QMCXML.xsd"> <PARAMETERS> <PARAMETER name="MEASURE[Green Function]">True</PARAMETER> <PARAMETER name="MEASURE[Correlations]">True</PARAMETER> <PARAMETER name="L">64</PARAMETER> <PARAMETER name="LATTICE">chain lattice</PARAMETER> <PARAMETER name="SEED">189185290</PARAMETER> <PARAMETER name="U">10.0</PARAMETER> <PARAMETER name="T">0.5</PARAMETER> <PARAMETER name="MODEL">boson Hubbard</PARAMETER> <PARAMETER name="Nmax">3</PARAMETER> <PARAMETER name="THERMALIZATION">1000</PARAMETER> <PARAMETER name="SWEEPS">50000</PARAMETER> <PARAMETER name="NONLOCAL">0</PARAMETER> <PARAMETER name="mu">0.4</PARAMETER> <PARAMETER name="t">0.5</PARAMETER> </PARAMETERS> </SIMULATION> The output I get is: <SCALAR_AVERAGE name="Density"> <COUNT>55257</COUNT> <MEAN method="simple">0.6790808</MEAN> <ERROR converged="yes" method="binning">0.000306</ERROR> <VARIANCE method="simple">0.00258</VARIANCE> <AUTOCORR method="binning">0.505</AUTOCORR> <VECTOR_AVERAGE name="Green's Function" nvalues="64"> <SCALAR_AVERAGE indexvalue="( 0 ) -- ( 0 )"> <COUNT>55257</COUNT> <MEAN method="simple">1.178037</MEAN> <ERROR converged="yes" method="binning">0.000505</ERROR> <VARIANCE method="simple">0.0117</VARIANCE> <AUTOCORR method="binning">0.103</AUTOCORR> Which is very different. Am I interpreting the output correctly or am I just not running it long enough? Also, what is the meaning of the last three fields in the xml output field? n,n^2 and n^4? Because my initial guess would be to interpret n as the number of particles but it is not consistent with the density and the lattice size. Thank you very much, Francisco Cordobes .

2 2

hey
by Majid Abedi 29 Jan '13

29 Jan '13

have a look at this http://bit.ly/VQiD6S Majid

1 0

alps build error?
by Mateusz Łącki 29 Jan '13

29 Jan '13

Dear All, I am compiling alps and I am stuct at 64%. The error message is at the very bottom. My cmake line is : cmake28 -DCMAKE_INSTALL_PREFIX=~/opt/alps -D Boost_ROOT_DIR:PATH=~/opt/alps-2.1.1-r6176-src-with-boost/boost/ -DCMAKE_CXX_COMPILER=/opt/intel/composer_xe_2011_sp1.7.256/bin/intel64/icpc -DCMAKE_C_COMPILER=/opt/intel/composer_xe_2011_sp1.7.256/bin/intel64/icc -D LAPACK_64_BIT=ON -D HDF5_LIBRARIES=/opt/hdf5-1.8.8-intel/lib/libhdf5.so -D HDF5_INCLUDE_DIR=/opt/hdf5-1.8.8-intel/include/ -DALPS_BUILD_FORTRAN=ON ~/opt/alps-2.1.1-r6176-src-with-boost/alps/ compliler lines: -- The C compiler identification is Intel 12.1.0.20111011 -- The CXX compiler identification is Intel 12.1.0.20111011 -- The Fortran compiler identification is Intel What could be wrong ? I am using version alps-2.1.1-r6176-src-with-boost (straing from the web page, with FindLapack taken from the alps-2.1.1-r6670-src-with-boost). The version alps-2.1.1-r6670-src-with-boost fails at exactly the same point with exactly the same error message. Regards, Mateusz Lacki [ 64%] Building CXX object example/ietl/CMakeFiles/jacobidavidson2.dir/jacobidavidson2.cpp.o /home/lkb/lacki/opt/alps-2.1.1-r6176-src-with-boost/alps/src/boost/numeric/bindings/lapack/computational/getrf.hpp(197): error: no instance of overloaded function "boost::numeric::bindings::lapack::detail::getrf" matches the argument list argument types are: (order, ptrdiff_t, ptrdiff_t, double *, ptrdiff_t, int *) return detail::getrf( order(), bindings::size_row(a), ^ detected during: instantiation of "ptrdiff_t={long} boost::numeric::bindings::lapack::getrf_impl<Value>::invoke(MatrixA &, VectorIPIV &) [with Value=double, MatrixA=boost::numeric::ublas::matrix<double, boost::numeric::ublas::column_major, boost::numeric::ublas::unbounded_array<double, std::allocator<double>>>, VectorIPIV=std::vector<int, std::allocator<int>>]" at line 220 instantiation of "ptrdiff_t={long} boost::numeric::bindings::lapack::getrf(MatrixA &, VectorIPIV &) [with MatrixA=boost::numeric::ublas::matrix<double, boost::numeric::ublas::column_major, boost::numeric::ublas::unbounded_array<double, std::allocator<double>>>, VectorIPIV=std::vector<int, std::allocator<int>>]" at line 285 of "/home/lkb/lacki/opt/alps-2.1.1-r6176-src-with-boost/alps/src/ietl/jd.h" instantiation of "void ietl::solver::left_prec_solver<SOLV, MATRIX, VS, PREC>::operator()(ietl::real_type<ietl::vectorspace_traits<VS>::scalar_type>::type, ietl::vectorspace_traits<VS>::vector_type &, ietl::vectorspace_traits<VS>::vector_type &, IT &) [with SOLV=ietl::gmres_wrapper, MATRIX=matrix_t, VS=ietl::vectorspace<vector_t>, PREC=jacobi_prec<matrix_t, double>, IT=ietl::jd_iteration<double>]" at line 611 of "/home/lkb/lacki/opt/alps-2.1.1-r6176-src-with-boost/alps/src/ietl/jd.h" instantiation of "void ietl::jd<MATRIX, VS>::jdqr(SOLVER &, IT &, GEN &, ietl::vectorspace_traits<VS>::size_type, bool) [with MATRIX=matrix_t, VS=ietl::vectorspace<vector_t>, SOLVER=ietl::solver::left_prec_solver<ietl::gmres_wrapper, matrix_t, ietl::vectorspace<vector_t>, jacobi_prec<matrix_t, double>>, IT=ietl::jd_iteration<double>, GEN=gen_t]" at line 384 of "/home/lkb/lacki/opt/alps-2.1.1-r6176-src-with-boost/alps/src/ietl/jd.h" instantiation of "void ietl::jd<MATRIX, VS>::eigensystem(IT &, GEN &, ietl::vectorspace_traits<VS>::size_type, PREC &, SOLV, bool) [with MATRIX=matrix_t, VS=ietl::vectorspace<vector_t>, SOLV=ietl::gmres_wrapper, IT=ietl::jd_iteration<double>, GEN=gen_t, PREC=jacobi_prec<matrix_t, double>]" at line 165 of "/home/lkb/lacki/opt/alps-2.1.1-r6176-src-with-boost/alps/example/ietl/jacobidavidson2.cpp" /home/lkb/lacki/opt/alps-2.1.1-r6176-src-with-boost/alps/src/boost/numeric/bindings/lapack/computational/getrs.hpp(222): error: no instance of overloaded function "boost::numeric::bindings::lapack::detail::getrs" matches the argument list argument types are: (order, trans, ptrdiff_t, ptrdiff_t, const double *, ptrdiff_t, const int *, double *, ptrdiff_t) return detail::getrs( order(), trans(), bindings::size_column_op(a, ^ detected during: instantiation of "ptrdiff_t={long} boost::numeric::bindings::lapack::getrs_impl<Value>::invoke(const MatrixA &, const VectorIPIV &, MatrixB &) [with Value=double, MatrixA=boost::numeric::ublas::matrix<double, boost::numeric::ublas::column_major, boost::numeric::ublas::unbounded_array<double, std::allocator<double>>>, VectorIPIV=std::vector<int, std::allocator<int>>, MatrixB=boost::numeric::ublas::vector<double, boost::numeric::ublas::unbounded_array<double, std::allocator<double>>>]" at line 247 instantiation of "ptrdiff_t={long} boost::numeric::bindings::lapack::getrs(const MatrixA &, const VectorIPIV &, MatrixB &) [with MatrixA=boost::numeric::ublas::matrix<double, boost::numeric::ublas::column_major, boost::numeric::ublas::unbounded_array<double, std::allocator<double>>>, VectorIPIV=std::vector<int, std::allocator<int>>, MatrixB=boost::numeric::ublas::vector<double, boost::numeric::ublas::unbounded_array<double, std::allocator<double>>>]" at line 294 of "/home/lkb/lacki/opt/alps-2.1.1-r6176-src-with-boost/alps/src/ietl/jd.h" instantiation of "void ietl::solver::left_prec_solver<SOLV, MATRIX, VS, PREC>::operator()(ietl::real_type<ietl::vectorspace_traits<VS>::scalar_type>::type, ietl::vectorspace_traits<VS>::vector_type &, ietl::vectorspace_traits<VS>::vector_type &, IT &) [with SOLV=ietl::gmres_wrapper, MATRIX=matrix_t, VS=ietl::vectorspace<vector_t>, PREC=jacobi_prec<matrix_t, double>, IT=ietl::jd_iteration<double>]" at line 611 of "/home/lkb/lacki/opt/alps-2.1.1-r6176-src-with-boost/alps/src/ietl/jd.h" instantiation of "void ietl::jd<MATRIX, VS>::jdqr(SOLVER &, IT &, GEN &, ietl::vectorspace_traits<VS>::size_type, bool) [with MATRIX=matrix_t, VS=ietl::vectorspace<vector_t>, SOLVER=ietl::solver::left_prec_solver<ietl::gmres_wrapper, matrix_t, ietl::vectorspace<vector_t>, jacobi_prec<matrix_t, double>>, IT=ietl::jd_iteration<double>, GEN=gen_t]" at line 384 of "/home/lkb/lacki/opt/alps-2.1.1-r6176-src-with-boost/alps/src/ietl/jd.h" instantiation of "void ietl::jd<MATRIX, VS>::eigensystem(IT &, GEN &, ietl::vectorspace_traits<VS>::size_type, PREC &, SOLV, bool) [with MATRIX=matrix_t, VS=ietl::vectorspace<vector_t>, SOLV=ietl::gmres_wrapper, IT=ietl::jd_iteration<double>, GEN=gen_t, PREC=jacobi_prec<matrix_t, double>]" at line 165 of "/home/lkb/lacki/opt/alps-2.1.1-r6176-src-with-boost/alps/example/ietl/jacobidavidson2.cpp" compilation aborted for /home/lkb/lacki/opt/alps-2.1.1-r6176-src-with-boost/alps/example/ietl/jacobidavidson2.cpp (code 2) make[2]: *** [example/ietl/CMakeFiles/jacobidavidson2.dir/jacobidavidson2.cpp.o] Erreur 2 make[1]: *** [example/ietl/CMakeFiles/jacobidavidson2.dir/all] Erreur 2 make: *** [all] Erreur 2

3 3

FindLapack.cmake
by Mateusz Łącki 25 Jan '13

25 Jan '13

Dear All I have found the following code in FindLapack.cmake trying to do installation of new alps: IF($ENV{MKL} MATCHES "mkl") MESSAGE(STATUS "Using intel/mkl library: $ENV{MKL}") #look for the path where the mkl is installed 43: STRING(REGEX MATCHALL "[-][L]([^ ;])+" MKL_PATH_WITH_PREFIX "$ENV{MKL}") STRING(REGEX REPLACE "[-][L]" "" MKL_PATHS ${MKL_PATH_WITH_PREFIX}) ENDIF($ENV{MKL} MATCHES "mkl") If I assume correctly this code assumes that there is an environment variable dalled MKL which contains something like "-L/path/to/mkl" on my platform it is set by default to "/path/to/mkl" If I understand it correctly at the line 43, the $MKL gets stripped of "-L" part and the ${MKL_PATH_WITH_PREFIX} then will contain just "/path/to/mkl". I think that perhaps some sort of conditional clause could be added so that both $MKL = /path/to/mkl and $MKL = -L/path/to/mkl worked. Regards, Mateusz Lacki

2 3

alps installation: lapack vs intel/mkl
by Mateusz Łącki 25 Jan '13

25 Jan '13

Dear All, I would like to ask about the lapack library used by alps: I have a following output: -- Using intel/mkl library: -L/opt/intel/composer_xe_2011_sp1.7.256/mkl MKL_PATH_WITH_PREFIX -L/opt/intel/composer_xe_2011_sp1.7.256/mkl -- Looking for intel/mkl library in /opt/intel/composer_xe_2011_sp1/mkl//lib/intel64;/opt/intel/composer_xe_2011_sp1/mkl//lib/emt;/opt/intel/composer_xe_2011_sp1/mkl//lib/;/opt/intel/composer_xe_2011_sp1/mkl//lib;/opt/intel/composer_xe_2011_sp1.7.256/mkl 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 not found. Please specify library location. (i have uncommented debug messages) so my mkl is located in: /opt/intel/composer_xe_2011_sp1/mkl//lib/intel64 which is btw the default localtion of mkl library in newer versions of the compiler in the FindLapack.cmake the lines 61-65 are: STRING(REGEX MATCH "[0-9]+\\.[0-9]+\\.[0-9]+" MKL_VERSION "$ENV{MKL_HOME}") SET(LINK_LAPACK_DEFAULT 1) IF(${MKL_VERSION} MATCHES "10\\.0\\.[0-2]") SET(LINK_LAPACK_DEFAULT 0) ENDIF(${MKL_VERSION} MATCHES "10\\.0\\.[0-2]") I think (but I have never ever used cmake in my project) that this lines assume that mkl s installed in a directory that that contains specific version: /opt/intel/Compiler/11.5.6/mkl (then MKL_Version becomes 11.5.6) which is often not the case. I will try to fix this cmake file as well as I can... Regards, Mateusz Lacki

1 0

Version of Boost in the ALPS+Boost download
by evanloon＠physnet.uni-hamburg.de 24 Jan '13

24 Jan '13

Dear all, It seems that the alps src-with-boost downloads still contain an old (1.40) version of the boost src, which will not allow to compile ALPS. Best regards, Erik

6 9

DMRG simulations of Fermi Hubbard models
by Giuliano Orso 24 Jan '13

24 Jan '13

Dear ALPS users, I want to run some DMRG simulations on 1D Fermi Hubbard models. I tried to work in the gran canonical ensemble, giving the chemical potential mu instead of the density and without specifying that Nup and Ndown are conserved quantities. The code is running but it produces strange results, in particular the calculated energy is varying a lot within the same sweep and the error oscillates from 10-15 to 10-2. Any ideas why ? Puzzled by that, I tried to do the same for the Bose Hubbard model, but there everything works perfectly. Can you help me on that ? My best regards, Giuliano

1 0

TEBD and custom lattices
by Deepak Iyer 21 Jan '13

21 Jan '13

Hi I'm trying to define a custom lattice to use with the tebd code. I've tested the lattice with sparsediag and it works as I want it to. However, TEBD seems to completely ignore any information about the lattice. It seems to look at only the model, and I guess assumes an open chain of length L with only nearest neighbor hopping. Is there a way to get around this? Thanks! Deepak

1 0

site and bond-operators in the models.xml file
by Deepak Iyer 20 Jan '13

20 Jan '13

Hi I'm confused as to how the system knows what i and j should be (like nearest neighbor, next nearest, etc) for a particular model, since the bond-operators just go from i to j. Also, when there are wild cards in parameters, e.g., mu# in the boson Hubbard Hamiltonian, where are the different types defined? Do I need to define a new lattice type to allow for spatially varying mu, e.g.? Thanks Deepak

2 2

Quenching 'L' in TEBD code
by Deepak Iyer 19 Jan '13

19 Jan '13

Hi I'm trying to quench 'L' in the TEBD code in the Bose-Hubbard model. I start with L=2 and quench to L=10. I'm interested in looking at what happens to the local occupation numbers. But the code only outputs 'Local Number' for 2 sites. I'm using an L=2 initial state to force the system into a particular initial state. I see that this is a bit ambiguous since its not clear what 'L' is midway in the quench, but is there any way I can get occupation for all 10 sites after the completion of the quench? I first quickly quench it and then hold the quench at a fixed U and L for sometime. Thanks! Deepak

2 4

Jump to page:

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

2007

2006

2005

2004

Comp-phys-alps-users