Add fix for empty matrices

laschuet · Dec 22, 2020 · 40aedf2 · 40aedf2
1 parent df6b5d2
commit 40aedf2
Showing 1 changed file with 8 additions and 0 deletions.
diff --git a/src/matrix.jl b/src/matrix.jl
@@ -65,6 +65,14 @@ respectively) of `A` or `B`.
 """
 function diff(A::AbstractMatrix, B::AbstractMatrix, ia::AbstractVector,
             ja::AbstractVector, ib::AbstractVector, jb::AbstractVector)
+    if size(A) == (0, 0) || size(B) == (0, 0)
+        T = promote_type(eltype(A), eltype(B))
+        modvals = sparse([], [], T[])
+        addvals = view(Vector{T}(undef, 0), :)
+        remvals = view(Vector{T}(undef, 0), :)
+        return MatrixDifference(modvals, addvals, remvals)
+    end
+
     iamap = Dict(zip(ia, 1:length(ia)))
     jamap = Dict(zip(ja, 1:length(ja)))
     ibmap = Dict(zip(ib, 1:length(ib)))