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()
