diff --git a/cvmatrix/__init__.py b/cvmatrix/__init__.py
index a6221b3..dab2cbb 100644
--- a/cvmatrix/__init__.py
+++ b/cvmatrix/__init__.py
@@ -1 +1 @@
-__version__ = '1.0.2'
+__version__ = '1.0.2.post1'
diff --git a/pyproject.toml b/pyproject.toml
index 6e81738..a30f974 100644
--- a/pyproject.toml
+++ b/pyproject.toml
@@ -1,6 +1,6 @@
 [tool.poetry]
 name = "cvmatrix"
-version = "1.0.2"
+version = "1.0.2.post1"
 description = "Fast computation of possibly centered/scaled training set kernel matrices in a cross-validation setting."
 authors = ["Sm00thix <oleemail@icloud.com>"]
 maintainers = ["Sm00thix <oleemail@icloud.com>"]
diff --git a/tests/naive_cvmatrix.py b/tests/naive_cvmatrix.py
index 805a39b..56678c5 100644
--- a/tests/naive_cvmatrix.py
+++ b/tests/naive_cvmatrix.py
@@ -139,7 +139,7 @@ def _training_matrices(
             X_train = X_train - X_train_mean
         if self.scale_X:
             X_train_std = X_train.std(axis=0, ddof=1, keepdims=True)
-            X_train_std[X_train_std == 0] = 1
+            X_train_std[X_train_std <= self.eps] = 1
             X_train = X_train / X_train_std
         if return_XTY:
             Y_train = self.Y_total[training_indices]
@@ -148,7 +148,7 @@ def _training_matrices(
                 Y_train = Y_train - Y_train_mean
             if self.scale_Y:
                 Y_train_std = Y_train.std(axis=0, ddof=1, keepdims=True)
-                Y_train_std[Y_train_std == 0] = 1
+                Y_train_std[Y_train_std <= self.eps] = 1
                 Y_train = Y_train / Y_train_std
         if return_XTX and return_XTY:
             return X_train.T @ X_train, X_train.T @ Y_train