def exp_4():
net = StateConstrainedNetwork.load('../networks/fpnet_with_demands.json')
net.objective = lambda: 0
l_9 = net.link_by_name('9')
l_3 = net.link_by_name('3')
l_8 = net.link_by_name('8')
for link in net.get_links():
if link in [l_9, l_3, l_8]:
link.set_cong_state(link.CongState.CONG)
else:
link.set_cong_state(link.CongState.FF)
net.get_program().cr_solve()
net.realize()
print net.total_travel_time()
net.dump('../networks/exps/exp4/net_cong.json')
net = StateConstrainedNetwork.load('../networks/exps/exp4/net_cong.json')
net.get_program().cr_solve()
net.realize()
print net.total_travel_time()
net = StateConstrainedComplacentNetwork.load('../networks/exps/exp4/net_cong.json')
net.scale = 1.1
net.get_program().cr_solve()
net.realize()
print net.total_travel_time()
net = StateConstrainedComplacentNetwork.load('../networks/exps/exp4/net_cong.json')
net.scale = 1000.0
net.get_program().cr_solve()
net.realize()
print net.total_travel_time()
net = StateConstrainedComplacentNetwork.load('../networks/exps/exp4/net_cong.json')
net.scale = 10.0
net.get_program().cr_solve()
net.realize()
print net.total_travel_time()
def exp_3_info():
net = StateConstrainedComplacentNetwork.load('../networks/exps/exp3/net_demand.json')
print 'ttt previous', net.total_travel_time()
net.get_program().cr_solve()
net.realize()
print 'ttt after complacence optimize', net.total_travel_time()
for route in net.all_routes():
print 'route', route
print 'heur tt', net.route_tt_heuristic(route).value
print 'actual tt', route.travel_time()
net = StateConstrainedNetwork.load('../networks/exps/exp3/net_demand.json')
net.get_program().cr_solve()
net.realize()
print 'ttt after non-comp optimize', net.total_travel_time()
def exp_3_info():
net = StateConstrainedComplacentNetwork.load('../networks/exps/exp3/net_demand.json')
for cls in net.__class__.__mro__:
print cls
print 'ttt previous', net.total_travel_time()
net.get_program().cr_print()
net.get_program().cr_solve(quiet=False)
for link in net.get_links():
print link
print link.v_dens.value
print link.l
net.realize()
print 'ttt after complacence optimize', net.total_travel_time()
for route in net.all_routes():
print 'route', route
print 'heur tt', net.route_tt_heuristic(route).value
print 'actual tt', net.route_travel_time(route)
net = StateConstrainedNetwork.load('../networks/exps/exp3/net_demand.json')
net.get_program().cr_solve()
net.realize()
print 'ttt after non-comp optimize', net.total_travel_time()
def exp_3_setup():
net = StateConstrainedComplacentNetwork()
l = 1; v = 1; w = 1; q_max = 1; rho_max = 2; r = 1;
fd = FundamentalDiagram.triangular(
v=v,q_max=q_max,rho_max=rho_max
)
l_rho = 1.372281323
r_rho = 3 - l_rho
source = StateConstrainedLink(
fd=fd,
name='source',
cong_state=StateConstrainedLink.CongState.ANY,
l=l,
flow=r,
density=fd.rho_ff(r)
)
left = StateConstrainedLink(
fd=fd,
name='left',
cong_state=StateConstrainedLink.CongState.CONG,
l=2*l,
flow=fd.flow_cong(l_rho),
density=l_rho
)
right = StateConstrainedLink(
fd=fd,
name='right',
cong_state=StateConstrainedLink.CongState.CONG,
l=l,
flow=fd.flow_cong(r_rho),
density=r_rho
)
sink = StateConstrainedLink(
fd=fd,
name='sink',
cong_state=StateConstrainedLink.CongState.ANY,
l=l,
flow=r,
density=fd.rho_ff(r)
)
js = [
Junction([source],[left,right]),
Junction([left,right],[sink]),
]
[net.add_junction(junction) for junction in js]
net.dump('../networks/exps/exp3/net.json')
left_demand = RouteDemand(
route=net.route_by_names(['source','left','sink']),
flow=.2
)
od_demand = ODDemand(
source=source,
sink=sink,
flow=.8
)
net.demands.append(left_demand)
net.demands.append(od_demand)
net.dump('../networks/exps/exp3/net_demand.json')
def exp_2():
"""
starting to figure out stateconstrained things
"""
net = StateConstrainedComplacentNetwork.load('../../networks/exps/exp2/lf_rc.json')
net.get_program().cr_print()
