Fix bug causing matrix dimension mismatch in state expansion when cus…

…tom constraints are used.
MunskyGroup · Oct 22, 2024 · d71de53 · d71de53
1 parent 776bd1a
commit d71de53
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/+ssit/@FspMatrixTerm/FspMatrixTerm.m b/src/+ssit/@FspMatrixTerm/FspMatrixTerm.m
@@ -304,9 +304,9 @@
             k = 1;
             for c = 1:state_set.numConstraints
                 ninsert = nnz(state_set.outboundTransitions(:, c + state_set.numConstraints*(ireaction-1)));
-                isinks(k:k+ninsert-1) = n_states + c;
-                jsinks(k:k+ninsert-1) = find(state_set.outboundTransitions(:,c + state_set.numConstraints*(ireaction-1)));
-                aijsinks(k:k+ninsert-1) = prop_val(jsinks(k:k+ninsert-1))./n_constrs_failed(jsinks(k:k+ninsert-1));
+                isinks(k:k+ninsert-1,1) = n_states + c;
+                jsinks(k:k+ninsert-1,1) = find(state_set.outboundTransitions(:,c + state_set.numConstraints*(ireaction-1)));
+                aijsinks(k:k+ninsert-1,1) = prop_val(jsinks(k:k+ninsert-1))./n_constrs_failed(jsinks(k:k+ninsert-1));
 
                 k = k + ninsert;
             end