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)
