Merge pull request #1156 from JuliaReach/mforets/promote_linear_map

Use promotion for linear map of polyhedra

src/AbstractPolyhedron_functions.jl

@@ -318,6 +318,18 @@ function linear_map(M::AbstractMatrix{N},
     end
 end
 
+# handle different numeric types
+function linear_map(M::AbstractMatrix{NM},
+                    P::AbstractPolyhedron{NP};
+                    kwargs...) where {NM<:Real, NP<:Real}
+    N = promote_type(NM, NP)
+    if N != NP
+        error("conversion between numeric types of polyhedra not implemented yet (see #1181)")
+    else
+        return linear_map(N.(M), P; kwargs...)
+    end
+end
+
 function _linear_map_vrep(M::AbstractMatrix{N}, P::AbstractPolyhedron{N}) where {N<:Real}
     if !isbounded(P)
         throw(ArgumentError("the linear map in vertex representation for an unbounded set is not defined"))

src/HPolytope.jl

@@ -55,7 +55,6 @@ HPolytope{N}(A::AbstractMatrix{N}, b::AbstractVector{N};
              check_boundedness::Bool=false) where {N<:Real} =
     HPolytope(A, b; check_boundedness=check_boundedness)
 
-
 # --- LazySet interface functions ---
 
 

test/unit_Polytope.jl

@@ -135,6 +135,11 @@ for N in [Float64, Rational{Int}, Float32]
     L4 = linear_map(M, P, check_invertibility=false)
     @test L4 isa VPolygon
 
+    # linear map for mixed types
+    M = [2 1; 0 1] # Int's
+    LM = linear_map(M, P)
+    @test LM isa HPolytope{N}
+
     # -----
     # V-rep
     # -----