ES_FIX_CONSTRAINT_VIOLATIONS Change metabolic state variables to enforce a feasible state [c,dmu] = es_fix_constraint_violations(v,c,dmu,constraints,N) 1. Threshold small and large concentrations 2. Recalculate reaction affinities accordingly 3. Threshold the resulting reaction affinities, giving up their proper relationship to concentrations