network_check.m

％The Metabolic Networks Toolbox contains functions to create, ％modify, display, and simulate bioche

function network = network_check(network)

%network = network_check(network)
%
%prints an error message if dimensions of network fields 
%do not match the number of reactions and metabolites.

n_S = length(network.metabolites);
n_A = length(network.actions);

if max(size(network.N) ~= [n_S,n_A]),        fprintf('network_check: error\n'); end
if max(size(network.reversible) ~= [n_A,1]), fprintf('network_check: error\n'); end

if isfield(network,'kinetics'),
  if max(size(network.kinetics.k_fwd) ~= [n_A,1]), fprintf('network_check: error\n'); end
  if max(size(network.kinetics.k_bwd) ~= [n_A,1]), fprintf('network_check: error\n'); end
end