Hello, I asked once but got little reply. So may be this time .... :-) If I want to get <n_i> in Bose-Hubbard (i-numbers sites) I put MEASURE_LOCAL[density]=n; (either in QMC or diagonalizations or DMRG).. Can I get easily e.g. <n_i^2> without recompiling applications ? Can I get both <n_i> and <n_i^2> from the same run ? Can I get e.g. <n_in_j> ? best wishes Kuba __________________________________________________________________