def network_simulation():
nw = load_network('data/queueing_params.mat')
target = get_availabilities(nw.station_names)
bal_rates, bal_routing = nw.balance()
nw.combine()
n = Network(nw.size, nw.rates, nw.travel_times, nw.routing, [20]* nw.size)
for i in range(100):
if i % 10 == 0:
print i
n.jump()
T()
def network_simulation():
nw = load_network('data/queueing_params.mat')
target = get_availabilities(nw.station_names)
bal_rates, bal_routing = nw.balance()
nw.combine()
n = Network(nw.size, nw.rates, nw.travel_times, nw.routing, [20] * nw.size)
for i in range(100):
if i % 10 == 0:
print i
n.jump()
T()
def cal_logo_experiment(adj):
nw = load_network('data/queueing_params.mat')
target = get_availabilities(nw.station_names)
bal_rates, bal_routing = nw.balance()
nw.combine()
res = []
for i in adj:
nw.update_adjacency(i)
att_rates, att_routing = nw.min_attack(target, full_adj=False)
T()
res.append(int(np.sum(att_rates)))
print 'Passenger Arrival Rate:', np.sum(nw.rates)
print 'Balance Cost: ', np.sum(bal_rates)
print 'Attack After Balance Cost (adjacency {}): {}'.format(adj, res)
return res
def cal_logo_experiment(adj):
nw = load_network('data/queueing_params.mat')
target = get_availabilities(nw.station_names)
bal_rates, bal_routing = nw.balance()
nw.combine()
res = []
for i in adj:
nw.update_adjacency(i)
att_rates, att_routing = nw.min_attack(target, full_adj=False)
T()
res.append(int(np.sum(att_rates)))
print 'Passenger Arrival Rate:', np.sum(nw.rates)
print 'Balance Cost: ', np.sum(bal_rates)
print 'Attack After Balance Cost (adjacency {}): {}'.format(adj, res)
return res
