Intruduce truncation of blocks

TODO: specify the value of DensityMatrixCutoff externally
j-otsuki · Feb 26, 2018 · c3a8d7c · c3a8d7c
1 parent 855c11a
commit c3a8d7c
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/c++/pomerol_ed.cpp b/c++/pomerol_ed.cpp
@@ -254,6 +254,8 @@ namespace pomerol2triqs {
       rho.reset(new Pomerol::DensityMatrix(*states_class, *matrix_h, beta));
       rho->prepare();
       rho->compute();
+      double DensityMatrixCutoff = 1e-10;  // TODO : specify externally
+      rho->truncateBlocks(DensityMatrixCutoff, verbose);
     }
   }