added constructor for complex interval matrices

src/matrix.jl

@@ -1,4 +1,4 @@
-import Base: similar, ∈, ⊆, ∩, ∪
+import Base: similar, ∈, ⊆, ∩, ∪, real, imag
 import Random: rand
 import IntervalArithmetic: ±, inf, sup, mid, diam, radius, midpoint_radius, hull
 
@@ -83,6 +83,10 @@ function IntervalMatrix(A::MT) where {T, IT<:AbstractInterval{T}, MT<:AbstractMa
     return IntervalMatrix{T, IT, MT}(A)
 end
 
+function IntervalMatrix(A::MT) where {T, IT<:Complex{Interval{T}}, MT <:AbstractMatrix{IT}}
+    return IntervalMatrix{T, IT, MT}(A)
+end
+
 # constructor from uniform scaling
 function IntervalMatrix(αI::UniformScaling{Interval{T}}, m::Integer, n::Integer=m) where {T}
     return IntervalMatrix(Matrix(αI, m, n))

src/operations/numops.jl

@@ -105,3 +105,7 @@ A pair `(C, S)` such that the entries of `C` are the central points and the
 entries of `S` are the (nonnegative) radii of the intervals in `A`.
 """
 midpoint_radius(A::IntervalMatrix) = (mid(A), radius(A))
+
+
+real(M::IntervalMatrix) = IntervalMatrix(real(M.mat))
+imag(M::IntervalMatrix) = IntervalMatrix(imag(M.mat))

test/constructors.jl

@@ -22,3 +22,21 @@ end
     @test c ≈ [-0.1 -4.; 4. -0.1]
     @test s ≈ [1. 0.1; 0.1 1.]
 end
+
+@testset "Complex interval matrices" begin
+    m1 = IntervalMatrix([(1..2)+im*(3..4) 1;2 3])
+    @test m1 isa IntervalMatrix{Float64}
+    @test eltype(m1) == Complex{Interval{Float64}}
+
+    rp = real(m1)
+    ip = imag(m1)
+    @test rp isa IntervalMatrix{Float64}
+    @test eltype(rp) == Interval{Float64}
+
+    @test ip isa IntervalMatrix{Float64}
+    @test eltype(ip) == Interval{Float64}
+
+    m2 = IntervalMatrix([1+im 2+im;3+im 4+im])
+    @test m2 isa IntervalMatrix{Float64}
+    @test eltype(m2) == Complex{Interval{Float64}}
+end