From 5805c7c3aa145e10fab19b3d7aeb58cc8f0251c3 Mon Sep 17 00:00:00 2001
From: Michel Ferrero <mferrero@cpht.polytechnique.fr>
Date: Sat, 8 Mar 2014 13:46:31 +0100
Subject: [PATCH] Fix in off-diagonal elements of density()

This should fix the bug noticed by Leonid in off-diagonal
elements of the density matrix.

  modified:   triqs/gfs/local/functions.cpp
---
 triqs/gfs/local/functions.cpp | 3 +--
 1 file changed, 1 insertion(+), 2 deletions(-)

diff --git a/triqs/gfs/local/functions.cpp b/triqs/gfs/local/functions.cpp
index 37e197df..cf9ff46e 100644
--- a/triqs/gfs/local/functions.cpp
+++ b/triqs/gfs/local/functions.cpp
@@ -55,9 +55,8 @@ namespace triqs { namespace gfs {
    }
   
   for (size_t n1=0; n1<N1;n1++) 
-   for (size_t n2=0; n2<N2;n2++) {
+   for (size_t n2=n1; n2<N2;n2++) {
     dens_part(n1,n2) = dens_part(n1,n2) + (G.mesh().positive_only()? ( real(dens_part(n2,n1)) - I * imag(dens_part(n2,n1)) ) : 0) + dens_tail(n1,n2);
-    // ?? STRANGE ??
     dens_part(n2,n1) = real(dens_part(n1,n2)) - I * imag(dens_part(n1,n2));
    }