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


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


+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