def __init__(self) : roadmap = ACC2014Square() super(ACC2014Distr,self).__init__(roadmap, 'length' ) rategraph = nx.DiGraph() rategraph.add_edge( 'N', 'E', rate=1./5 ) rategraph.add_edge( 'W', 'E', rate=1./5 ) rategraph.add_edge( 'W', 'S', rate=3./5 ) self._meancarry = roadcplx.carryMileageRate( roadmap, rategraph ) self._meanfetch = roadcplx.fetchMileageRate( roadmap, rategraph ) # for sampler weights = { (r1,r2) : data.get( 'rate' ) for r1,r2, data in rategraph.edges_iter( data=True ) } #weights[ ('N','E') ] = 1./5 #weights[ ('W','E') ] = 1./5 #weights[ ('W','S') ] = 3./5 self.wset = roadprob.WeightedSet(weights)
probe.horizon = 10000. complexity_computed = [] complexity_estimated = [] def showresults() : plt.scatter( complexity_computed, complexity_estimated ) for t in range(1) : roadnet, rategraph = get_sim_setting( mu=2. ) R = totalrate( rategraph ) n_rategraph = scale_rates( rategraph, 1. / R ) enroute_velocity = carryMileageRate( roadnet, n_rategraph ) balance_velocity = fetchMileageRate( roadnet, n_rategraph ) complexity = enroute_velocity + balance_velocity #complexity = moverscomplexity( roadnet, n_rategraph ) for k in range(1) : numveh = np.random.randint(1,5+1) rate_predicted = convert_complexity_and_servicerate( complexity, numveh ) simrates = scale_rates( n_rategraph, .99 * rate_predicted ) probe.scenario( roadnet, simrates, numveh, 1. ) probe.sim_init() probe.run_sim() rate_observed = probe.rate_observed()