improve auto-epsilon for low-res data

JuliaDynamics · Apr 1, 2021 · db6ed63 · db6ed63
1 parent 565e261
commit db6ed63
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 4 deletions.
diff --git a/Project.toml b/Project.toml
@@ -1,7 +1,7 @@
 name = "ChaosTools"
 uuid = "608a59af-f2a3-5ad4-90b4-758bdf3122a7"
 repo = "https://github.com/JuliaDynamics/ChaosTools.jl.git"
-version = "1.25.1"
+version = "1.25.2"
 
 [deps]
 Combinatorics = "861a8166-3701-5b0c-9a16-15d98fcdc6aa"

diff --git a/src/dimensions/linear_regions.jl b/src/dimensions/linear_regions.jl
@@ -163,8 +163,14 @@ function estimate_boxsizes(
 
     mi, ma = minmaxima(A)
     max_d = maximum(ma - mi)
-
     min_d, _ = minimum_pairwise_distance(A)
+    if min_d == 0
+        @warn(
+        "Minimum distance in the dataset is zero! Probably because of having data "*
+        "with low resolution, or duplicate data points. Setting to `d₊/base^4` for now.")
+        min_d = max_d/(base^4)
+    end
+
     lower = log(base, min_d)
     upper = log(base, max_d)
 
@@ -175,12 +181,10 @@ function estimate_boxsizes(
         "Automatic boxsize determination was inappropriate: `lower+w` was found ≥ than "*
         "`upper+z`. Returning `base .^ range(lower, upper; length = k)`. "*
         "Please adjust keywords or provide a bigger dataset.")
-
         εs = float(base) .^ range(lower, upper; length = k)
     else
         εs = float(base) .^ range(lower+w, upper+z; length = k)
     end
-    @assert issorted(εs)
     return εs
 end