def checkBlocked(model_path, reacs=None):
lp = metabolism.Metabolism(util.ImportCplex(model_path))
lp.simplex()
print lp
if reacs:
result = analyseBlockedReactions(lp, reacs)
else:
result = analyseBlockedReactions(lp)
return result
def determineFreeReactions(lp, setCarbonSources, freeReacsSet=None):
if not freeReacsSet:
allReactions = lp.getReactions()
freeReacsSet = set(allReactions)
for carbSource in setCarbonSources:
print 'length of freeReactionSet', len(freeReacsSet)
print carbSource
lp.modifyColumnBounds({carbSource:(0,20)})
print 'normal FBA: ', lp.fba()[biom]
blockedReactions = set(bc.analyseBlockedReactions(lp, freeReacsSet))
lp.initialize()
freeReacsSet.difference_update(blockedReactions)
print 5*'\n'
return freeReacsSet
