master_constraints.py

def masterConstraints(M, d, x, y, m, n, t, q, Q, f, gb):
    '''constraint 2'''
    for j in range(n):
        M.addConstr(gb.quicksum(x[i,j] for i in range(m)) ==  d[j])
    M.update()
    '''constraint 2'''
    for j in range(n):
        M.addConstr(gb.quicksum(gb.quicksum( y[i, j + 1, k] for k in range(t)) for i in range(m)) == f[j+1] * d[j])
    M.update()
    '''constraint 3'''
    for i in range(m):
        for j in range(n):
            for k in range(t):
                M.addConstr(y[i, j + 1, k] <= x[i, j])
    M.update()
    '''constraint 4'''
    for i in range(m):
        for j in range(n):
            if q[j+1] != Q[i]:
                M.addConstr(x[i, j] == 0)
    M.update()

    '''constraint 5'''
    for i in range(m):
        for k in range(t):
            M.addConstr(y[i, 0, k] == 1)  # j = 0 for depot location
            M.addConstr(y[i, 0, k] == 1)
    M.update()

    '''constraint 6'''
    for i in range(m):
        for j in range(n):
            for k in range(t):
                for tou in range(1, (4 - (f[j+1]) + 1)):
                    if (k + tou) <= 4:
                        M.addConstr(y[i, j + 1, k] + y[i, j + 1, k + tou] <= 1)
    M.update()