added constraints_list for CartesianProduct

src/CartesianProduct.jl

@@ -151,6 +151,23 @@ function ∈(x::AbstractVector{<:Real}, cp::CartesianProduct)::Bool
            ∈(view(x, n1+1:length(x)), cp.Y)
 end
 
+"""
+    constraints_list(cp::CartesianProduct{N})::Vector{LinearConstraint{N}} where N<:Real
+
+Return the list of constraints of a (polytopic) Cartesian product.
+
+### Input
+
+- `cp` -- Cartesian product
+
+### Output
+
+A list of constraints.
+"""
+function constraints_list(cp::CartesianProduct{N})::Vector{LinearConstraint{N}} where N<:Real
+    return constraints_list(CartesianProductArray([cp.X, cp.Y]))
+end
+
 """
     vertices_list(cp::CartesianProduct{N})::Vector{Vector{N}} where N<:Real
 
@@ -353,6 +370,53 @@ function ∈(x::AbstractVector{N}, cpa::CartesianProductArray{N, <:LazySet{N}}
     return true
 end
 
+"""
+    constraints_list(cpa::CartesianProductArray{N})::Vector{LinearConstraint{N}} where N<:Real
+
+Return the list of constraints of a (polytopic) Cartesian product.
+
+### Input
+
+- `cpa` -- Cartesian product
+
+### Output
+
+A list of constraints.
+
+"""
+function constraints_list(cpa::CartesianProductArray{N})::Vector{LinearConstraint{N}} where N<:Real
+    # collect low-dimensional constraints lists
+    clist_low = []
+    n = 0
+    for c_low in array(cpa)
+        if c_low isa LinearConstraint
+            push!(clist_low, [c_low])
+            n += 1
+        else
+            constraints = constraints_list(c_low)
+            push!(clist_low, constraints)
+            n += length(constraints)
+        end
+    end
+
+    clist = Vector{LinearConstraint{N}}()
+    sizehint!(clist, n)
+    prev_step = 1
+    # create high-dimensional constraints list
+    for c_low in clist_low
+        if !isempty(c_low)
+            indices = prev_step : (dim(c_low[1]) + prev_step - 1)
+        end
+        for constr in c_low
+            new_constr = LinearConstraint(sparsevec(prev_step : (dim(constr) + prev_step-1), constr.a), constr.b)
+            push!(clist, new_constr)
+        end
+        prev_step += dim(c_low[1])
+    end
+
+    return clist
+end
+
 """
     vertices_list(cpa::CartesianProductArray{N})::Vector{Vector{N}} where N<:Real