Allow to construct H-rep from mixed numeric types

src/Convert/convert.jl

@@ -80,7 +80,7 @@ for T in subtypes(AbstractHPolygon, true)
             @assert dim(X) == 2 "set must be two-dimensional for conversion, but it " *
                                 "has dimension $(dim(X))"
             if check_boundedness && !isbounded(X)
-                throw(ArgumentError("expected a bounded set for conversion to `$T`"))
+                throw(ArgumentError("expected a bounded set for conversion to `$($T)`"))
             end
             return $T(constraints_list(X); prune=prune)
         end

src/Sets/HalfSpace/HalfSpaceModule.jl

@@ -104,9 +104,14 @@ function is_tighter_same_dir_2D(c1::HalfSpace,
 end
 
 # TODO: after #2032, #2041 remove use of this function
-_normal_Vector(c::HalfSpace) = HalfSpace(convert(Vector, c.a), c.b)
-_normal_Vector(C::Vector{<:HalfSpace}) = [_normal_Vector(c) for c in C]
 _normal_Vector(P::LazySet) = _normal_Vector(constraints_list(P))
+_normal_Vector(c::HalfSpace) = HalfSpace(convert(Vector, c.a), c.b)
+_normal_Vector(C::Vector{<:HalfSpace{N}}) where {N} = [_normal_Vector(c) for c in C]
+
+function _normal_Vector(C::Vector{<:HalfSpace})
+    N = promote_type([eltype(c) for c in C]...)
+    return [HalfSpace(convert(Vector{N}, c.a), N(c.b)) for c in C]
+end
 
 include("init.jl")
 

test/Sets/Polygon.jl

@@ -738,6 +738,10 @@ rand(HPolygonOpt)
 @test HPolygonOpt() isa HPolygonOpt{Float64,Vector{Float64}}
 @test VPolygon() isa VPolygon{Float64}
 
+# construction from constraints of mixed numeric types
+P = HPolygon([HalfSpace([1.0, 1.0], 1.0), HalfSpace(Float16[1, 1], Float16(1))])
+@test P isa HPolygon{Float64}
+
 # isoperationtype
 @test !isoperationtype(HPolygon)
 @test !isoperationtype(HPolygonOpt)