fix 'vertices_list' of zonotopes with immutable center vector

Author: schillic

src/Interfaces/AbstractZonotope.jl

@@ -426,9 +426,10 @@ function _vertices_list_zonotope_2D_positive(c::AbstractVector{N}, G::AbstractMa
     return vlist
 end
 
-function _vertices_list_zonotope_iterative(c::VN, G::AbstractMatrix{N};
-                                           apply_convex_hull::Bool) where {N,VN<:AbstractVector{N}}
+function _vertices_list_zonotope_iterative(c::AbstractVector{N}, G::AbstractMatrix{N};
+                                           apply_convex_hull::Bool) where {N}
     p = size(G, 2)
+    VN = vector_type(typeof(c))
     vlist = Vector{VN}()
     sizehint!(vlist, 2^p)
 
@@ -440,9 +441,9 @@ function _vertices_list_zonotope_iterative(c::VN, G::AbstractMatrix{N};
 end
 
 # special case 2D zonotope of order 1/2
-function _vertices_list_zonotope_2D_order_one_half(c::VN, G::AbstractMatrix{N};
-                                                   apply_convex_hull::Bool) where {N,
-                                                                                   VN<:AbstractVector{N}}
+function _vertices_list_zonotope_2D_order_one_half(c::AbstractVector{N}, G::AbstractMatrix{N};
+                                                   apply_convex_hull::Bool) where {N}
+    VN = vector_type(typeof(c))
     vlist = Vector{VN}(undef, 2)
     g = view(G, :, 1)
     @inbounds begin
@@ -453,9 +454,9 @@ function _vertices_list_zonotope_2D_order_one_half(c::VN, G::AbstractMatrix{N};
 end
 
 # special case 2D zonotope of order 1
-function _vertices_list_zonotope_2D_order_one(c::VN, G::AbstractMatrix{N};
-                                              apply_convex_hull::Bool) where {N,
-                                                                              VN<:AbstractVector{N}}
+function _vertices_list_zonotope_2D_order_one(c::AbstractVector{N}, G::AbstractMatrix{N};
+                                              apply_convex_hull::Bool) where {N}
+    VN = vector_type(typeof(c))
     vlist = Vector{VN}(undef, 4)
     a = [one(N), one(N)]
     b = [one(N), -one(N)]

test/Sets/Zonotope.jl

@@ -407,6 +407,10 @@ for N in @tN([Float64, Float32, Rational{Int}])
                             [N[2, 3, 4], N[2, 3, -2], N[2, -1, 4], N[2, -1, -2],
                              N[0, 3, 4], N[0, 3, -2], N[0, -1, 4], N[0, -1, -2]])
     end
+    # SubArray
+    Z2 = Zonotope(view(c, 1:2), genmat(Z))
+    vlistZ = vertices_list(Z2)
+    @test vlistZ == vertices_list(Z)
 
     # vertices
     res = collect(vertices(Z))