color_dict = colors.make_core_cmap(overlaps, cmap='tab20c', seed=-1) frame_list = MOD.loops[this_simname].frame_list[0:2] CHT.plot_2d(htool, core_list=overlaps, frames=frame_list, accumulate=True, label_cores=[0], prefix="S%02d" % nset, color_dict=color_dict) import supersets reload(supersets) if 0: for this_simname in sim_list: new_super = supersets.superset(MOD.loops[this_simname], ht_set[this_simname]) new_super.find() #supersets = st_set[this_simname] #supersets.supersets for nset, this_superset in enumerate(new_super.supersets): print(this_superset) #if 124 not in this_superset: # continue color_dict = colors.make_core_cmap(this_superset, cmap='tab20c', seed=-1) frame_list = MOD.loops[this_simname].frame_list[0:1] CHT.plot_2d(ht_set[this_simname], core_list=this_superset, frames=frame_list, accumulate=True,
if 0: import convex_hull_tools as CHT reload(CHT) if 'ht' not in dir(): ht = {} for this_simname in sim_list: ht[this_simname] = CHT.hull_tool(set_looper) ht[this_simname].make_hulls() ht[this_simname].make_overlaps() import supersets reload(supersets) if 'st' not in dir(): st = {} for this_simname in sim_list: st[this_simname] = supersets.superset(set_looper, ht[this_simname]) st[this_simname].find() do_supersets = True else: do_supersets = False import colors import core_proj reload(core_proj) if do_supersets: stuff = st[this_simname].supersets for nset, superset in enumerate(stuff): if nset != 1: continue core_list = list(superset)
ht[this_simname] = CHT.hull_tool(TL.loops[this_simname]) ht[this_simname].make_hulls() ht[this_simname].make_overlaps() if 'ct' not in dir(): ct = {} for this_simname in sim_list: ct[this_simname] = close_tool.close_tool(TL.loops[this_simname]) ct[this_simname].make_distance() import supersets reload(supersets) if 'st' not in dir(): st = {} for this_simname in sim_list: st[this_simname] = supersets.superset(TL.loops[this_simname], ht[this_simname]) st[this_simname].find() if 'overlap_dict' not in dir(): overlap_dict = {} overlap_numb = {} particles = {} ratio_matrix = {} ratio_matrixb = {} particle_matrix = {} for ns, this_simname in enumerate(sim_list): htool = ht[this_simname] overlap_dict[this_simname] = np.zeros([len(htool.cores_used)] * 2) - 1 overlap_numb[this_simname] = np.zeros([len(htool.cores_used)] * 2) - 1 ratio_matrix[this_simname] = np.zeros([len(htool.cores_used)] * 2)