Properly increment/decrement during pattern check

JuliaSparse · Dec 14, 2021 · 65db5e1 · 65db5e1 · rayegun · Dec 14, 2021
1 parent 2dbb52d
commit 65db5e1
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/KLU.jl b/src/KLU.jl
@@ -560,7 +560,11 @@ end
 
 function klu!(K::KLUFactorization{U}, S::SparseMatrixCSC{U}) where {U}
  size(K) == size(S) || throw(ArgumentError("Sizes of K and S must match."))
+ increment!(K.colptr)
+ increment!(K.rowval)
  K.colptr == S.colptr && K.rowval == S.rowval || throw(ArgumentError("The pattern of the original matrix must match the pattern of the refactor."))
+ decrement!(K.colptr)
+ decrement!(K.rowval)
  return klu!(K, S.nzval)
 end
 #B is the modified argument here. To match with the math it should be (klu, B). But convention is