implement isconvex for Pyramid

src/Meshes.jl

@@ -226,6 +226,7 @@ export
   angles,
   innerangles,
   normal,
+  base,
   ≗,
 
   # multi-geometries

src/geometries/polytopes/pyramid.jl

@@ -11,6 +11,11 @@ A pyramid with points `p1`, `p2`, `p3`, `p4`, `p5`.
 
 nvertices(::Type{<:Pyramid}) = 5
 
+function base(p::Pyramid)
+  a, b, c, d, o = p.vertices
+  Quadrangle(a, b, c, d)
+end
+
 ==(p₁::Pyramid, p₂::Pyramid) = p₁.vertices == p₂.vertices
 
 Base.isapprox(p₁::Pyramid, p₂::Pyramid; atol=atol(lentype(p₁)), kwargs...) =
@@ -20,7 +25,7 @@ function (pyramid::Pyramid)(u, v, w)
   if (u < 0 || u > 1) || (v < 0 || v > 1) || (w < 0 || w > 1)
     throw(DomainError((u, v, w), "pyramid(u, v, w) is not defined for u, v, w outside [0, 1]³."))
   end
-  a, b, c, d, o = vertices(pyramid)
+  a, b, c, d, o = p.vertices
   q = Quadrangle(a, b, c, d)
   s = Segment(q(u, v), o)
   s(w)

src/predicates/isconvex.jl

@@ -66,6 +66,8 @@ _isconvex(p::Polygon, ::Val{3}) = isconvex(proj2D(p))
 
 isconvex(h::Hexahedron) = all(isconvex, boundary(h))
 
+isconvex(p::Pyramid) = isconvex(base(p))
+
 isconvex(m::Multi) = isapproxequal(measure(convexhull(m)), measure(m))
 
 # --------------

test/predicates.jl

@@ -328,6 +328,22 @@ end
   @test isconvex(h0)
   @test isconvex(h1)
   @test isconvex(h2)
+  p = Pyramid(
+    cart(0, 0, 0),
+    cart(1, 0, 0),
+    cart(1, 1, 0),
+    cart(0, 1, 0),
+    cart(0, 0, 1)
+  )
+  @test isconvex(p)
+  p = Pyramid(
+    cart(0, 0, 0),
+    cart(1, 0, 0),
+    cart(0.25, 0.25, 0),
+    cart(0, 1, 0),
+    cart(0, 0, 1)
+  )
+  @test !isconvex(p)
 end
 
 @testitem "issubset" setup = [Setup] begin