fix: readYAMLmodel correct rxnEnzMat

SysBioChalmers · Oct 17, 2024 · 58f4ef2 · 58f4ef2
1 parent 0e45b13
commit 58f4ef2
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/doc/io/readYAMLmodel.html b/doc/io/readYAMLmodel.html
@@ -682,7 +682,7 @@ <h2><a name="_source"></a>SOURCE CODE <a href="#_top"><img alt="^" border="0" sr
 0625     [~,enzIdx]          = ismember(enzStoich(:,2),model.ec.enzymes);
 0626     coeffs              = cell2mat(enzStoich(:,3));
 0627     model.ec.rxnEnzMat  = zeros(numel(model.ec.rxns), numel(model.ec.genes));
-0628     linearIndices       = sub2ind([max(rxnIdx), max(enzIdx)], rxnIdx, enzIdx);
+0628     linearIndices       = sub2ind([numel(model.ec.rxns), numel(model.ec.genes)], rxnIdx, enzIdx);
 0629     model.ec.rxnEnzMat(linearIndices) = coeffs;
 0630     <span class="comment">%Parse ec-codes</span>
 0631     <span class="keyword">if</span> ~isempty(ecGecko)

diff --git a/io/readYAMLmodel.m b/io/readYAMLmodel.m
@@ -625,7 +625,7 @@
     [~,enzIdx]          = ismember(enzStoich(:,2),model.ec.enzymes);
     coeffs              = cell2mat(enzStoich(:,3));
     model.ec.rxnEnzMat  = zeros(numel(model.ec.rxns), numel(model.ec.genes));
-    linearIndices       = sub2ind([max(rxnIdx), max(enzIdx)], rxnIdx, enzIdx);
+    linearIndices       = sub2ind([numel(model.ec.rxns), numel(model.ec.genes)], rxnIdx, enzIdx);
     model.ec.rxnEnzMat(linearIndices) = coeffs;
     %Parse ec-codes
     if ~isempty(ecGecko)