diff --git a/python/w2dyn_cthyb/solver.py b/python/w2dyn_cthyb/solver.py
index 1c39704..7614d7f 100644
--- a/python/w2dyn_cthyb/solver.py
+++ b/python/w2dyn_cthyb/solver.py
@@ -394,7 +394,7 @@ def solve(self, **params_kw):
                     s1 = spins[0]
                     s2 = spins[1]
 
-                    all_zeros = not np.any(ftau[:, b1, s1, b2, s2] > 1e-5)
+                    all_zeros = not np.any(np.abs(ftau[:, b1, s1, b2, s2]) > 1e-5)
 
                     if not all_zeros:
                         components = np.append(components, comp_ind)