add constructor of matrix power given a target index

JuliaReach · Feb 11, 2020 · d67d1c0 · d67d1c0
1 parent 16edd54
commit d67d1c0
Showing 1 changed file with 9 additions and 0 deletions.
diff --git a/src/power.jl b/src/power.jl
@@ -65,6 +65,15 @@ mutable struct IntervalMatrixPower{T}
         @assert k >= 1 "matrix powers must be positive"
         return new{T}(M, Mᵏ, k)
     end
+
+    function IntervalMatrixPower(M::IntervalMatrix{T}, k::Int) where {T}
+        @assert k >= 1 "matrix powers must be positive"
+        pow = IntervalMatrixPower(M, M, 1)
+        @inbounds for i in 1:k
+            increment!(pow)
+        end
+        return pow
+    end
 end
 
 function copy(pow::IntervalMatrixPower)