Skip to content

Commit 0c3713b

Browse files
FarnazHFanwangM
authored andcommitted
Update script to make plot folders for each selection method
1 parent 3cd0009 commit 0c3713b

File tree

1 file changed

+9
-3
lines changed

1 file changed

+9
-3
lines changed

notebooks/generate_data_pes_4well.py

Lines changed: 9 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -185,7 +185,7 @@ def select_subsets(points, percentages, methods, seed=42):
185185
indices = Medoid().select(points, size=size)
186186
elif method.startswith("GridPartition"):
187187
grid_method = method.split("-")[1]
188-
indices = GridPartition(numb_bins_axis=5, grid_method=grid_method).select(
188+
indices = GridPartition(nbins_axis=5, bin_method=grid_method).select(
189189
points, size=size
190190
)
191191
else:
@@ -291,7 +291,7 @@ def plot_data_2d(points, values, highlight=None, title="", fname=None):
291291
"GridPartition-equifrequent_dependent",
292292
"GridPartition-equifrequent_independent",
293293
]
294-
percentages = [2**i for i in range(6)]
294+
percentages = [2**i for i in range(7)]
295295

296296
# Make an instance of potential for visualization
297297
# -----------------------------------------------
@@ -323,12 +323,18 @@ def plot_data_2d(points, values, highlight=None, title="", fname=None):
323323
# indices are returned as a nested dictionary {method: {percentage: indices}}
324324
selected_indices = select_subsets(points_iter, percentages, selectors, seed=42)
325325
database_indices[fname_iter] = selected_indices
326+
print("DONE SELECTING")
326327
# Plot Subsets
327328
# ------------
328329
for method, values in selected_indices.items():
330+
print("PLOTTING.....")
331+
os.makedirs(f"{folder}/plots/{method}", exist_ok=True)
329332
for p, indices in values.items():
330333
p = int(p)
331-
fname_plot = f"{folder}/plots/subset_{method}_p{p:02d}_{fname_iter}.png"
334+
fname_plot = (
335+
f"{folder}/plots/{method}/subset_{method}_p{p:02d}_{fname_iter}.png"
336+
)
337+
print(f"PLOT {fname_plot}")
332338
plot_data_2d(
333339
points_plot,
334340
values_plot,

0 commit comments

Comments
 (0)