/usr/share/SBMLToolbox/AccessModel/GetStoichiometrySparse.m is in sbmltoolbox 4.1.0-3.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 | function S=GetStoichiometrySparse(model)
% S = GetStoichiometrySparse(SBMLModel)
%
% Takes
%
% 1. SBMLModel, an SBML Model structure
%
% Returns
%
% 1. a sparse stoichiometry matrix produced from the reactions/species
%
% *NOTE:* This function was contributed by: Arsen Batagov (2006)
%<!---------------------------------------------------------------------------
% This file is part of SBMLToolbox. Please visit http://sbml.org for more
% information about SBML, and the latest version of SBMLToolbox.
%
% Copyright 2005-2007 California Institute of Technology.
% Copyright 2002-2005 California Institute of Technology and
% Japan Science and Technology Corporation.
%
% This library is free software; you can redistribute it and/or modify it
% under the terms of the GNU Lesser General Public License as published by
% the Free Software Foundation. A copy of the license agreement is provided
% in the file named "LICENSE.txt" included with this software distribution.
% and also available online as http://sbml.org/software/sbmltoolbox/license.html
%----------------------------------------------------------------------- -->
% check input is an SBML model
if (~isValidSBML_Model(model))
error('GetStoichiometrySparse(model)\n%s', 'input must be an SBMLModel structure');
end;
if (model.SBML_level > 1)
SpeciesID={model.species.id};
ReactionID={model.reaction.id};
else
SpeciesID={model.species.name};
ReactionID={model.reaction.name};
end;
NS=length(SpeciesID);
NR=length(ReactionID);
LNS=zeros(NS,1);
S=sparse(zeros(NS,NR));
for i=1:NS
LNS(i)=length(SpeciesID{i});
end
% Create a hash-table for the future search
MaxLNS=max(LNS);
InitialApprox=cell(MaxLNS);
LIA=zeros(MaxLNS);
for i=min(LNS):MaxLNS
InitialApprox{i}=find(LNS==i);
LIA(i)=length(InitialApprox{i});
end
% Search
for i=1:NR
t={model.reaction(i).reactant.species};
if ((model.SBML_level == 2) && (~isempty(model.reaction(i).reactant.stoichiometryMath)))
error('GetStoichiometrySparse(SBMLModel)\n%s', 'stoichiometry has been entered as a formula');
end;
if ((model.SBML_level == 2 && model.SBML_version > 1) ...
|| model.SBML_level == 3)
denominator = 1.0;
else
denominator = [model.reaction(i).reactant.denominator];
end;
s=[model.reaction(i).reactant.stoichiometry]./double(denominator);
rl=length(t);
if ~isempty(s)
if rl>1
for j=1:rl
lt=length(t{j});
for k=1:LIA(lt)
if cell2mat(SpeciesID(InitialApprox{lt}(k)))==t{j}
S(InitialApprox{lt}(k),i)=-s(j);
end
end
end
else
t=cell2mat(t);
lt=length(t);
for k=1:LIA(lt)
if cell2mat(SpeciesID(InitialApprox{lt}(k)))==t
S(InitialApprox{lt}(k),i)=-s;
end
end
end
end
t={model.reaction(i).product.species};
if ((model.SBML_level == 2) && (~isempty(model.reaction(i).product.stoichiometryMath)))
error('GetStoichiometrySparse(SBMLModel)\n%s', 'stoichiometry has been entered as a formula');
end;
if ((model.SBML_level == 2 && model.SBML_version > 1) ...
|| model.SBML_level == 3)
denominator = 1.0;
else
denominator = [model.reaction(i).product.denominator];
end;
s=[model.reaction(i).product.stoichiometry]./double(denominator);
rl=length(t);
if ~isempty(s)
if rl>1
for j=1:rl
lt=length(t{j});
for k=1:LIA(lt)
if cell2mat(SpeciesID(InitialApprox{lt}(k)))==t{j}
S(InitialApprox{lt}(k),i)=s(j);
end
end
end
else
t=cell2mat(t);
lt=length(t);
for k=1:LIA(lt)
if cell2mat(SpeciesID(InitialApprox{lt}(k)))==t
S(InitialApprox{lt}(k),i)=s;
end
end
end
end
end
|