diff --git a/mdopt/mps/canonical.py b/mdopt/mps/canonical.py
index ebeab6c..ecd0e10 100644
--- a/mdopt/mps/canonical.py
+++ b/mdopt/mps/canonical.py
@@ -272,7 +272,7 @@ def dense(
 
         if renormalise:
             norm = float(np.linalg.norm(dense, ord=norm))
-            if norm > 1e-12:
+            if norm > 1e-17:
                 dense /= norm
 
         return dense
@@ -315,7 +315,7 @@ def density_mpo(self) -> List[np.ndarray]:
 
         return list(mpo)
 
-    def entanglement_entropy(self, tolerance: float = float(1e-12)) -> np.ndarray:
+    def entanglement_entropy(self, tolerance: float = float(1e-17)) -> np.ndarray:
         """
         Returns the entanglement entropy for bipartitions at each of the bonds.
         """
@@ -702,7 +702,7 @@ def compress_bond(
         self,
         bond: int,
         chi_max: int = int(1e4),
-        cut: float = float(1e-12),
+        cut: float = float(1e-17),
         renormalise: bool = False,
         strategy: str = "svd",
         return_truncation_error: bool = False,
@@ -873,7 +873,7 @@ def compress_bond(
     def compress(
         self,
         chi_max: int = int(1e4),
-        cut: float = float(1e-12),
+        cut: float = float(1e-17),
         renormalise: bool = False,
         strategy: str = "svd",
         return_truncation_errors: bool = False,
diff --git a/mdopt/mps/explicit.py b/mdopt/mps/explicit.py
index 0bcad96..2f093d6 100644
--- a/mdopt/mps/explicit.py
+++ b/mdopt/mps/explicit.py
@@ -292,7 +292,7 @@ def dense(
 
         if renormalise:
             norm = float(np.linalg.norm(dense, ord=norm))
-            if norm > 1e-12:
+            if norm > 1e-17:
                 dense /= norm
 
         return dense
@@ -551,7 +551,7 @@ def compress_bond(
         self,
         bond: int,
         chi_max: int = int(1e4),
-        cut: float = float(1e-12),
+        cut: float = float(1e-17),
         renormalise: bool = False,
         return_truncation_error: bool = False,
     ) -> Tuple["ExplicitMPS", Optional[float]]:
@@ -621,7 +621,7 @@ def compress_bond(
     def compress(
         self,
         chi_max: int = int(1e4),
-        cut: float = float(1e-12),
+        cut: float = float(1e-17),
         renormalise: bool = False,
         return_truncation_errors: bool = False,
     ) -> Tuple["ExplicitMPS", List[Optional[float]]]: