November 2011 - Comp-phys-alps-users

Structure factor
by "Michał Maik" 22 Nov '11

22 Nov '11

Hello, I had another question regarding the structure factor of the 2D triangular lattice. I am currently using the SSE and worms code (interchangably) in order to study spins or hardcore bosons on a triangular lattice in finite temperature. Using these codes I am able to get the spin or density structure factor but only for the Q=(0, 0) case. I am interested in the Q=(4*Pi/3, 0) case in order to search for supersolid phases. Is there a way to specify the different wave vectors in these programs in order to get the correct structure factors? Thank you for the help. Best Regards Michal

1 0

cmake error
by Mateusz Łącki 22 Nov '11

22 Nov '11

Dear All I wish to report the following error that occured while trying to reconfigure the cmake scripts: CMake Error at config/FindLapack.cmake:43 (STRING): string sub-command REGEX, mode REPLACE needs at least 6 arguments total to command. Call Stack (most recent call first): CMakeLists.txt:315 (find_package) this shows in ccmake after pressing "c". Therefore option to "generate" does not show. and this:-- MPI compiler was /opt/openmpi-1.5-intel/bin/mpic++ -- Using intel/mkl library: /opt/intel/composerxe-2011.1.107/mkl CMake Error at config/FindLapack.cmake:43 (STRING): string sub-command REGEX, mode REPLACE needs at least 6 arguments total to command. Call Stack (most recent call first): CMakeLists.txt:315 (find_package) during make Regards, Mateusz Lacki

3 7

Re: [ALPS-users] BONDOPERATOR measurement question
by Sebastian Huber 22 Nov '11

22 Nov '11

