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
|
function SBMLReaction = Reaction_addReactant(SBMLReaction, SBMLReactant)
% SBMLReaction = Reaction_addReactant(SBMLReaction, SBMLReactant)
%
% Takes
%
% 1. SBMLReaction, an SBML Reaction structure
% 2. SBMLReactant, an SBML Reactant structure
%
% Returns
%
% 1. the SBML Reaction structure with the SBML Reactant structure added
%
%<!---------------------------------------------------------------------------
% 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.
%----------------------------------------------------------------------- -->
%get level and version and check the input arguments are appropriate
[level, version] = GetLevelVersion(SBMLReaction);
[reactant_level, reactant_version] = GetLevelVersion(SBMLReactant);
if level ~= reactant_level
error('mismatch in levels');
elseif version ~= reactant_version
error('mismatch in versions');
end;
if isfield(SBMLReaction, 'reactant')
index = length(SBMLReaction.reactant);
if index == 0
SBMLReaction.reactant = SBMLReactant;
else
SBMLReaction.reactant(index+1) = SBMLReactant;
end;
else
error('reactant not an element on SBML L%dV%d Reaction', level, version);
end;
|
