work around equality bounds in coupling

(as produced by mgPipe)
LCSB-BioCore · May 18, 2024 · b9eceed · b9eceed
1 parent 50f6ffa
commit b9eceed
Showing 1 changed file with 5 additions and 0 deletions.
diff --git a/src/base/types/MATModel.jl b/src/base/types/MATModel.jl
@@ -108,6 +108,11 @@ function coupling_bounds(m::MATModel)
             sparse(fill(-Inf, nc)),
             sparse(reshape(m.mat["b"], length(m.mat["b"]))[n_reactions(m)+1:end]),
         )
+    elseif haskey(m.mat, "d")
+        (
+            sparse(reshape(get(m.mat, "d", fill(-Inf, nc, 1)), nc)),
+            sparse(reshape(get(m.mat, "d", fill(Inf, nc, 1)), nc)),
+        )
     else
         (
             sparse(reshape(get(m.mat, "cl", fill(-Inf, nc, 1)), nc)),