def main(): # initialize the timer print_current_time() # the data set name case_name = '6node' # load activity data load_activity(case_name) # load network data load_network(case_name) # generate path sets gen_path_set() find_shortest_path() # generate utils gen_solo_activity_util() print '\n DATA LOADED' print_current_time() # run multiple scenarios corr_list = [corr / 10.0 for corr in range(10, -1, -5)] run_multi_scenarios(case_name, corr_list) # export multi-run data export_multi_run_data(case_name)
def find_fixed_point(N, case_name, corr): " Find the equilibrium flows using method of successive average (MSA). " # iterate demand and supply sides init_flow_variables() init_socio_activity_util() for iter_num in xrange(N): print "\n ### interation %d ###" % iter_num print '\n [initialization]' init_step_variables() init_state_util() init_transition_choice_prob() print_current_time() print '\n [dynamic choice]' calc_state_util() build_choice_model() calc_inclusive_values() calc_commodity_steps() print_current_time() print '\n [traffic loading]' init_state_flows(0.0) init_transition_flows(0.0) calc_state_flows() print_current_time() print '\n [intra-household interaction]' init_activity_choice_prob() calc_activity_choice_prob() calc_socio_activity_util() print_current_time() print '\n [update flows]' update_movement_flows(iter_num) update_commodity_flows(iter_num) print_current_time() print "\n [export data]" export_data(case_name + '_r' + str(corr) + '_n' + str(iter_num)) print_current_time()
def find_fixed_point(N, case_name, corr): " Find the equilibrium flows using method of successive average (MSA). " # iterate demand and supply sides init_flow_variables() init_socio_activity_util() for iter_num in xrange(N): print "\n ### interation %d ###" % iter_num print '\n [initialization]' init_step_variables() init_state_util() init_transition_choice_prob() print_current_time() print '\n [dynamic choice]' calc_state_util() build_choice_model() calc_inclusive_values() calc_commodity_steps() print_current_time() print '\n [traffic loading]' init_state_flows(0.0) init_transition_flows(0.0) calc_state_flows() print_current_time() print '\n [intra-household interaction]' init_activity_choice_prob() calc_activity_choice_prob() calc_socio_activity_util() print_current_time() print '\n [update flows]' update_movement_flows(iter_num) update_commodity_flows(iter_num) print_current_time() print "\n [export data]" export_data(case_name+'_r'+str(corr)+'_n'+str(iter_num)) print_current_time()