Merge pull request #3674 from JuliaReach/schillic/copy

Type-preserving `copy` for `Universe` and `ZeroSet`

src/Sets/Universe/UniverseModule.jl

@@ -15,6 +15,7 @@ using ReachabilityBase.Require: require
                         translate, translate!, cartesian_product, intersection
 @reexport import ..LazySets: constrained_dimensions, linear_map_inverse,
                              tosimplehrep
+import Base: copy
 @reexport using ..API
 
 export Universe
@@ -25,6 +26,7 @@ include("an_element.jl")
 include("complement.jl")
 include("constraints.jl")
 include("constraints_list.jl")
+include("copy.jl")
 include("diameter.jl")
 include("dim.jl")
 include("isbounded.jl")

src/Sets/Universe/copy.jl

@@ -0,0 +1,4 @@
+# since there is only a single instance (per dimensionality), we do not copy
+function copy(U::Universe)
+    return U
+end

src/Sets/ZeroSet/ZeroSetModule.jl

@@ -10,12 +10,14 @@ using ReachabilityBase.Require: require
 @reexport import ..API: dim, isoperationtype, rand, rectify, reflect, ∈,
                         linear_map, scale, scale!, ρ, translate, minkowski_sum
 @reexport import ..LazySets: element
+import Base: copy
 @reexport using ..API
 
 export ZeroSet
 
 include("ZeroSet.jl")
 
+include("copy.jl")
 include("dim.jl")
 include("isoperationtype.jl")
 include("rand.jl")

src/Sets/ZeroSet/copy.jl

@@ -0,0 +1,4 @@
+# since there is only a single instance (per dimensionality), we do not copy
+function copy(Z::ZeroSet)
+    return Z
+end

test/Sets/Universe.jl

@@ -30,6 +30,10 @@ for N in [Float64, Rational{Int}, Float32]
     # dim
     @test dim(U) == 2
 
+    # copy
+    U2 = copy(U)
+    @test U2 isa Universe{N} && dim(U2) == 2
+
     # support function and support vector
     @test ρ(N[0, 1], U) == N(Inf)
     @test ρ(N[0, 0], U) == zero(N)

test/Sets/ZeroSet.jl

@@ -5,6 +5,10 @@ for N in [Float64, Rational{Int}, Float32]
     Z = ZeroSet{N}(2)
     B = BallInf(ones(N, 2), N(1))
 
+    # copy
+    Z2 = copy(Z)
+    @test Z2 isa ZeroSet{N} && dim(Z2) == 2
+
     # isoperationtype
     @test !isoperationtype(ZeroSet)