Attached the input files. If I remeber correctly, with alps 2.0.0 I did not encounter this problem. Sebastian -- Sebastian Huber Department of Condensed Matter Physics The Weizmann Institute of Science 76100 Rehovot Web: http://www.weizmann.ac.il/~huber Tel: +972 8 934 4258 On 22.11.2011, at 09:36, Matthias Troyer wrote: > Can you please send the full input files as attachment and I will take a look? > > Matthias > > On Nov 22, 2011, at 12:17 AM, Sebastian Huber wrote: > >> Didn't help. I played a bit with operator names. But the fact that MEASURE_LOCAL does perform something made me believe that it is not the cure. >> >> Thanks for the quick answer, Sebastian >> -- >> Sebastian Huber >> Department of Condensed Matter Physics >> The Weizmann Institute of Science >> 76100 Rehovot >> >> Web: http://www.weizmann.ac.il/~huber >> Tel: +972 8 934 4258 >> >> On 22.11.2011, at 01:11, Matthias Troyer wrote: >> >>> Try not to use a minus (-) in the operator name, but instead an underscore or similar >>> >>> >>> On 22 Nov 2011, at 00:09, Sebastian Huber wrote: >>> >>>> Dear all, >>>> >>>> I have an inhomogeneous lattice problem and I try to measure a "local" bondoperator, e.g., the current over one bond only. I get the following error from the sparsediag application >>>> >>>> Did not parse to end of string '(one-bond)/4(i)' >>>> >>>> when I try to do a "MEASURE_AVERAGE". I get no error if I use "MEASURE_LOCAL". I am a bit puzzled that "MEASURE_LOCAL" does something all together as my operator contains more than one site, which is in violation to a condition on local measurements on the alps webpage. I paste a minimal (silly) set of configuration files which produces this error below. Help would be greatly appreciated. >>>> >>>> Sebastian >>>> >>>> lanczosconf: >>>> >>>> MODEL="simple"; >>>> LATTICE="2x2" >>>> LATTICE_LIBRARY="graph.xml" >>>> MODEL_LIBRARY="model.xml" >>>> COMPLEX=true; >>>> t=-1; >>>> Nup_total=1; >>>> Ndown_total=1; >>>> CONSERVED_QUANTUMNUMBERS="N,Sz"; >>>> NUMBER_EIGENVALUES=2; >>>> TRANSLATION_SYMMETRY=false; >>>> MEASURE_AVERAGE[one-bond]=one-bond; >>>> {U=1;} >>>> >>>> >>>> graph.xml: >>>> >>>> <LATTICES> >>>> >>>> <GRAPH name="2x2" dimension="2" vertices="4" edges="9"> >>>> <VERTEX id="1" type="0"><COORDINATE>0 0</COORDINATE></VERTEX> >>>> <VERTEX id="2" type="0"><COORDINATE>1 0</COORDINATE></VERTEX> >>>> <VERTEX id="3" type="0"><COORDINATE>0 1</COORDINATE></VERTEX> >>>> <VERTEX id="4" type="0"><COORDINATE>1 1</COORDINATE></VERTEX> >>>> <EDGE source="1" target="2" id="1" type="1" vector="1 0"/> >>>> <EDGE source="2" target="1" id="2" type="2" vector="1 0"/> >>>> <EDGE source="3" target="4" id="3" type="3" vector="1 0"/> >>>> <EDGE source="4" target="3" id="4" type="4" vector="1 0"/> >>>> <EDGE source="1" target="3" id="5" type="5" vector="0 1"/> >>>> <EDGE source="2" target="4" id="6" type="6" vector="0 1"/> >>>> <EDGE source="3" target="1" id="7" type="7" vector="0 1"/> >>>> <EDGE source="4" target="2" id="8" type="8" vector="0 1"/> >>>> <EDGE source="1" target="4" id="9" type="9" vector="1 1"/> >>>> </GRAPH> >>>> >>>> </LATTICES> >>>> >>>> >>>> model.xml >>>> >>>> <MODELS> >>>> <SITEBASIS name="fermion"> >>>> <QUANTUMNUMBER name="Nup" min="0" max="1" type="fermionic"/> >>>> <QUANTUMNUMBER name="Ndown" min="0" max="1" type="fermionic"/> >>>> <OPERATOR name="Splus" matrixelement="1/2"> >>>> <CHANGE quantumnumber="Nup" change="1"/> >>>> <CHANGE quantumnumber="Ndown" change="-1"/> >>>> </OPERATOR> >>>> <OPERATOR name="Sminus" matrixelement="1/2"> >>>> <CHANGE quantumnumber="Nup" change="-1"/> >>>> <CHANGE quantumnumber="Ndown" change="+1"/> >>>> </OPERATOR> >>>> <OPERATOR name="Sz" matrixelement="(Nup-Ndown)/2"/> >>>> <OPERATOR name="Nup" matrixelement="Nup"/> >>>> <OPERATOR name="Ndown" matrixelement="Ndown"/> >>>> <OPERATOR name="c_down" matrixelement="1"> >>>> <CHANGE quantumnumber="Ndown" change="-1"/> >>>> </OPERATOR> >>>> <OPERATOR name="cdag_down" matrixelement="1"> >>>> <CHANGE quantumnumber="Ndown" change="1"/> >>>> </OPERATOR> >>>> <OPERATOR name="c_up" matrixelement="1"> >>>> <CHANGE quantumnumber="Nup" change="-1"/> >>>> </OPERATOR> >>>> <OPERATOR name="cdag_up" matrixelement="1"> >>>> <CHANGE quantumnumber="Nup" change="1"/> >>>> </OPERATOR> >>>> <OPERATOR name="n" matrixelement="Nup+Ndown"/> >>>> <OPERATOR name="n_up" matrixelement="Nup"/> >>>> <OPERATOR name="n_down" matrixelement="Ndown"/> >>>> </SITEBASIS> >>>> >>>> <BASIS name="fermion"> >>>> <SITEBASIS ref="fermion"/> >>>> <CONSTRAINT quantumnumber="Nup" value="Nup_total"/> >>>> <CONSTRAINT quantumnumber="Ndown" value="Ndown_total"/> >>>> </BASIS> >>>> >>>> <BONDOPERATOR name="one-bond" source="1" target="2"> >>>> I*(cdag_up(1)*c_up(2)-cdag_up(2)*c_up(1)) >>>> </BONDOPERATOR> >>>> >>>> <HAMILTONIAN name="simple"> >>>> <PARAMETER name="U" default="1"/> >>>> <BASIS ref="fermion"/> >>>> <SITETERM site="i"> >>>> U*n_up(i)*n_down(i) >>>> </SITETERM> >>>> </HAMILTONIAN> >>>> >>>> </MODELS> >>>> >>>> -- >>>> Sebastian Huber >>>> Department of Condensed Matter Physics >>>> The Weizmann Institute of Science >>>> 76100 Rehovot >>>> >>>> Web: http://www.weizmann.ac.il/~huber >>>> Tel: +972 8 934 4258 >>>> >>> >> >

