fix bug on HyperRectangle affecting DyadicSplitter

Author: maltezfaria

src/Interpolation/hyperrectangle.jl

@@ -194,8 +194,10 @@ function HyperRectangle(els,cube=false)
     if cube # fit a square/cube instead
         w  = maximum(ub-lb)
         xc = (ub + lb) / 2
-        lb = xc .- w/2
-        ub = xc .+ w/2
+        # min/max below helps avoid floating point issues with results with the
+        # cube not being exactly contained in the original rectangle
+        lb = min.(xc .- w/2,lb)
+        ub = max.(xc .+ w/2,ub)
         # TODO: return HyperCube instead
     end
     return HyperRectangle(lb,ub)

test/Trees/clustertree_test.jl

@@ -56,7 +56,7 @@ end
         @test test_cluster_tree(clt)
         splitter = WPB.DyadicSplitter(nmax=32)
         clt = WPB.ClusterTree(points,splitter)
-        @test test_cluster_tree(clt) skip=true
+        @test test_cluster_tree(clt)
     end
 
     @testset "3d" begin
@@ -75,7 +75,7 @@ end
         @test test_cluster_tree(clt)
         splitter   = WPB.DyadicSplitter(nmax=32)
         clt = WPB.ClusterTree(points,splitter)
-        @test test_cluster_tree(clt) skip=true
+        @test test_cluster_tree(clt)
     end
 
     @testset "3d + threads" begin
@@ -95,6 +95,6 @@ end
         @test test_cluster_tree(clt)
         splitter   = WPB.DyadicSplitter(nmax=32)
         clt = WPB.ClusterTree(points,splitter;threads)
-        @test test_cluster_tree(clt) skip=true
+        @test test_cluster_tree(clt)
     end
 end