If you define f over q as a dynamic set (and not an alias), I think you can do it as follows:
Code: Select all
set f(q) dynamic set;
...
additionalrestriction.. sum(f, K(f))=e=0;
...
loop (q,
solve Model9 MAXIMIZING z using MIP;
if (z.l gt 0,
f(q)$k.l(q) = yes;
);
);
...
Renger