diff --git a/Project.toml b/Project.toml
index 123bd7a..efcc893 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "GeoStatsFunctions"
 uuid = "6771c435-bc22-4842-b0c3-41852a255103"
 authors = ["Júlio Hoffimann <julio.hoffimann@gmail.com> and contributors"]
-version = "0.4.1"
+version = "0.4.2"
 
 [deps]
 Bessels = "0e736298-9ec6-45e8-9647-e4fc86a2fe38"
@@ -33,7 +33,7 @@ GeoTables = "1.21"
 InteractiveUtils = "1.9"
 LinearAlgebra = "1.9"
 Makie = "0.21"
-Meshes = "0.46"
+Meshes = "0.47"
 NearestNeighbors = "0.4"
 Optim = "1.7"
 Printf = "1.9"
diff --git a/test/sampling.jl b/test/sampling.jl
index 3c54e2d..bb35c71 100644
--- a/test/sampling.jl
+++ b/test/sampling.jl
@@ -2,13 +2,13 @@
   γ = GaussianVariogram()
   seg = Segment((0.0, 0.0), (1.0, 1.0))
   ps = variosample(γ, seg) |> collect
-  @test all(p -> Point(0.0, 0.0) ⪯ p ⪯ Point(1.0, 1.0), ps)
+  @test all(p -> Point(0.0, 0.0) ≤ p ≤ Point(1.0, 1.0), ps)
   @test length(ps) == 3
 
   γ = GaussianVariogram()
   quad = Quadrangle((0.0, 0.0), (1.0, 0.0), (1.0, 1.0), (0.0, 1.0))
   ps = variosample(γ, quad) |> collect
-  @test all(p -> Point(0.0, 0.0) ⪯ p ⪯ Point(1.0, 1.0), ps)
+  @test all(p -> Point(0.0, 0.0) ≤ p ≤ Point(1.0, 1.0), ps)
   @test length(ps) == 3 * 3
 
   γ = GaussianVariogram()
@@ -23,7 +23,7 @@
     (0.0, 1.0, 1.0)
   )
   ps = variosample(γ, hex) |> collect
-  @test all(p -> Point(0.0, 0.0, 0.0) ⪯ p ⪯ Point(1.0, 1.0, 1.0), ps)
+  @test all(p -> Point(0.0, 0.0, 0.0) ≤ p ≤ Point(1.0, 1.0, 1.0), ps)
   @test length(ps) == 3 * 3 * 3
 
   # regularization properties