Merge pull request #243 from JuliaReach/mforets/fixes

Fix dim for constrained reset map method + overapprox with an HPolytope
JuliaReach · Jun 25, 2020 · 18a53ae · 18a53ae
2 parents d97d1db + 5288647
commit 18a53ae
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 0 deletions.
diff --git a/src/Flowpipes/setops.jl b/src/Flowpipes/setops.jl
@@ -453,6 +453,12 @@ end
     return c, r
 end
 
+# overapproximate a hyperrectangular set with a polytope
+function _overapproximate(H::AbstractHyperrectangle, ::Type{<:HPolytope})
+    P = overapproximate(H, Hyperrectangle)
+    HPolytope([HalfSpace(Vector(c.a), c.b) for c in constraints_list(H)])
+end
+
 # ==================================
 # Zonotope splitting methods
 # ==================================

diff --git a/src/Hybrid/transitions.jl b/src/Hybrid/transitions.jl
@@ -74,6 +74,9 @@ _affine_term(m::MathematicalSystems.AffineMap{N, <:LazySet}) where {N} = m.c
 _affine_term(m::ConstrainedAffineMap{N, MT, VT, <:AbstractVector}) where {N, MT, VT} = Singleton(m.c)
 _affine_term(m::ConstrainedAffineMap{N, MT, VT, <:LazySet}) where {N, MT, VT} = Singleton(m.c)
 
+# TODO see LazySets#2198
+LazySets.dim(c::ConstrainedResetMap) = c.dim
+
 function _state_matrix(m::Union{<:MathematicalSystems.ResetMap{N},
                                 <:MathematicalSystems.ConstrainedResetMap{N}}) where {N}
     n = dim(m)