Fixed bug in index counter.HEADmaster

1 files changed, 2 insertions, 2 deletions

diff --git a/collectMorseData.hpp b/collectMorseData.hpp
index eddb825..fdee480 100644
--- a/collectMorseData.hpp
+++ b/collectMorseData.hpp
@@ -34,7 +34,7 @@ void collectMorseData(std::string tag, unsigned N, Generator& r, Real ε, Real Î
   unsigned indx = 0;
 
   for (Real λ : eigenS.eigenvalues().real()) {
-    if (λ < 0) indx++;
+    if (λ < 0 && std::abs(λ) > 1e-14) indx++;
   }
 
   pSpinModel M = ReM;
@@ -64,7 +64,7 @@ void collectMorseData(std::string tag, unsigned N, Generator& r, Real ε, Real Î
     indz = 0;
 
     for (Real λ : eigenSz.eigenvalues().real()) {
-      if (λ < 0) indz++;
+      if (λ < 0 && std::abs(λ) > 1e-14) indz++;
     }
   }