yeild objective

[Daniel Farkas]

Dear cobrapy community,

I am trying to set the objective function to be the ATP production yield. The below code does not work, because it is not a linear function. Do you know how to get around this problem?


ATPyield_objective = model.problem.Objective(
    model.reactions.ATPS4rpp.flux_expression/model.reactions.EX_glc__D_e.flux_expression,
    direction='max')
model.objective = ATPyield_objective
solution = model.optimize()

Thank you kindly,

Daniel