def schedule_all(self, timelimit=None):
if not self.reservation_list:
return self.schedule_dict
self.build_data_structures()
# allocate A & b
# find the row size of A:
# first find the number of reservations participating in oneofs
oneof_reservation_num = 0
for c in self.oneof_constraints:
oneof_reservation_num += len(c)
A_numrows = len(self.reservation_list) + len(self.aikt) + len(
self.oneof_constraints) - oneof_reservation_num
A_rows = []
A_cols = []
A_data = []
# try:
# A = numpy.zeros((A_numrows, len(self.Yik)), dtype=numpy.int)
# except ValueError:
# print("Number of A rows: {}".format(A_numrows))
b = numpy.zeros(A_numrows, dtype=numpy.int16)
# build A & b
row = 0
# constraint 5: oneof
for c in self.oneof_constraints:
for r in c:
for entry in r.Yik_entries:
A_rows.append(row)
A_cols.append(entry)
A_data.append(1)
# A[row,entry] = 1
r.skip_constraint2 = True
b[row] = 1
row += 1
# constraint 2: each res should have one start:
# optimization:
# if the reservation participates in a oneof, then this is
# redundant with the oneof constraint added above, so don't add it.
for r in self.reservation_list:
if hasattr(r, 'skip_constraint2'):
continue
for entry in r.Yik_entries:
A_rows.append(row)
A_cols.append(entry)
A_data.append(1)
# A[row,entry] = 1
b[row] = 1
row += 1
# constraint 3: each slice should only have one sched. reservation:
for s in self.aikt.keys():
for entry in self.aikt[s]:
A_rows.append(row)
A_cols.append(entry)
A_data.append(1)
# A[row,entry] = 1
b[row] = 1
row += 1
A = coo_matrix((A_data, (A_rows, A_cols)),
shape=(A_numrows, len(self.Yik)))
# constraint 6: and
# figure out size of constraint matrix
if not self.and_constraints:
Aeq = []
beq = []
else:
Aeq_numrows = 0
for c in self.and_constraints:
Aeq_numrows += len(c) - 1
# allocate Aeq and beq
# Aeq = numpy.zeros((Aeq_numrows, len(self.Yik)), dtype=numpy.int)
Aeq_rows = []
Aeq_cols = []
Aeq_data = []
beq = numpy.zeros(Aeq_numrows, dtype=numpy.int16)
row = 0
for c in self.and_constraints:
constraint_size = len(c)
left_idx = 0
right_idx = 1
while right_idx < constraint_size:
left_r = c[left_idx]
right_r = c[right_idx]
for entry in left_r.Yik_entries:
# Aeq[row, entry] = 1
Aeq_rows.append(row)
Aeq_cols.append(entry)
Aeq_data.append(1)
for entry in right_r.Yik_entries:
Aeq_rows.append(row)
Aeq_cols.append(entry)
Aeq_data.append(-1)
# Aeq[row, entry] = -1
left_idx += 1
right_idx += 1
row += 1
# print(Aeq_numrows)
Aeq = coo_matrix((Aeq_data, (Aeq_rows, Aeq_cols)),
shape=(Aeq_numrows, len(self.Yik)))
# bounds:
lb = numpy.zeros(len(self.Yik), dtype=numpy.int16)
ub = numpy.ones(len(self.Yik), dtype=numpy.int16)
# objective function:
f = numpy.zeros(len(self.Yik))
f = numpy.zeros(len(self.Yik), dtype=numpy.int16)
row = 0
for entry in self.Yik:
f[row] = entry[2] # priority
row += 1
dump_matrix_sizes(f, A, Aeq, b, beq, lb, ub,
len(self.compound_reservation_list))
p = LP(f=f, A=A, Aeq=Aeq, b=b, beq=beq, lb=lb, ub=ub)
# r = p.minimize('pclp')
r = p.maximize('glpk', iprint=-1)
# r = p.minimize('lpsolve')
return self.unpack_result(r)
