def hmm_graph(): global hmm_graph_ if hmm_graph_ is None: if not os.path.exists(hmm_graph_fname): if graph_type == 'simple': hmm_graph_ = model.createHMMGraphFromNetwork(net, mode_counts=mode_counts()) else: # Complex model not implemented assert False else: tic("Reading completed hmm graph from %s"%hmm_graph_fname) hmm_graph_ = pickle.load(open(hmm_graph_fname,'r')) return hmm_graph_
def hmm_graph(): global hmm_graph_ if hmm_graph_ is None: if not os.path.exists(hmm_graph_fname): if graph_type == 'simple': tic("creating empty hmm graph", experiment_name) hmm_graph_ = model.createHMMGraphFromNetwork(net, mode_counts=mode_counts()) tic("done creating empty hmm graph", experiment_name) tic("saving hmm graph pickle", experiment_name) pickle.dump(hmm_graph(),open(hmm_graph_fname,'w')) tic("done saving hmm graph pickle", experiment_name) else: # Complex model not implemented assert False else: tic("Reading completed hmm graph from %s"%hmm_graph_fname) hmm_graph_ = pickle.load(open(hmm_graph_fname,'r')) tic("done reading completed hmm graph from %s"%hmm_graph_fname) return hmm_graph_
traj_conv_param = experiment_design['trajectory_conversion']['params'] traj_conv = createTrajectoryConversion(graph_type=graph_type, process=experiment_design['trajectory_conversion']['process'], params=traj_conv_param, network=net, max_nb_mixture=traj_conv_param['max_n_modes']) traj_conv_one_mode = createTrajectoryConversion(graph_type=graph_type, process='mixture_auto', params=traj_conv_param, network=net, max_nb_mixture=1) # mode_counts = dict([(link_id,1) for link_id in net.keys()]) if graph_type == 'simple': hmm_graph = model.createHMMGraphFromNetwork(net, mode_counts=traj_conv.modeCounts()) hmm_graph_one_mode = model.createHMMGraphFromNetwork(net, mode_counts=traj_conv_one_mode.modeCounts()) # hmm_graph = model.createHMMGraphFromNetwork(net, mode_counts=mode_counts) else: # Complex model not implemented assert False tt_graph = createTravelTimeGraph(hmm_graph, radius=2e-4) tt_graph.checkInvariants() tt_graph_one_mode = createTravelTimeGraph(hmm_graph_one_mode, radius=2e-4) tt_graph_one_mode.checkInvariants() gmrf = emptyValues(tt_graph) gmrf_one_mode = emptyValues(tt_graph)