/usr/share/SBMLToolbox/AccessModel/GetVaryingParameters.m is in sbmltoolbox 4.1.0-2.
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 | function varargout = GetVaryingParameters(SBMLModel)
% [names, values] = GetVaryingParameters(SBMLModel)
%
% Takes
%
% 1. SBMLModel, an SBML Model structure
%
% Returns
%
% 1. an array of strings representing the identifiers of any non-constant parameters
% within the model
% 2. an array of the values of each of these parameter
%
% *NOTE:* the value returned will be (in order)
%
% - determined from assignmentRules/initialAssignments where appropriate
% - the attribute 'value' for the given parameter
% - NaN; if the value is not specified in any way within the model
%<!---------------------------------------------------------------------------
% 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(SBMLModel))
error('GetVaryingParameters(SBMLModel)\n%s', 'input must be an SBMLModel structure');
end;
%------------------------------------------------------------
% determine the number of parameters within the model
NumParams = length(SBMLModel.parameter);
count = 1;
%------------------------------------------------------------
% loop through the list of parameters
for i = 1:NumParams
%determine the name or id of the parameter
if (SBMLModel.SBML_level == 1)
name = SBMLModel.parameter(i).name;
else
if (isempty(SBMLModel.parameter(i).id))
name = SBMLModel.parameter(i).name;
else
name = SBMLModel.parameter(i).id;
end;
end;
% if the parameter is not constant add to arrays
if SBMLModel.parameter(i).constant == 0
% save into an array of character names
CharArray{count} = name;
% put the value into the array
Values(count) = SBMLModel.parameter(i).value;
% might be set by assignment rule
AR = Model_getAssignmentRuleByVariable(SBMLModel, name);
if (~isempty(AR))
newSBMLModel = SBMLModel;
newSBMLModel.parameter(i) = [];
for fd = 1:Model_getNumFunctionDefinitions(SBMLModel)
newFormula = SubstituteFunction(AR.formula, Model_getFunctionDefinition(SBMLModel, fd));
if (~isempty(newFormula))
AR.formula = newFormula;
end;
end;
Values(count) = Substitute(AR.formula, newSBMLModel);
end;
% might be an initial assignment in l2v2
if ((SBMLModel.SBML_level == 2 && SBMLModel.SBML_version > 1) ...
|| SBMLModel.SBML_level == 3)
IA = Model_getInitialAssignmentBySymbol(SBMLModel, name);
if (~isempty(IA))
% remove this from the substtution
newSBMLModel = SBMLModel;
newSBMLModel.parameter(i) = [];
for fd = 1:Model_getNumFunctionDefinitions(SBMLModel)
newFormula = SubstituteFunction(IA.math, Model_getFunctionDefinition(SBMLModel, fd));
if (~isempty(newFormula))
IA.math = newFormula;
end;
end;
Values(count) = Substitute(IA.math, newSBMLModel);
end;
end;
count = count + 1;
end;
end;
%--------------------------------------------------------------------------
% assign output
if (count ~= 1)
varargout{1} = CharArray;
varargout{2} = Values;
else
varargout{1} = [];
varargout{2} = [];
end;
|