Add fallback method for Base.in predicate

JuliaGeometry · Dec 26, 2023 · 3e3a95a · 3e3a95a
1 parent 6d56ab4
commit 3e3a95a
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 1 deletion.
diff --git a/src/predicates/in.jl b/src/predicates/in.jl
@@ -7,7 +7,7 @@
 
 Tells whether or not the `point` is in the `geometry`.
 """
-function Base.in(::Point, ::Geometry) end
+Base.in(p::Point, g::Geometry) = sideof(p, boundary(g)) == IN
 
 Base.in(p₁::Point, p₂::Point) = p₁ == p₂
 

diff --git a/test/predicates.jl b/test/predicates.jl
@@ -27,6 +27,14 @@
     @test isconvex(b)
   end
 
+  @testset "in" begin
+    h = first(CartesianGrid{T}(10, 10, 10))
+    @test P3(0, 0, 0) ∈ h
+    @test P3(0.5, 0.5, 0.5) ∈ h
+    @test P3(-1, 0, 0) ∉ h
+    @test P3(0, 2, 0) ∉ h
+  end
+
   @testset "issubset" begin
     point = P2(0.5, 0.5)
     box = Box(P2(0, 0), P2(1, 1))