Io

homogeniseCouplingConstraints(model)[source]

Homogenise the coupling constraints of a model so that every retained constraint couples exactly a pair of reactions

Coupling constraints with a single non-zero entry are removed, constraints with three non-zero entries are reduced by removing the coupling associated with reaction sIEC_biomass_reactionIEC01b_trtr, and constraints whose two remaining non-zero entries are both negative are sign-flipped to be both positive (updating the constraint sense accordingly)

USAGE:

model = homogeniseCouplingConstraints (model)

INPUTS:

model – COBRA model structure with coupling-constraint fields:

  • .C - ctrs x n coupling-constraint matrix

  • .d - ctrs x 1 right-hand side values for the coupling constraints

  • .ctrs - ctrs x 1 coupling-constraint identifiers

  • .dsense - ctrs x 1 coupling-constraint senses (E, L, G)

  • .rxns - n x 1 reaction identifiers

OUTPUTS:

model – COBRA model with homogenised coupling constraints, updating the fields:

  • .C - coupling-constraint matrix with single-entry and triple-entry rows resolved and negative pairs sign-flipped

  • .d - right-hand side values for the retained constraints

  • .ctrs - identifiers of the retained coupling constraints

  • .dsense - senses of the retained coupling constraints

Note

Assumes sIEC_biomass_reactionIEC01b_trtr and sIEC_biomass_reactionIEC01b are the only pair of doubly coupled reactions and will error if not