November 2017 - Comp-phys-alps-users

measuring truncation errors in mps_optim
by Robertson Esperanza 07 Dec '17

07 Dec '17

Thank you for all the help so far, now I'm trying to check for convergence of my results. First, is it enough to just check how the truncation error changes with sweeps (at a certain MAXSTATE value) or I need to both check the iteration history of the ground state energy and the truncation error? Second, checking the supplementary materials provided (Computer Physics Communications 185 (2014) 3430–3440) for the Hubbard ladder, why are there 381 values of energy and truncation weights for every iteration? Why is the total truncation error for that iteration is the sum of all the 381 truncation weight values (based on iteration.py)? How is truncation error defined (e.g. is it 1 - sum of all reduced density eigenvalues corresponding to all the picked states?) ? As far as I understand the DMRG algorithm every iteration/sweep has a corresponding bond dimension in which as you progress through iterations (until you reach the "SWEEPS"th iteration), bond dimension increases until you reach "MAXSTATES" at the "SWEEPS"th iteration. Regards, Robertson Esperanza

4 10

mpirun
by Tadeusz Wasiutyński 22 Nov '17

22 Nov '17

I have installed alps with Conda. Now I can not figure out how to run it mpi mode: mpirun not found. Previously I always build from source and put some parameters in cmake command. What can I do now? -- Tadeusz Wasiutyński

2 1

Questions about MEASURE_LOCAL_AT in MPS-DMRG
by Amir M.Aghaei 22 Nov '17

22 Nov '17

Hi, This is the documentation for 'MEASURE_LOCAL_AT' MEASURE_LOCAL_AT[NAME] > Syntax for the value: > [image: local_at_value.svg] > Defines a new measurement called NAME where the sequence of operators > op_1:...:op_n is applied to all tuples of indices (of length *n* like the > number of operators) listed after the vertical bar symbol. > I have three question: Q1: Are the indices, (i_11, ...i_1n), the vertex ids or vertex coordinates or the index of the vertex in some array, or sth else? Q2: I understand that the restriction i_1 < i_2 < i_3 ... is probably there to simplify the MPS contraction but is there a simple way to go around it? Q3: There is no limit on the number of measurement or the number of operators that can be defined, right? Amir

2 1

compiling error in version 2.3.0
by benoît grémaud 22 Nov '17

22 Nov '17

Dear Developers, I am trying to install alps 2.3.0 on a centos 7 from the sources (including boost). cmake -D ALPS_ENABLE_OPENMP=ON -D Boost_ROOT_DIR:PATH=/home/gremaud/alps/alps-2.3.0-src-with-boost/boost/ -D CMAKE_INSTALL_PREFIX=/home/gremaud/alps/alps-2.3.0-install ../alps-2.3.0-src-with-boost/alps the file CMakeCache.txt is attached I obtain the following error: In file included from /usr/lib64/python2.7/site-packages/numpy/core/include/numpy/ndarraytypes.h:1728:0, from /usr/lib64/python2.7/site-packages/numpy/core/include/numpy/ndarrayobject.h:17, from /usr/lib64/python2.7/site-packages/numpy/core/include/numpy/arrayobject.h:15, from /home/gremaud/alps/alps-2.3.0-src-with-boost/alps/src/alps/python/numpy_import.hpp:40, from /home/gremaud/alps/alps-2.3.0-src-with-boost/alps/src/alps/ngs/detail/get_numpy_type.hpp:38, from /home/gremaud/alps/alps-2.3.0-src-with-boost/alps/src/alps/ngs/detail/paramvalue_reader.hpp:37, from /home/gremaud/alps/alps-2.3.0-src-with-boost/alps/src/alps/ngs/detail/paramvalue.hpp:34, from /home/gremaud/alps/alps-2.3.0-src-with-boost/alps/src/alps/ngs/params.hpp:33, from /home/gremaud/alps/alps-2.3.0-src-with-boost/alps/src/alps/ngs/api.hpp:32, from /home/gremaud/alps/alps-2.3.0-src-with-boost/alps/src/alps/ngs/lib/api.cpp:28: /usr/lib64/python2.7/site-packages/numpy/core/include/numpy/npy_deprecated_api.h:11:2: warning: #warning "Using deprecated NumPy API, disable it by #defining NPY_NO_DEPRECATED_API NPY_1_7_API_VERSION" [-Wcpp] #warning "Using deprecated NumPy API, disable it by #defining NPY_NO_DEPRECATED_API NPY_1_7_API_VERSION" ^ /usr/lib64/python2.7/site-packages/numpy/core/include/numpy/npy_deprecated_api.h:23:2: error: #error Should never include npy_deprecated_api directly. #error Should never include npy_deprecated_api directly. ^ make[2]: *** [src/alps/CMakeFiles/alps.dir/ngs/lib/api.cpp.o] Error 1 The version 2.2.b5 is compiling nicely on the same machine, even though there are a lot of similar warnings about "Using deprecated NumPy API", but no error. Thank you for your help, best regards, Benoît -------------------- Researcher at CNRS, Laboratoire Kastler Brossel (Paris, France) and UMI MajuLab (Singapore) Visiting Associate Professor at the Centre for Quantum Technologies (NUS, Singapore)

