sim.make_simu(storymode=False)
sim.queue = post_process_queue(sim.queue)
G = sim.G
results_queue = extract_aggregate_values_on_queue(sim.queue, paras['par'])
queue_res = [f for f in sim.queue if G.G_nav.has_node(f.source) and G.G_nav.has_node(f.destination) and G.G_nav.node[f.source]['sec']!=paras['STS'] and G.G_nav.node[f.destination]['sec']!=paras['STS']]
results_queue_res = extract_aggregate_values_on_queue(queue_res, paras['par'])
results_G = extract_aggregate_values_on_network(sim.G)
for met in results_G:
results[met]=results_G[met]
for met in results_queue_res:
results_queue[met + '_res'] = results_queue_res[met]
for met in results_queue:
results[met]={tuple(p):[] for p in paras['par']}
for p in paras['par']:
results[met][tuple(p)]=results_queue[met][tuple(p)]
del sim
return results
if __name__=='__main__':
if yes('Ready?'):
results = iter_sim(read_paras_iter(), do=do_shocks, build_pat=build_path)
print 'Done.'
# Used for debugging.
# if 0:
# with open('state.pic','w') as f:
# pickle.dump(getstate(),f)
# else:
# with open('state.pic','r') as f:
# setstate(pickle.load(f))
print header(paras)
if paras['fixnetwork']:
G = paras['G']
else:
G = None
loop({p:paras[p + '_iter'] for p in paras['paras_to_loop']}, paras['paras_to_loop'], \
paras, thing_to_do=average_sim, paras=paras, G=G, do=do, build_pat=build_pat, save=save, rep=rep)
if __name__=='__main__':
"""
Manual entry
"""
paras_file = None if len(sys.argv)==1 else sys.argv[1]
paras = read_paras_iter(paras_file=paras_file)
if yes('Ready?'):
results = iter_sim(paras, save=1)
print 'Done.'
