def simplex_rule(m, l):
    return (pe.quicksum((m.q[(i, l)] for (i, ll) in m.q if (ll == l))) == 1)
def simplex_rule(m, l):
    return (pe.quicksum((m.q[(i, l)] for (i, ll) in m.q if (ll == l))) == 1)