def test_rateLoss_monotonicity(self):
input = get_general_input()
subset_amount_of_routes = input["subset_amount_of_routes"]
sum_r = sum(input["arrival_rates"])
C = input["capacity"]
A = input["requirements_of_routes"]
arrRates = input["arrival_rates"]
serRates = input["service_rates"]
block_prob = bp_v1(sum_r, C, A)
lossRate = rl_v1(block_prob, arrRates, serRates)
for i in range(subset_amount_of_routes):
A_new = A[0:i]
A_new += A[i + 1:]
arrRates_new = arrRates[0:i]
arrRates_new += arrRates[i + 1:]
serRates_new = serRates[0:i]
serRates_new += serRates[i + 1:]
block_prob_new = bp_v1(sum(arrRates_new), C, A_new)
lossRate_new = rl_v1(block_prob_new, arrRates_new, serRates_new)
message = "By removing route " + str(
i) + " the rate of loss becomes " + str(lossRate_new)
message += ". Its original value was " + str(lossRate)
self.assertTrue(bool(lossRate_new <= lossRate), str(message))
def test_simple_cases(self):
p1 = ["capacity", "int", [1.0], [10.0], 1]
p2 = ["amount_of_routes", "int", [3.0], [8.0], 1]
p3 = ["arrival_rates", "float", [0.0], [0.0], "amount_of_routes"]
p4 = ["service_rates", "float", [1.0], [1.0], "amount_of_routes"]
p5 = [
"subset_amount_of_routes", "int", [2.0],
["amount_of_routes", "-", 1], 1
]
p6 = [
"requirements_of_routes", "int", [1.0], ["capacity", "/", 2],
"amount_of_routes"
]
input = generate_input([p1, p2, p3, p4, p5, p6])
bp = bp_v1(sum(input['arrival_rates']), input['capacity'],
input['requirements_of_routes'])
for i in range(input['amount_of_routes']):
assert (not input['arrival_rates'][i]
== 0) or bp[i] == 0, "arrival rate is " + str(
input['arrival_rates'][i]) + "blocking prob is " + str(
bp[i])
assert (not input['requirements_of_routes'][i]
== 0) or bp[i] == 0, "requirement is " + str(
input['requirements_of_routes']
[i]) + "blocking prob is " + str(bp[i])
def test_blockProbs_small_vs_large_requirements(self):
input = get_same_rate_input()
sum_r = sum(input["arrival_rates"])
C = input["capacity"]
A = input["requirements_of_routes"]
block_probs = bp_v1(sum_r, C, A)
# make sure that block_prob_i <= block_prob_j => A_i <= A_j
# and block_prob >= block_prob_j => A_i >= A_j
for i, p1 in enumerate(block_probs[:-1]):
for j in range(i + 1, len(block_probs)):
p2 = block_probs[j]
# mystery: changing the inner for loop from below to above resulted in correct test runs
# for j, p2 in enumerate(block_probs[i+1:]):
if p1 < p2:
message = "Route " + str(
i) + " has smaller block prob than route " + str(
j) + " "
message += "(" + str(p1) + " and " + str(
p2) + " respectively),"
message += "but greater resource requirements "
message += "(" + str(A[i]) + " and " + str(
A[j]) + " respectively),"
value = bool(A[i] <= A[j])
self.assertTrue(value, str(message))
elif p1 > p2:
message = "Route " + str(
i) + " has greater block prob than route " + str(
j) + " "
message += "(" + str(p1) + " and " + str(
p2) + " respectively),"
message += "but smaller resource requirements "
message += "(" + str(A[i]) + " and " + str(
A[j]) + " respectively),"
value = bool(A[i] >= A[j])
self.assertTrue(value, str(message))
def test_redundant_lossRate_v2_v3(self):
p1 = ["capacity", "int", [1.0], [10.0], 1]
p2 = ["amount_of_routes", "int", [3.0], [8.0], 1]
p3 = ["arrival_rates", "float", [0.0], [10.0], "amount_of_routes"]
p4 = ["service_rates", "float", [1.0], [1.0], "amount_of_routes"]
p5 = [
"subset_amount_of_routes", "int", [2.0],
["amount_of_routes", "-", 1], 1
]
p6 = [
"requirements_of_routes", "int", [1.0], ["capacity", "/", 2],
"amount_of_routes"
]
input = generate_input([p1, p2, p3, p4, p5, p6])
bp1 = bp_v1(sum(input['arrival_rates']), input['capacity'],
input['requirements_of_routes'])
bp2 = bp_v2(sum(input['arrival_rates']), input['capacity'],
input['requirements_of_routes'])
bp3 = bp_v3(sum(input['arrival_rates']), input['capacity'],
input['requirements_of_routes'])
lossRate1 = rl_v1(bp1, input["arrival_rates"], input["service_rates"])
lossRate2 = rl_v2(bp1, input["arrival_rates"], input["service_rates"])
lossRate3 = rl_v3(bp3, input["arrival_rates"], input["service_rates"])
assert lossRate1 == lossRate3, str(input) + "lossRate1: " + str(
lossRate1) + ", lossRate3: " + str(lossRate3)
assert lossRate1 == lossRate2, str(input) + "lossRate1: " + str(
lossRate1) + ", lossRate2: " + str(lossRate2)
def test_permutations_bp(self):
input = get_general_input()
bp0 = bp_v1(sum(input['arrival_rates']), input['capacity'],
input['requirements_of_routes'])
list_of_tuples = []
for i in range(len(input['requirements_of_routes'])):
list_of_tuples.append((input['requirements_of_routes'][i], bp0[i]))
list_perm = list(itertools.permutations(list_of_tuples))
for i in range(len(list_perm)):
r_r = [elem[0] for elem in list_perm[i]]
bp_temp = bp_v1(sum(input['arrival_rates']), input['capacity'],
r_r)
bp = [elem[1] for elem in list_perm[i]]
self.assertEqual(bp, bp_temp,
"BP: " + str(bp) + "new BP" + str(bp_temp))
def test_variations_lossRate(self):
input = get_general_input()
bp1 = bp_v1(sum(input['arrival_rates']), input['capacity'],
input['requirements_of_routes'])
bp2 = bp_v2(sum(input['arrival_rates']), input['capacity'],
input['requirements_of_routes'])
bp3 = bp_v3(sum(input['arrival_rates']), input['capacity'],
input['requirements_of_routes'])
lossRate1 = rl_v1(bp1, input["arrival_rates"], input["service_rates"])
lossRate2 = rl_v2(bp1, input["arrival_rates"], input["service_rates"])
lossRate3 = rl_v3(bp3, input["arrival_rates"], input["service_rates"])
assert lossRate1 == lossRate3, str(input) + "lossRate1: " + str(
lossRate1) + ", lossRate3: " + str(lossRate3)
assert lossRate1 == lossRate2, str(input) + "lossRate1: " + str(
lossRate1) + ", lossRate2: " + str(lossRate2)
def test_permutations_lossRate(self):
input = get_general_input()
print(str(input))
bp0 = bp_v1(sum(input['arrival_rates']), input['capacity'],
input['requirements_of_routes'])
loss_rate0 = rl_v1(bp0, input['arrival_rates'], input['service_rates'])
list_of_tuples = []
for i in range(len(bp0)):
list_of_tuples.append(
(bp0[i], input['arrival_rates'][i], input['service_rates'][i]))
list_perm = list(itertools.permutations(list_of_tuples))
for i in range(len(list_perm)):
bp = [elem[0] for elem in list_perm[i]]
arr_rates = [elem[1] for elem in list_perm[i]]
ser_rates = [elem[2] for elem in list_perm[i]]
loss_rate = rl_v1(bp, arr_rates, ser_rates)
self.assertEqual(
loss_rate0, loss_rate, "lossrate0: " + str(loss_rate0) +
"\nlossrate: " + str(loss_rate))
def test_redundant_bp_v3(self):
p1 = ["capacity", "int", [1.0], [10.0], 1]
p2 = ["amount_of_routes", "int", [3.0], [8.0], 1]
p3 = ["arrival_rates", "float", [0.0], [10.0], "amount_of_routes"]
p4 = ["service_rates", "float", [1.0], [1.0], "amount_of_routes"]
p5 = [
"subset_amount_of_routes", "int", [2.0],
["amount_of_routes", "-", 1], 1
]
p6 = [
"requirements_of_routes", "int", [1.0], ["capacity", "/", 2],
"amount_of_routes"
]
input = generate_input([p1, p2, p3, p4, p5, p6])
bp1 = bp_v1(sum(input['arrival_rates']), input['capacity'],
input['requirements_of_routes'])
bp3 = bp_v3(sum(input['arrival_rates']), input['capacity'],
input['requirements_of_routes'])
#assert bp1 == bp2, str(input) + "bp1: " + str(bp1) + "bp2: " + str(bp2)
assert bp1 == bp3, str(input) + "bp1: " + str(bp1) + "bp3: " + str(bp3)