Fix bug in lattice tools dos function

TRIQS · Nov 11, 2019 · 3e0ed62 · 3e0ed62
1 parent 2ede41a
commit 3e0ed62
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/triqs/lattice/tight_binding.cpp b/triqs/lattice/tight_binding.cpp
@@ -130,7 +130,7 @@ namespace triqs {
         for (int j = 0; j < grid.size(); j++) {
           int a = int((eval(l, j) - epsmin) / deps);
           if (a == int(neps)) a = a - 1;
-          for (int k = 0; k < norb; k++) { rho(a, l) += real(conj(evec(l, k, j)) * evec(l, k, j)); }
+          for (int k = 0; k < norb; k++) { rho(a, k) += real(conj(evec(l, k, j)) * evec(l, k, j)); }
         }
       }
       rho /= grid.size() * deps;