1 0

Complex hopping intensity on some graph
by Volodymyr Derzhko 22 Nov '11

22 Nov '11

Hi, I introduced complex hopping parameter as follows (in the models.xml file): <BONDOPERATOR name="fermion_hop_l" source="x" target="y"> cdag_up(x)*c_up(y)+cdag_down(x)*c_down(y) </BONDOPERATOR> <BONDOPERATOR name="fermion_hop_r" source="x" target="y"> cdag_up(y)*c_up(x)+cdag_down(y)*c_down(x) </BONDOPERATOR> and later, as a bound term, <BONDTERM source="i" target="j"> <PARAMETER name="tl#" default="0"/> <PARAMETER name="tr#" default="0"/> -tl#*fermion_hop_l(i,j) -tr#*fermion_hop_r(i,j) </BONDTERM> In the parameters file, I put tr=tl*, so it seems to me, that such definition always lead to hermitian matrix, and do not depend on lattice\graph. Nevertheless, it does not work on a graph ('can not convert complex number into real one' error) but does work on predefined square lattice. What is the problem? Is the definition of the Hamiltonian is correct? Thank you in advance, Volodymyr Derzhko

2 5

Re: [ALPS-users] BONDOPERATOR measurement question
by Sebastian Huber 22 Nov '11

22 Nov '11

