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)
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)