use improved high numerical precision linspace from Base

KernelDensity.jl rolled its own range generator which suffers from similar flaws as Base's linspace prior to being fixed in JuliaLang/julia#18777. This closes JuliaStats#39.
tlnagy · Jun 3, 2018 · 8a6b97c · 8a6b97c
1 parent bf4191d
commit 8a6b97c
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 2 deletions.
diff --git a/src/univariate.jl b/src/univariate.jl
@@ -66,8 +66,7 @@ function kde_range(boundary::Tuple{Real,Real}, npoints::Int)
  lo, hi = boundary
  lo < hi || error("boundary (a,b) must have a < b")
 
- step = (hi - lo) / (npoints-1)
- lo:step:hi
+ linspace(lo, hi, npoints)
 end
 
 struct UniformWeights{N} end

diff --git a/test/univariate.jl b/test/univariate.jl
@@ -13,6 +13,8 @@ end
 
 r = kde_range((-2.0,2.0), 128)
 @test step(r) > 0
+r = kde_range((0.12698109160784082, 0.9785547869337731), 256)
+@test length(r) == 256
 
 for X in ([0.0], [0.0,0.0], [0.0,0.5], [-0.5:0.1:0.5;])
  w = default_bandwidth(X)