Didn't help. I played a bit with operator names. But the fact that MEASURE_LOCAL does perform something made me believe that it is not the cure. Thanks for the quick answer, Sebastian -- Sebastian Huber Department of Condensed Matter Physics The Weizmann Institute of Science 76100 Rehovot Web: http://www.weizmann.ac.il/~huber Tel: +972 8 934 4258 On 22.11.2011, at 01:11, Matthias Troyer wrote: > Try not to use a minus (-) in the operator name, but instead an underscore or similar > > > On 22 Nov 2011, at 00:09, Sebastian Huber wrote: > >> Dear all, >> >> I have an inhomogeneous lattice problem and I try to measure a "local" bondoperator, e.g., the current over one bond only. I get the following error from the sparsediag application >> >> Did not parse to end of string '(one-bond)/4(i)' >> >> when I try to do a "MEASURE_AVERAGE". I get no error if I use "MEASURE_LOCAL". I am a bit puzzled that "MEASURE_LOCAL" does something all together as my operator contains more than one site, which is in violation to a condition on local measurements on the alps webpage. I paste a minimal (silly) set of configuration files which produces this error below. Help would be greatly appreciated. >> >> Sebastian >> >> lanczosconf: >> >> MODEL="simple"; >> LATTICE="2x2" >> LATTICE_LIBRARY="graph.xml" >> MODEL_LIBRARY="model.xml" >> COMPLEX=true; >> t=-1; >> Nup_total=1; >> Ndown_total=1; >> CONSERVED_QUANTUMNUMBERS="N,Sz"; >> NUMBER_EIGENVALUES=2; >> TRANSLATION_SYMMETRY=false; >> MEASURE_AVERAGE[one-bond]=one-bond; >> {U=1;} >> >> >> graph.xml: >> >> <LATTICES> >> >> <GRAPH name="2x2" dimension="2" vertices="4" edges="9"> >> <VERTEX id="1" type="0"><COORDINATE>0 0</COORDINATE></VERTEX> >> <VERTEX id="2" type="0"><COORDINATE>1 0</COORDINATE></VERTEX> >> <VERTEX id="3" type="0"><COORDINATE>0 1</COORDINATE></VERTEX> >> <VERTEX id="4" type="0"><COORDINATE>1 1</COORDINATE></VERTEX> >> <EDGE source="1" target="2" id="1" type="1" vector="1 0"/> >> <EDGE source="2" target="1" id="2" type="2" vector="1 0"/> >> <EDGE source="3" target="4" id="3" type="3" vector="1 0"/> >> <EDGE source="4" target="3" id="4" type="4" vector="1 0"/> >> <EDGE source="1" target="3" id="5" type="5" vector="0 1"/> >> <EDGE source="2" target="4" id="6" type="6" vector="0 1"/> >> <EDGE source="3" target="1" id="7" type="7" vector="0 1"/> >> <EDGE source="4" target="2" id="8" type="8" vector="0 1"/> >> <EDGE source="1" target="4" id="9" type="9" vector="1 1"/> >> </GRAPH> >> >> </LATTICES> >> >> >> model.xml >> >> <MODELS> >> <SITEBASIS name="fermion"> >> <QUANTUMNUMBER name="Nup" min="0" max="1" type="fermionic"/> >> <QUANTUMNUMBER name="Ndown" min="0" max="1" type="fermionic"/> >> <OPERATOR name="Splus" matrixelement="1/2"> >> <CHANGE quantumnumber="Nup" change="1"/> >> <CHANGE quantumnumber="Ndown" change="-1"/> >> </OPERATOR> >> <OPERATOR name="Sminus" matrixelement="1/2"> >> <CHANGE quantumnumber="Nup" change="-1"/> >> <CHANGE quantumnumber="Ndown" change="+1"/> >> </OPERATOR> >> <OPERATOR name="Sz" matrixelement="(Nup-Ndown)/2"/> >> <OPERATOR name="Nup" matrixelement="Nup"/> >> <OPERATOR name="Ndown" matrixelement="Ndown"/> >> <OPERATOR name="c_down" matrixelement="1"> >> <CHANGE quantumnumber="Ndown" change="-1"/> >> </OPERATOR> >> <OPERATOR name="cdag_down" matrixelement="1"> >> <CHANGE quantumnumber="Ndown" change="1"/> >> </OPERATOR> >> <OPERATOR name="c_up" matrixelement="1"> >> <CHANGE quantumnumber="Nup" change="-1"/> >> </OPERATOR> >> <OPERATOR name="cdag_up" matrixelement="1"> >> <CHANGE quantumnumber="Nup" change="1"/> >> </OPERATOR> >> <OPERATOR name="n" matrixelement="Nup+Ndown"/> >> <OPERATOR name="n_up" matrixelement="Nup"/> >> <OPERATOR name="n_down" matrixelement="Ndown"/> >> </SITEBASIS> >> >> <BASIS name="fermion"> >> <SITEBASIS ref="fermion"/> >> <CONSTRAINT quantumnumber="Nup" value="Nup_total"/> >> <CONSTRAINT quantumnumber="Ndown" value="Ndown_total"/> >> </BASIS> >> >> <BONDOPERATOR name="one-bond" source="1" target="2"> >> I*(cdag_up(1)*c_up(2)-cdag_up(2)*c_up(1)) >> </BONDOPERATOR> >> >> <HAMILTONIAN name="simple"> >> <PARAMETER name="U" default="1"/> >> <BASIS ref="fermion"/> >> <SITETERM site="i"> >> U*n_up(i)*n_down(i) >> </SITETERM> >> </HAMILTONIAN> >> >> </MODELS> >> >> -- >> Sebastian Huber >> Department of Condensed Matter Physics >> The Weizmann Institute of Science >> 76100 Rehovot >> >> Web: http://www.weizmann.ac.il/~huber >> Tel: +972 8 934 4258 >> >

