From: Computation of elementary modes: a unifying framework and the new binary approach
N ← reconfigured stoichiometry matrix [N - N Rev ] A ← [NT-NTI]T | Reconfiguration |
---|---|
A q ← q independent rows of A R ← A q -1 | Initialization |
for each unprocessed row Ai•of A do J+ ← {j ∈ J : Ai•rj> 0} J0 ← {j ∈ J : Ai•rj= 0} J- ← {j ∈ J : Ai•rj< 0} R' ← {rj: j ∈ J0 ∪ J+} For (j+, j-) ∈ J+ × J- do | Processing of constraints in a given order |
If and adjacent in R then | Adjacency test |
end if end for R ← R' End for | Gaussian combination step |
R ← R \ { futile two-cycles } R ← back-configuration of R | Back-configuration |