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.'