2 1

BONDOPERATOR measurement question
by Sebastian Huber 22 Nov '11

22 Nov '11

Dear all, I have an inhomogeneous lattice problem and I try to measure a "local" bondoperator, e.g., the current over one bond only. I get the following error from the sparsediag application Did not parse to end of string '(one-bond)/4(i)' when I try to do a "MEASURE_AVERAGE". I get no error if I use "MEASURE_LOCAL". I am a bit puzzled that "MEASURE_LOCAL" does something all together as my operator contains more than one site, which is in violation to a condition on local measurements on the alps webpage. I paste a minimal (silly) set of configuration files which produces this error below. Help would be greatly appreciated. Sebastian lanczosconf: MODEL="simple"; LATTICE="2x2" LATTICE_LIBRARY="graph.xml" MODEL_LIBRARY="model.xml" COMPLEX=true; t=-1; Nup_total=1; Ndown_total=1; CONSERVED_QUANTUMNUMBERS="N,Sz"; NUMBER_EIGENVALUES=2; TRANSLATION_SYMMETRY=false; MEASURE_AVERAGE[one-bond]=one-bond; {U=1;} graph.xml: <LATTICES> <GRAPH name="2x2" dimension="2" vertices="4" edges="9"> <VERTEX id="1" type="0"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX id="2" type="0"><COORDINATE>1 0</COORDINATE></VERTEX> <VERTEX id="3" type="0"><COORDINATE>0 1</COORDINATE></VERTEX> <VERTEX id="4" type="0"><COORDINATE>1 1</COORDINATE></VERTEX> <EDGE source="1" target="2" id="1" type="1" vector="1 0"/> <EDGE source="2" target="1" id="2" type="2" vector="1 0"/> <EDGE source="3" target="4" id="3" type="3" vector="1 0"/> <EDGE source="4" target="3" id="4" type="4" vector="1 0"/> <EDGE source="1" target="3" id="5" type="5" vector="0 1"/> <EDGE source="2" target="4" id="6" type="6" vector="0 1"/> <EDGE source="3" target="1" id="7" type="7" vector="0 1"/> <EDGE source="4" target="2" id="8" type="8" vector="0 1"/> <EDGE source="1" target="4" id="9" type="9" vector="1 1"/> </GRAPH> </LATTICES> model.xml <MODELS> <SITEBASIS name="fermion"> <QUANTUMNUMBER name="Nup" min="0" max="1" type="fermionic"/> <QUANTUMNUMBER name="Ndown" min="0" max="1" type="fermionic"/> <OPERATOR name="Splus" matrixelement="1/2"> <CHANGE quantumnumber="Nup" change="1"/> <CHANGE quantumnumber="Ndown" change="-1"/> </OPERATOR> <OPERATOR name="Sminus" matrixelement="1/2"> <CHANGE quantumnumber="Nup" change="-1"/> <CHANGE quantumnumber="Ndown" change="+1"/> </OPERATOR> <OPERATOR name="Sz" matrixelement="(Nup-Ndown)/2"/> <OPERATOR name="Nup" matrixelement="Nup"/> <OPERATOR name="Ndown" matrixelement="Ndown"/> <OPERATOR name="c_down" matrixelement="1"> <CHANGE quantumnumber="Ndown" change="-1"/> </OPERATOR> <OPERATOR name="cdag_down" matrixelement="1"> <CHANGE quantumnumber="Ndown" change="1"/> </OPERATOR> <OPERATOR name="c_up" matrixelement="1"> <CHANGE quantumnumber="Nup" change="-1"/> </OPERATOR> <OPERATOR name="cdag_up" matrixelement="1"> <CHANGE quantumnumber="Nup" change="1"/> </OPERATOR> <OPERATOR name="n" matrixelement="Nup+Ndown"/> <OPERATOR name="n_up" matrixelement="Nup"/> <OPERATOR name="n_down" matrixelement="Ndown"/> </SITEBASIS> <BASIS name="fermion"> <SITEBASIS ref="fermion"/> <CONSTRAINT quantumnumber="Nup" value="Nup_total"/> <CONSTRAINT quantumnumber="Ndown" value="Ndown_total"/> </BASIS> <BONDOPERATOR name="one-bond" source="1" target="2"> I*(cdag_up(1)*c_up(2)-cdag_up(2)*c_up(1)) </BONDOPERATOR> <HAMILTONIAN name="simple"> <PARAMETER name="U" default="1"/> <BASIS ref="fermion"/> <SITETERM site="i"> U*n_up(i)*n_down(i) </SITETERM> </HAMILTONIAN> </MODELS> -- Sebastian Huber Department of Condensed Matter Physics The Weizmann Institute of Science 76100 Rehovot Web: http://www.weizmann.ac.il/~huber Tel: +972 8 934 4258

