network_choose.m

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

%network=network_choose(network,indm,indr)
% 
%choose subnetwork containing metabolites/reactions 
% with the (old) indices indm (for metabolites) and 
% indr (for reactions)

function network = network_choose(network,indm,indr)

 dummi               = zeros(1,length(network.metabolites));
 network.metabolites = network.metabolites(indm);
 network.N           = network.N(indm,indr);
 network.reversible  = network.reversible(indr);
 network.actions     = network.actions(indr);
 network.EC          = network.EC(indr);
 network.formulae    = network.formulae(indr);
 if isfield(network,'index'), network.index=network.index(indr); end
  dummi(find(network.external))=1;
 network.external   = bit_vector(find(dummi(indm)),length(network.metabolites));
 if isfield(network,'graphics_par'), network=netgraph_make_graph(network); end
 
 if isfield(network,'kinetics'),
   fprintf('Warning: updating the kinetics is not yet supported.\n');
 end