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
|
function y = Parameter_isInAlgebraicRule(SBMLParameter, SBMLRules)
% y = Parameter_isInAlgebraicRule(SBMLParameter, SBMLRules)
%
% Takes
%
% 1. SBMLParameter, an SBML Parameter structure
% 2. SBMLRules; the array of rules from an SBML Model structure
%
% Returns
%
% y =
% - an array of the indices of any algebraicRules the id of the Parameter appears in
% - 0 if the Parameter appears in no algebraicRules
%
%<!---------------------------------------------------------------------------
% This file is part of SBMLToolbox. Please visit http://sbml.org for more
% information about SBML, and the latest version of SBMLToolbox.
%
% Copyright (C) 2009-2012 jointly by the following organizations:
% 1. California Institute of Technology, Pasadena, CA, USA
% 2. EMBL European Bioinformatics Institute (EBML-EBI), Hinxton, UK
%
% Copyright (C) 2006-2008 jointly by the following organizations:
% 1. California Institute of Technology, Pasadena, CA, USA
% 2. University of Hertfordshire, Hatfield, UK
%
% Copyright (C) 2003-2005 jointly by the following organizations:
% 1. California Institute of Technology, Pasadena, CA, USA
% 2. Japan Science and Technology Agency, Japan
% 3. University of Hertfordshire, Hatfield, UK
%
% SBMLToolbox 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.
%----------------------------------------------------------------------- -->
y = 0;
%-------------------------------------------------------------------
% check input arguments are as expected
if (~isstruct(SBMLParameter))
error(sprintf('%s', ...
'argument must be an SBML Parameter structure'));
end;
[sbmlLevel, sbmlVersion] = GetLevelVersion(SBMLParameter);
if (~isSBML_Parameter(SBMLParameter, sbmlLevel, sbmlVersion))
error('Parameter_isInAlgebraicRule(SBMLParameter, SBMLRules)\n%s', ...
'first argument must be an SBMLParameter structure');
end;
NumRules = length(SBMLRules);
if (NumRules < 1)
error('Parameter_isInAlgebraicRule(SBMLParameter, SBMLRules)\n%s', ...
'SBMLRule structure is empty');
else
for i = 1:NumRules
if (~isSBML_Rule(SBMLRules(i), sbmlLevel, sbmlVersion))
error('Parameter_isInAlgebraicRule(SBMLParameter, SBMLRules)\n%s', ...
'second argument must be an array of SBMLRule structures');
end;
end;
end;
%--------------------------------------------------------------------------
% loop through each rule and check whether the Parameter occurs
%determine the name or id of the Parameter
if (sbmlLevel == 1)
name = SBMLParameter.name;
else
if (isempty(SBMLParameter.id))
name = SBMLParameter.name;
else
name = SBMLParameter.id;
end;
end;
y = [];
for i = 1:NumRules
index = matchName(SBMLRules(i).formula, name);
if (~isempty(index) && strcmp(SBMLRules(i).typecode, 'SBML_ALGEBRAIC_RULE'))
y = [y;i];
end;
end;
if isempty(y)
y = 0;
end;
|