2 1

Vector average
by Michal Maik 21 Nov '11

21 Nov '11

Hi all I was wondering how can I extract a single part of a vector average, for example a single value from the Structure Factor. I am looking along the lines of how would I write the xml file to get the results out. Below I have a simple plot file for Magnetization. How could I change it to get out one part of a vector average, such as from the structure factor? <?xml version="1.0" encoding="UTF-8"?> <plot name="Magnetization Density"> <legend show="true"/> <xaxis label="chemical potential" type="PARAMETER" name="mu"/> <yaxis label="Magnetization Density" type="SCALAR_AVERAGE"/> <set label="2D triangular lattice"/> </plot> In advance thank you for the help. Best Regards Michal Maik

3 3

Re: [ALPS-users] Fwd: alps 2.0.2 - failed tests
by Mateusz Łącki 20 Nov '11

20 Nov '11

I am using gcc version 4.6.1 20110908 (Red Hat 4.6.1-9) (GCC) for c++ and ifort for fortran Version 12.0.0 (trief gfortran before but i got errors that some hdf5 modules are not gfortran compatible) Regards, Mateusz >This all seems to boil down to one issue "Illegal term in expression", a runtime error in the parser. Which compiler are >you using? > >Matthias

2 1

Re: [ALPS-users] alps 2.0.2 - compile error
by Mateusz Łącki 20 Nov '11

20 Nov '11

I attach the requested file (i changed the name) Regards Mateusz Lacki

2 5

alps 2.0.2 - compile error
by Mateusz Łącki 19 Nov '11

19 Nov '11

