#2453 - Vertices list of IntervalArithmetic types (#2457)

* vertices list of IA types * typo * Update unit_Interval.jl
JuliaReach · Dec 22, 2020 · 3910423 · 3910423
1 parent 7a12227
commit 3910423
Show file tree

Hide file tree

Showing 2 changed files with 24 additions and 0 deletions.
diff --git a/src/Sets/Interval.jl b/src/Sets/Interval.jl
@@ -699,3 +699,19 @@ end
         return result
     end
 end
+
+function vertices_list(x::IntervalArithmetic.Interval{N}) where {N}
+    a = IntervalArithmetic.inf(x)
+    b = IntervalArithmetic.sup(x)
+    ST = IntervalArithmetic.SVector{1, N}
+    if _isapprox(a, b)
+        vlist = [ST(a)]
+    else
+        vlist = [ST(a), ST(b)]
+    end
+    return vlist
+end
+
+function vertices_list(H::IntervalArithmetic.IntervalBox)
+    return vertices_list(convert(Hyperrectangle, H))
+end
diff --git a/test/unit_Interval.jl b/test/unit_Interval.jl
@@ -208,4 +208,12 @@ for N in Ns
     @test rectify(x) == Interval(N(0), N(2))
     x = Interval(N(1), N(2))
     @test rectify(x) == x
+
+    # list of vertices of IA types
+    b = IA.IntervalBox(IA.Interval(0 , 1), IA.Interval(0, 1))
+    vlistIB = vertices_list(b)
+    @test is_cyclic_permutation(vlistIB, [SA[N(1), N(1)], SA[N(0), N(1)], SA[N(1), N(0)], SA[N(0), N(0)]])
+
+    vlistI = vertices_list(b[1])
+    @test is_cyclic_permutation(vlistI, [SA[N(0)], SA[N(1)]])
 end