0-distances bug fxed

docs/source/pecuzal_embedding.py

@@ -952,6 +952,9 @@ def uzal_cost_pecuzal(Y, Y_trial, Tw, K = 3, theiler = 1, norm = 'euclidean', ec
         # increased
 
         dist = compute_L_decrease(E2, E2_trial, eps2, eps2_trial, cnt+1, NN)
+        if math.isnan(dist):
+            raise Exception('Computed 0-distances. You might use model-data, thus try to add minimal additive noise to the signal you wish to embed and try again.')
+
         if dist > dist_former and dist_former<0:
             break
         else:

setup.py

@@ -5,7 +5,7 @@
 
 setuptools.setup(
     name="pecuzal_embedding", 
-    version="1.3.0",
+    version="1.3.1",
     author="K.H.Kraemer",
     author_email="[email protected]",
     description="PECUZAL automatic embedding of uni- and multivariate time series",

src/pecuzal_embedding.py

@@ -952,6 +952,9 @@ def uzal_cost_pecuzal(Y, Y_trial, Tw, K = 3, theiler = 1, norm = 'euclidean', ec
         # increased
 
         dist = compute_L_decrease(E2, E2_trial, eps2, eps2_trial, cnt+1, NN)
+        if math.isnan(dist):
+            raise Exception('Computed 0-distances. You might use model-data, thus try to add minimal additive noise to the signal you wish to embed and try again.')
+
         if dist > dist_former and dist_former<0:
             break
         else:

tests/pecuzal_embedding.py

@@ -952,6 +952,9 @@ def uzal_cost_pecuzal(Y, Y_trial, Tw, K = 3, theiler = 1, norm = 'euclidean', ec
         # increased
 
         dist = compute_L_decrease(E2, E2_trial, eps2, eps2_trial, cnt+1, NN)
+        if math.isnan(dist):
+            raise Exception('Computed 0-distances. You might use model-data, thus try to add minimal additive noise to the signal you wish to embed and try again.')
+
         if dist > dist_former and dist_former<0:
             break
         else: