diff --git a/src/geometries/polytopes/wedge.jl b/src/geometries/polytopes/wedge.jl index 2ef0a62c2..579dbf080 100644 --- a/src/geometries/polytopes/wedge.jl +++ b/src/geometries/polytopes/wedge.jl @@ -23,8 +23,9 @@ function (wedge::Wedge)(u, v, w) throw(DomainError((u, v, w), "wedge(u, v, w) is not defined for u, v, w outside [0, 1]³.")) end a1, a2, a3, b1, b2, b3 = vertices(wedge) - a = Triangle(a1, a2, a3) - b = Triangle(b1, b2, b3) - s = Segment(a(T(u), T(v)), b(T(u), T(v))) + a = Quadrangle(a1, b1, b2, a2) + b = Quadrangle(a1, b1, b3, a3) + uv = T(u), T(v) + s = Segment(a(uv...), b(uv...)) s(T(w)) -end \ No newline at end of file +end diff --git a/test/polytopes.jl b/test/polytopes.jl index 554635cae..3e5da0c4a 100644 --- a/test/polytopes.jl +++ b/test/polytopes.jl @@ -978,7 +978,7 @@ end @test crs(w) <: Cartesian{NoDatum} @test Meshes.lentype(w) == ℳ @test volume(w) ≈ T(1 / 2) * u"m^3" - @test w(T(0), T(0), T(1)) == vertices(w)[4] + @test w(T(1), T(1), T(1)) == vertices(w)[6] @test_throws DomainError w(T(0), T(0), T(1.5)) m = boundary(w) @test m isa Mesh