2 1

HDF5 runtime error when using mps_optim application with 'optimization' parameter set 'singlesite' on CentOS Linux
by Rongyang Sun 19 Nov '17

19 Nov '17

Dear all, The application mps_optim works well on OS X system which is installed by dmg package. But on Linux, it threw the following error when I set the 'optimization' parameter as 'singlesite': ********************************************************************************* terminate called after throwing an instance of 'alps::hdf5::archive_error' what(): HDF5 error: -1 #0 ../../hdf5-1.10.1/src/H5D.c line 145 in H5Dcreate2(): unable to create dataset #1 ../../hdf5-1.10.1/src/H5Dint.c line 490 in H5D__create_named(): unable to create and link to dataset #2 ../../hdf5-1.10.1/src/H5L.c line 1695 in H5L_link_object(): unable to create new link to object #3 ../../hdf5-1.10.1/src/H5L.c line 1939 in H5L_create_real(): can't insert link #4 ../../hdf5-1.10.1/src/H5Gtraverse.c line 867 in H5G_traverse(): internal path traversal failed #5 ../../hdf5-1.10.1/src/H5Gtraverse.c line 639 in H5G_traverse_real(): traversal operator failed #6 ../../hdf5-1.10.1/src/H5L.c line 1731 in H5L_link_cb(): name already exists In /tmp/alps-2.3.0-src-with-boost/alps/src/alps/hdf5/archive.cpp on 133 in resource /opt/alps/lib/libalps.so.2(_ZN4alps4hdf56detail11check_errorEl+0x121) [0x7f44c5184db1] /opt/alps/lib/libalps.so.2(_ZNK4alps4hdf57archive5writeESsPKdSt6vectorImSaImEES6_S6_+0x168c) [0x7f44c526ddec] mps_optim(_ZN4alps4hdf54saveIdSt6vectorIdSaIdEEEEvRNS0_7archiveERKSsRKNS_7numeric6matrixIT_T0_EES2_ImSaImEESH_SH_+0x2bb) [0x51c7bb] mps_optim(_ZN4alps4hdf54saveINS_7numeric6matrixIdSt6vectorIdSaIdEEEESaIS7_EEEvRNS0_7archiveERKSsRKS4_IT_T0_ES4_ImSaImEESJ_SJ_+0x2e1) [0x51ce01] mps_optim(_ZN4alps4hdf5lsIRKSt6vectorINS_7numeric6matrixIdS2_IdSaIdEEEESaIS7_EEEEN5boost10disable_ifINS0_20has_complex_elementsINS_6detail10remove_cvrIT_E4typeEEERNS0_7archiveEE4typeESM_RKNS0_6detail14make_pvp_proxyISH_EE+0x6e) [0x51d3ae] mps_optim(_ZNK12block_matrixIN4alps7numeric6matrixIdSt6vectorIdSaIdEEEE12NU1_templateILi6EiEE4saveINS0_4hdf57archiveEEEvRT_+0x251) [0x571241] mps_optim(_ZN4alps4hdf54saveI12block_matrixINS_7numeric6matrixIdSt6vectorIdSaIdEEEE12NU1_templateILi6EiEEEEvRNS0_7archiveERKSsRKT_S5_ImSaImEESK_SK_+0xb2) [0x571492] mps_optim(_ZN4alps4hdf5lsIRK12block_matrixINS_7numeric6matrixIdSt6vectorIdSaIdEEEE12NU1_templateILi6EiEEEEN5boost10disable_ifINS0_20has_complex_elementsINS_6detail10remove_cvrIT_E4typeEEERNS0_7archiveEE4typeESO_RKNS0_6detail14make_pvp_proxyISJ_EE+0x6e) [0x57188e] mps_optim(_ZNK9MPSTensorIN4alps7numeric6matrixIdSt6vectorIdSaIdEEEE12NU1_templateILi6EiEE4saveINS0_4hdf57archiveEEEvRT_+0x158) [0x571a58] mps_optim(_ZN4alps4hdf54saveI9MPSTensorINS_7numeric6matrixIdSt6vectorIdSaIdEEEE12NU1_templateILi6EiEEEEvRNS0_7archiveERKSsRKT_S5_ImSaImEESK_SK_+0xb2) [0x571be2] mps_optim(_ZN4alps4hdf5lsIRK9MPSTensorINS_7numeric6matrixIdSt6vectorIdSaIdEEEE12NU1_templateILi6EiEEEEN5boost10disable_ifINS0_20has_complex_elementsINS_6detail10remove_cvrIT_E4typeEEERNS0_7archiveEE4typeESO_RKNS0_6detail14make_pvp_proxyISJ_EE+0x6e) [0x571fde] mps_optim() [0x47da4b] /usr/lib64/libgomp.so.1(+0x16435) [0x7f44c249e435] /usr/lib64/libpthread.so.0(+0x7e25) [0x7f44c389ae25] /usr/lib64/libc.so.6(clone+0x6d) [0x7f44c1fa734d] ********************************************************************************* The error is always thrown after the first sweep procedure,i.e., the beginning of the second sweep procedure. I built the ALPS from source code without met any error. So anyone has any ideas? Thanks a lot! Best regards, Rongyang Sun

