def test_add_mode(self):
"""
to test add_mode method
:return:
"""
m = TransportModeManager()
name = "train"
bya = 1
co = 1
c1 = 1
c2 = 1
v = 20
t = 0.2
fmax = 100
kmax = 50
theta = 0.5
tat = 1
d = 1
fini = 28
mode_obj = TransportMode(name, bya, co, c1, c2, v, t, fmax, kmax, theta, tat, d, fini)
m.add_mode(mode_obj)
self.assertEqual(len(m.get_modes()), 3)
with self.assertRaises(exceptions.AddModeExceptions):
mode_obj = TransportMode("bus", bya, co, c1, c2, v, t, fmax, kmax, theta, tat, d, fini)
m.add_mode(mode_obj)
with self.assertRaises(exceptions.ModeIsNotValidException):
m.add_mode("mode_obj")
def test_fmax_constrains(self):
graph_obj = Graph.build_from_parameters(2, 10, 1, 2)
network_obj = TransportNetwork(graph_obj)
[bus_obj, metro_obj] = TransportMode.get_default_modes()
radial_bus = network_obj.get_radial_routes(bus_obj)
feeder_bus = network_obj.get_feeder_routes(bus_obj)
radial_metro = network_obj.get_radial_routes(metro_obj, short=True)
for route in radial_bus:
network_obj.add_route(route)
for route in feeder_bus:
network_obj.add_route(route)
for route in radial_metro:
network_obj.add_route(route)
f = defaultdict(float)
for route in network_obj.get_routes():
f[route.id] = 28
ineq_f = Constrains.fmax_constrains(graph_obj,
network_obj.get_routes(),
network_obj.get_modes(), f)
self.assertEqual(ineq_f, [
-136.0, -40.0, -136.0, -40.0, -143.0, -12.0, -143.0, -12.0, -150.0,
-40.0, -150.0, -40.0, -136.0, -40.0, -136.0, -40.0, -143.0, -12.0,
-143.0, -12.0
])
def setUp(self) -> None:
graph_obj = Graph.build_from_parameters(n=2, l=10, g=0.5, p=2)
demand_obj = Demand.build_from_parameters(graph_obj=graph_obj, y=100, a=0.5, alpha=1 / 3, beta=1 / 3)
passenger_obj = Passenger.get_default_passenger()
[bus_obj, metro_obj] = TransportMode.get_default_modes()
self.network_obj = TransportNetwork(graph_obj=graph_obj)
feeder_routes_metro = self.network_obj.get_feeder_routes(mode_obj=metro_obj)
radial_routes_bus = self.network_obj.get_radial_routes(mode_obj=bus_obj)
for route in feeder_routes_metro:
self.network_obj.add_route(route_obj=route)
for route in radial_routes_bus:
self.network_obj.add_route(route_obj=route)
self.extended_graph_obj = ExtendedGraph(graph_obj=graph_obj, routes=self.network_obj.get_routes(),
TP=passenger_obj.pt, frequency_routes=None)
hyperpath_obj = Hyperpath(extended_graph_obj=self.extended_graph_obj, passenger_obj=passenger_obj)
self.hyperpaths, self.labels, self.successors, self.frequency, self.Vij = hyperpath_obj.get_all_hyperpaths(
OD_matrix=demand_obj.get_matrix())
self.OD_assignment = Assignment.get_assignment(hyperpaths=self.hyperpaths, labels=self.labels, p=2,
vp=passenger_obj.va, spa=passenger_obj.spa,
spv=passenger_obj.spv)
def test_most_loaded_section_constrains(self):
graph_obj = Graph.build_from_parameters(2, 10, 1, 2)
network_obj = TransportNetwork(graph_obj)
[bus_obj, metro_obj] = TransportMode.get_default_modes()
radial_bus = network_obj.get_radial_routes(bus_obj)
feeder_bus = network_obj.get_feeder_routes(bus_obj)
radial_metro = network_obj.get_radial_routes(metro_obj, short=True)
for route in radial_bus:
network_obj.add_route(route)
for route in feeder_bus:
network_obj.add_route(route)
for route in radial_metro:
network_obj.add_route(route)
most_loaded_section = defaultdict(float)
for route in network_obj.get_routes():
most_loaded_section[route.id] = 10
ineq_k = Constrains.most_loaded_section_constrains(
network_obj.get_routes(), most_loaded_section)
self.assertEqual(len(ineq_k), 6)
self.assertEqual(ineq_k, [-150, -150, -150, -150, -1430, -1430])
def test_get_mode_network_distance(self):
graph_obj = Graph.build_from_parameters(2, 10, 1, 2)
network_obj = TransportNetwork(graph_obj)
[bus_obj, metro_obj] = TransportMode.get_default_modes()
CI_obj = InfrastructureCost()
circular = network_obj.get_circular_routes(mode_obj=metro_obj)
radial = network_obj.get_radial_routes(mode_obj=bus_obj)
for route in circular:
network_obj.add_route(route)
for route in radial:
network_obj.add_route(route)
f = defaultdict(float)
for route in network_obj.get_routes():
f[route.id] = 28
mode_distance = CI_obj.get_mode_network_distance(graph_obj,
network_obj,
f=f)
self.assertEqual(mode_distance[bus_obj], 160)
self.assertEqual(mode_distance[metro_obj], 20)
def test_get_operators_cost(self):
"""
to test operators cost method
:return:
"""
graph_obj = Graph.build_from_parameters(n=2, l=10, g=0.5, p=2)
[bus_obj, metro_obj] = TransportMode.get_default_modes()
network_obj = TransportNetwork(graph_obj=graph_obj)
feeder_routes_metro = network_obj.get_feeder_routes(mode_obj=metro_obj)
radial_routes_bus = network_obj.get_radial_routes(mode_obj=bus_obj)
for route in feeder_routes_metro:
network_obj.add_route(route)
for route in radial_routes_bus:
network_obj.add_route(route)
f = defaultdict(float)
k = defaultdict(float)
line_cycle_time = defaultdict(float)
for route in network_obj.get_routes():
f[route.id] = 28
k[route.id] = 100
line_cycle_time[route.id] = 10
cost = OperatorsCost.get_operators_cost(network_obj.get_routes(),
line_cycle_time, f, k)
self.assertEqual(cost, 75331.2)
def test_lines_travel_time(self):
"""
test lines_travel_time method of class operators_cost
:return:
"""
graph_obj = Graph.build_from_parameters(n=2, l=10, g=0.5, p=2)
[bus_obj, metro_obj] = TransportMode.get_default_modes()
network_obj = TransportNetwork(graph_obj=graph_obj)
feeder_routes_metro = network_obj.get_feeder_routes(mode_obj=metro_obj)
radial_routes_bus = network_obj.get_radial_routes(mode_obj=bus_obj)
for route in feeder_routes_metro:
network_obj.add_route(route_obj=route)
for route in radial_routes_bus:
network_obj.add_route(route_obj=route)
lines_travel_time = OperatorsCost.lines_travel_time(
routes=network_obj.get_routes(),
edge_distance=graph_obj.get_edges_distance())
self.assertEqual(lines_travel_time["R_bus_1"], 1.5)
self.assertEqual(lines_travel_time["F_metro_1"], 0.25)
def test_resources_consumer(self):
graph_obj = Graph.build_from_parameters(n=2, l=10, g=0.5, p=2)
demand_obj = Demand.build_from_parameters(graph_obj=graph_obj,
y=100,
a=0.5,
alpha=1 / 3,
beta=1 / 3)
passenger_obj = Passenger.get_default_passenger()
[bus_obj, metro_obj] = TransportMode.get_default_modes()
network_obj = TransportNetwork(graph_obj=graph_obj)
feeder_routes_metro = network_obj.get_feeder_routes(mode_obj=metro_obj)
radial_routes_bus = network_obj.get_radial_routes(mode_obj=bus_obj)
for route in feeder_routes_metro:
network_obj.add_route(route_obj=route)
for route in radial_routes_bus:
network_obj.add_route(route_obj=route)
extended_graph_obj = ExtendedGraph(graph_obj=graph_obj,
routes=network_obj.get_routes(),
TP=passenger_obj.pt,
frequency_routes=None)
hyperpath_obj = Hyperpath(extended_graph_obj=extended_graph_obj,
passenger_obj=passenger_obj)
hyperpaths, labels, successors, frequency, Vij = hyperpath_obj.get_all_hyperpaths(
OD_matrix=demand_obj.get_matrix())
OD_assignment = Assignment.get_assignment(hyperpaths=hyperpaths,
labels=labels,
p=2,
vp=passenger_obj.va,
spa=passenger_obj.spa,
spv=passenger_obj.spv)
f = defaultdict(float)
z = defaultdict(lambda: defaultdict(lambda: defaultdict(float)))
v = defaultdict(lambda: defaultdict(lambda: defaultdict(float)))
for route in network_obj.get_routes():
f[route.id] = 28
CU_obj = UsersCost()
ta, te, tv, t = CU_obj.resources_consumer(hyperpaths, Vij,
OD_assignment, successors,
extended_graph_obj,
passenger_obj.va, f, z, v)
self.assertEqual(round(ta, 4), 3.069)
self.assertEqual(round(te, 4), 13.9597)
self.assertEqual(round(tv, 4), 74.685)
self.assertEqual(round(t, 4), 41.6667)
def test_raises_mode_exceptions(self):
"""
to test mode class exceptions
:return:
"""
name = "bus"
bya = 1
co = 1
c1 = 1
c2 = 1
v = 20
t = 0.2
fmax = 100
kmax = 50
theta = 0.5
tat = 1
d = 1
fini = 28
with self.assertRaises(exceptions.NameIsNotValidExceptions):
TransportMode(None, bya, co, c1, c2, v, t, fmax, kmax, theta, tat,
d, fini)
with self.assertRaises(exceptions.ByaIsNotValidExceptions):
TransportMode(name, -1, co, c1, c2, v, t, fmax, kmax, theta, tat,
d, fini)
with self.assertRaises(exceptions.CoIsNotValidExceptions):
TransportMode(name, bya, -2, c1, c2, v, t, fmax, kmax, theta, tat,
d, fini)
with self.assertRaises(exceptions.C1IsNotValidExceptions):
TransportMode(name, bya, co, None, c2, v, t, fmax, kmax, theta,
tat, d, fini)
with self.assertRaises(exceptions.C2IsNotValidExceptions):
TransportMode(name, bya, co, c1, -3, v, t, fmax, kmax, theta, tat,
d, fini)
with self.assertRaises(exceptions.VIsNotValidExceptions):
TransportMode(name, bya, co, c1, c2, None, t, fmax, kmax, theta,
tat, d, fini)
with self.assertRaises(exceptions.TIsNotValidExceptions):
TransportMode(name, bya, co, c1, c2, v, -3, fmax, kmax, theta, tat,
d, fini)
with self.assertRaises(exceptions.FmaxIsNotValidExceptions):
TransportMode(name, bya, co, c1, c2, v, t, -200, kmax, theta, tat,
d, fini)
with self.assertRaises(exceptions.KmaxIsNotValidExceptions):
TransportMode(name, bya, co, c1, c2, v, t, fmax, -100, theta, tat,
d, fini)
with self.assertRaises(exceptions.ThetaIsNotValidExceptions):
TransportMode(name, bya, co, c1, c2, v, t, fmax, kmax, 2, tat, d, fini)
with self.assertRaises(exceptions.TatIsNotValidExceptions):
TransportMode(name, bya, co, c1, c2, v, t, fmax, kmax, theta, None,
d, fini)
with self.assertRaises(exceptions.FiniIsNotValidException):
TransportMode(name, bya, co, c1, c2, v, t, fmax,
kmax, theta, tat, d, -1)
spv = 2.74
spw = 2 * 2.74
spa = 3 * 2.74
spt = 16
passenger_obj = Passenger(va=va,
pv=pv,
pw=pw,
pa=pa,
pt=pt,
spv=spv,
spw=spw,
spa=spa,
spt=spt)
[bus, metro] = TransportMode.get_default_modes()
bus.c2 = 92.39
# dc4_metro = network_obj.get_diametral_routes(metro, 4, True)
# tc1_bus = network_obj.get_tangencial_routes(bus, 1, True)
# tc2_bus = network_obj.get_tangencial_routes(bus, 2, True)
# dc3_bus = network_obj.get_diametral_routes(bus, 3, True)
# dc4_bus = network_obj.get_diametral_routes(bus, 4, True)
# rc_bus = network_obj.get_radial_routes(bus, True)
# r_bus = network_obj.get_radial_routes(bus)
# f_bus = network_obj.get_feeder_routes(bus)
# t1_bus = network_obj.get_tangencial_routes(bus, 1)
# t2_bus = network_obj.get_tangencial_routes(bus, 2)
# d3_bus = network_obj.get_diametral_routes(bus, 3)
# d4_bus = network_obj.get_diametral_routes(bus, 4)
# c_bus = network_obj.get_circular_routes(bus)
def test_get_cycle_time(self):
"""
test get_cycle_time method of class operators_cost
:return:
"""
graph_obj = Graph.build_from_parameters(n=2, l=10, g=0.5, p=2)
demand_obj = Demand.build_from_parameters(graph_obj=graph_obj,
y=100,
a=0.5,
alpha=1 / 3,
beta=1 / 3)
passenger_obj = Passenger.get_default_passenger()
[bus_obj, metro_obj] = TransportMode.get_default_modes()
network_obj = TransportNetwork(graph_obj=graph_obj)
feeder_routes_metro = network_obj.get_feeder_routes(mode_obj=metro_obj)
radial_routes_bus = network_obj.get_radial_routes(mode_obj=bus_obj)
for route in feeder_routes_metro:
network_obj.add_route(route_obj=route)
for route in radial_routes_bus:
network_obj.add_route(route_obj=route)
extended_graph_obj = ExtendedGraph(graph_obj=graph_obj,
routes=network_obj.get_routes(),
TP=passenger_obj.pt,
frequency_routes=None)
hyperpath_obj = Hyperpath(extended_graph_obj=extended_graph_obj,
passenger_obj=passenger_obj)
hyperpaths, labels, successors, frequency, Vij = hyperpath_obj.get_all_hyperpaths(
OD_matrix=demand_obj.get_matrix())
OD_assignment = Assignment.get_assignment(hyperpaths=hyperpaths,
labels=labels,
p=2,
vp=passenger_obj.va,
spa=passenger_obj.spa,
spv=passenger_obj.spv)
f = defaultdict(float)
for route in network_obj.get_routes():
f[route.id] = 28
z, v, loaded_section_route = Assignment.get_alighting_and_boarding(
Vij=Vij,
hyperpaths=hyperpaths,
successors=successors,
assignment=OD_assignment,
f=f)
lines_travel_time = OperatorsCost.lines_travel_time(
routes=network_obj.get_routes(),
edge_distance=graph_obj.get_edges_distance())
line_cycle_time = OperatorsCost.get_cycle_time(
z, v, network_obj.get_routes(), lines_travel_time)
self.assertEqual(round(line_cycle_time["F_metro_1"], 7), 0.2500082)
self.assertEqual(round(line_cycle_time["R_bus_1"], 7), 1.5032756)
self.assertEqual(line_cycle_time["l1"], 0)