Hi Dominique,
I think the only problem is that you are using an old version of the code. The last “official” ALPS release was published before we finalized our paper. To get a more recent version you can download the sources attached to the paper itself, or, better, simply rely on the latest /trunk version: svn co https://alps.comp-phys.org/svn/alps1/trunk/alps https://alps.comp-phys.org/svn/alps1/trunk/alps
Regarding the additional output. Yes, the actual parameter included more measurements that we eventually didn’t use in the paper. To enable those just add the following lines to parms.py p['MEASURE_LOCAL[Local density up]'] = 'n_up' p['MEASURE_LOCAL[Local density down]'] = 'n_down' p['MEASURE_HALF_CORRELATIONS[dens corr down-down]'] = 'n_down:n_down' Note that Energy^2 is actually coming together with the energy variance.
Cheers, Michele
-- ETH Zurich Michele Dolfi Institute for Theoretical Physics HIT G 32.4 Wolfgang-Pauli-Str. 27 8093 Zurich Switzerland
dolfim@phys.ethz.ch www.itp.phys.ethz.ch
+41 44 633 78 56 phone +41 44 633 11 15 fax
On Jul 9, 2015, at 12:06 PM, Dominique Delande Dominique.Delande@lkb.upmc.fr wrote:
Dear all,
I am trying to use mps_optim and have serious problems. I try to run Example 1 (itinerant fermions in the Hubbard ladder) presented in M. Dolfi et al, CPC 185, 3430 (2014). I downloaded the mmcs1.zip file which should contain all scripts and results. When I run the python script for Example 1, my results are quite different from what they should be:
- No "EnergyVariance" variable in the result files sim.taskX.out.h5
- All "pair field" results are zero, or vanishingly small (of the order of 1.e-300).
- In the result files extracted from mmcs1.zip (e.g. sim.task1.out.h5),
there are additional results such as "Energy^2", "Local density down", "dens corr down-down", etc.
- In the result files extracted from mmcs1.zip (e.g. sim.task1.out.h5),
in the "parameters" section, there are things like MEASURE_HALF_NN_CORRELATIONS which are not defined in the simulation parameters (most probably it it the same as MEASURE_HALF_CORRELATIONS).
I thus have the impression that the results files provided in mmcs1.zip have NOT been obtained using the provided python scripts. More importantly, how to obtain correct "EnergyVariance" and "pair field" results ?
I am using alps-2.2.b3-r7462, compiled from source on a Fedora 21 x86_64 machine, using the GNU compilers and Intel MKL libraries.
Thanks in advance for your help.
Dominique
-- Dominique Delande (Dominique.Delande@lkb.upmc.fr) Laboratoire Kastler-Brossel - Case 74 - Universite P. et M. Curie 4, place Jussieu, F-75252 Paris Cedex 05, FRANCE Phone : +33 1 44 27 27 97 - Fax : +33 1 44 27 38 45 Acces : Tour 13, Couloir 12-13, 3eme etage - Bureau 316