1 1

How to count dimension of graph defined in Lattice xml library
by Sun Rongyang 10 Nov '17

10 Nov '17

Dear all, The definition method of a lattice in ALPS is delicate. But I meet a problem about dimension counting of a pure graph. On the one hand, I define an pure anisotropic graph as: ************************************************************************************************* <GRAPH name="anisotropic open ladder graph test">  <EDGE tyep="0" source="1" target="2"/> <EDGE tyep="0" source="2" target="3"/> <EDGE tyep="0" source="3" target="4"/> <EDGE tyep="0" source="4" target="5"/> <EDGE tyep="0" source="6" target="7"/> <EDGE tyep="0" source="7" target="8"/> <EDGE tyep="0" source="8" target="9"/> <EDGE tyep="0" source="9" target="10"/> <EDGE tyep="1" source="1" target="6"/> <EDGE tyep="1" source="2" target="7"/> <EDGE tyep="1" source="3" target="8"/> <EDGE tyep="1" source="4" target="9"/> <EDGE tyep="1" source="5" target="10"/> </GRAPH> ************************************************************************************************* On the other hand, I use the default LATTICEGRAPH "open ladder", I copy here: ************************************************************************************************* <LATTICE name="square lattice" dimension="2"> <PARAMETER name="a" default="1"/> <BASIS><VECTOR>a 0</VECTOR><VECTOR>0 a</VECTOR></BASIS> <RECIPROCALBASIS><VECTOR>2*pi/a 0</VECTOR><VECTOR>0 2*pi/a</VECTOR></RECIPROCALBASIS> </LATTICE> <UNITCELL name="anisotropic2d" dimension="2"> <VERTEX/> <EDGE type="0"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="1 0"/></EDGE> <EDGE type="1"><SOURCE vertex="1" offset="0 0"/><TARGET vertex="1" offset="0 1"/></EDGE> </UNITCELL> <LATTICEGRAPH name = "open ladder" vt_name="OpenLadderLattice"> <FINITELATTICE> <LATTICE ref="square lattice"/> <PARAMETER name="W" default="2"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY dimension="1" type="open"/> <BOUNDARY dimension="2" type="open"/> </FINITELATTICE> <UNITCELL ref="anisotropic2d"/> </LATTICEGRAPH> ************************************************************************************************* with parameters L=2, W=5. Will they be treated as the same in ALPS? Or the graph equals as L=5,W=2 open ladder? Thanks a lot! Best regards, Rongyang Sun

2 7