Dear All, I am trying to compile the newest alps release. I managed to get up to 62% of the compilation (on two machines - one running Fedora, the other ubuntu 10.04) up to exactly the same stage: Scanning dependencies of target jacobidavidson2 [ 62%] Building CXX object example/ietl/CMakeFiles/jacobidavidson2.dir/jacobidavidson2.cpp.o In file included from /home/stefan/alps/alps-2.0.2-r5790-src-with-boost/alps/src/boost/numeric/bindings/lapack/computational.hpp:39:0, from /home/stefan/alps/alps-2.0.2-r5790-src-with-boost/alps/src/ietl/jd.h:51, from /home/stefan/alps/alps-2.0.2-r5790-src-with-boost/alps/example/ietl/jacobidavidson2.cpp:52: /home/stefan/alps/alps-2.0.2-r5790-src-with-boost/alps/src/boost/numeric/bindings/lapack/computational/getrs.hpp: In static member function ‘static std::ptrdiff_t boost::numeric::bindings::lapack::getrs_impl<Value>::invoke(const MatrixA&, const VectorIPIV&, MatrixB&) [with MatrixA = boost::numeric::ublas::matrix<double, boost::numeric::ublas::basic_column_major<>, boost::numeric::ublas::unbounded_array<double, std::allocator<double> > >, VectorIPIV = std::vector<int>, MatrixB = boost::numeric::ublas::vector<double>, Value = double, std::ptrdiff_t = long int]’: /home/stefan/alps/alps-2.0.2-r5790-src-with-boost/alps/src/boost/numeric/bindings/lapack/computational/getrs.hpp:247:51: instantiated from ‘std::ptrdiff_t boost::numeric::bindings::lapack::getrs(const MatrixA&, const VectorIPIV&, MatrixB&) [with MatrixA = boost::numeric::ublas::matrix<double, boost::numeric::ublas::basic_column_major<>, boost::numeric::ublas::unbounded_array<double, std::allocator<double> > >, VectorIPIV = std::vector<int>, MatrixB = boost::numeric::ublas::vector<double>, std::ptrdiff_t = long int]’ /home/stefan/alps/alps-2.0.2-r5790-src-with-boost/alps/src/ietl/jd.h:313:9: instantiated from ‘void ietl::solver::left_prec_solver<SOLV, MATRIX, VS, PREC>::operator()(ietl::solver::left_prec_solver<SOLV, MATRIX, VS, PREC>::real_type, ietl::solver::left_prec_solver<SOLV, MATRIX, VS, PREC>::vector_type&, ietl::solver::left_prec_solver<SOLV, MATRIX, VS, PREC>::vector_type&, IT&) [with IT = ietl::jd_iteration<double>, SOLV = ietl::gmres_wrapper, MATRIX = boost::numeric::ublas::symmetric_matrix<double, boost::numeric::ublas::basic_lower<> >, VS = ietl::vectorspace<boost::numeric::ublas::vector<double> >, PREC = jacobi_prec<boost::numeric::ublas::symmetric_matrix<double, boost::numeric::ublas::basic_lower<> >, double>, ietl::solver::left_prec_solver<SOLV, MATRIX, VS, PREC>::real_type = double, ietl::solver::left_prec_solver<SOLV, MATRIX, VS, PREC>::vector_type = boost::numeric::ublas::vector<double>]’ /home/stefan/alps/alps-2.0.2-r5790-src-with-boost/alps/src/ietl/jd.h:633:13: instantiated from ‘void ietl::jd<MATRIX, VS>::jdqr(SOLVER&, IT&, GEN&, ietl::jd<MATRIX, VS>::size_type, bool) [with SOLVER = ietl::solver::left_prec_solver<ietl::gmres_wrapper, boost::numeric::ublas::symmetric_matrix<double, boost::numeric::ublas::basic_lower<> >, ietl::vectorspace<boost::numeric::ublas::vector<double> >, jacobi_prec<boost::numeric::ublas::symmetric_matrix<double, boost::numeric::ublas::basic_lower<> >, double> >, IT = ietl::jd_iteration<double>, GEN = boost::random::lagged_fibonacci_01_engine<double, 48, 44497u, 21034u>, MATRIX = boost::numeric::ublas::symmetric_matrix<double, boost::numeric::ublas::basic_lower<> >, VS = ietl::vectorspace<boost::numeric::ublas::vector<double> >, ietl::jd<MATRIX, VS>::size_type = long unsigned int]’ /home/stefan/alps/alps-2.0.2-r5790-src-with-boost/alps/src/ietl/jd.h:403:17: instantiated from ‘void ietl::jd<MATRIX, VS>::eigensystem(IT&, GEN&, ietl::jd<MATRIX, VS>::size_type, PREC&, SOLV, bool) [with SOLV = ietl::gmres_wrapper, IT = ietl::jd_iteration<double>, GEN = boost::random::lagged_fibonacci_01_engine<double, 48, 44497u, 21034u>, PREC = jacobi_prec<boost::numeric::ublas::symmetric_matrix<double, boost::numeric::ublas::basic_lower<> >, double>, MATRIX = boost::numeric::ublas::symmetric_matrix<double, boost::numeric::ublas::basic_lower<> >, VS = ietl::vectorspace<boost::numeric::ublas::vector<double> >, ietl::jd<MATRIX, VS>::size_type = long unsigned int]’ /home/stefan/alps/alps-2.0.2-r5790-src-with-boost/alps/example/ietl/jacobidavidson2.cpp:136:54: instantiated from here /home/stefan/alps/alps-2.0.2-r5790-src-with-boost/alps/src/boost/numeric/bindings/lapack/computational/getrs.hpp:225:69: error: no matching function for call to ‘getrs(order, trans, boost::numeric::bindings::result_of::size_row_op<const boost::numeric::ublas::matrix<double, boost::numeric::ublas::basic_column_major<>, boost::numeric::ublas::unbounded_array<double, std::allocator<double> > >, boost::numeric::bindings::tag::no_transpose>::type, boost::numeric::bindings::result_of::size_column<boost::numeric::ublas::vector<double> >::type, boost::numeric::bindings::result_of::begin_value<const boost::numeric::ublas::matrix<double, boost::numeric::ublas::basic_column_major<>, boost::numeric::ublas::unbounded_array<double, std::allocator<double> > > >::type, boost::numeric::bindings::result_of::stride_major<const boost::numeric::ublas::matrix<double, boost::numeric::ublas::basic_column_major<>, boost::numeric::ublas::unbounded_array<double, std::allocator<double> > > >::type, boost::numeric::bindings::result_of::begin_value<const std::vector<int> >::type, boost::numeric::bindings::result_of::begin_value<boost::numeric::ublas::vector<double> >::type, boost::numeric::bindings::result_of::stride_major<boost::numeric::ublas::vector<double> >::type)’ /home/stefan/alps/alps-2.0.2-r5790-src-with-boost/alps/src/boost/numeric/bindings/lapack/computational/getrs.hpp:225:69: note: candidates are: /home/stefan/alps/alps-2.0.2-r5790-src-with-boost/alps/src/boost/numeric/bindings/lapack/computational/getrs.hpp:120:23: note: template<class Order, class Trans> std::ptrdiff_t boost::numeric::bindings::lapack::detail::getrs(Order, Trans, fortran_int_t, fortran_int_t, const float*, fortran_int_t, const fortran_int_t*, float*, fortran_int_t) /home/stefan/alps/alps-2.0.2-r5790-src-with-boost/alps/src/boost/numeric/bindings/lapack/computational/getrs.hpp:136:23: note: template<class Order, class Trans> std::ptrdiff_t boost::numeric::bindings::lapack::detail::getrs(Order, Trans, fortran_int_t, fortran_int_t, const double*, fortran_int_t, const fortran_int_t*, double*, fortran_int_t) /home/stefan/alps/alps-2.0.2-r5790-src-with-boost/alps/src/boost/numeric/bindings/lapack/computational/getrs.hpp:152:23: note: template<class Order, class Trans> std::ptrdiff_t boost::numeric::bindings::lapack::detail::getrs(Order, Trans, fortran_int_t, fortran_int_t, const std::complex<float>*, fortran_int_t, const fortran_int_t*, std::complex<float>*, fortran_int_t) /home/stefan/alps/alps-2.0.2-r5790-src-with-boost/alps/src/boost/numeric/bindings/lapack/computational/getrs.hpp:169:23: note: template<class Order, class Trans> std::ptrdiff_t boost::numeric::bindings::lapack::detail::getrs(Order, Trans, fortran_int_t, fortran_int_t, const std::complex<double>*, fortran_int_t, const fortran_int_t*, std::complex<double>*, fortran_int_t) make[2]: *** [example/ietl/CMakeFiles/jacobidavidson2.dir/jacobidavidson2.cpp.o] Error 1 make[1]: *** [example/ietl/CMakeFiles/jacobidavidson2.dir/all] Error 2 it looks like a typical compile error. Am I doing something wrong? Is it missing dependency? I did not do anything nonstandard - I followed instructions from the webpage. Regards, Mateusz Lacki

2 1