77 lines
2.5 KiB
Python
77 lines
2.5 KiB
Python
import matplotlib.pyplot as plt
|
|
import sys
|
|
|
|
vector_width = 8 # 8 doubles per zmm vector
|
|
|
|
# Filter condition of which data to plot
|
|
def plot_filter(atoms_per_unit_cell):
|
|
#return atoms_per_unit_cell < 2048
|
|
return True
|
|
|
|
filename = sys.argv[1]
|
|
plot_output_file = filename.replace(".txt", ".pdf")
|
|
raw_output_file = filename.replace(".txt", ".csv")
|
|
fig = plt.figure()
|
|
ax = plt.axes()
|
|
plot_data = {}
|
|
|
|
with open(filename, 'r') as fp:
|
|
for line in fp.readlines():
|
|
steps, unit_cells, atoms_per_unit_cell, total_atoms, total_vol, atoms_vol, neigh_vol, time, atom_upds_per_sec, cy_per_atom, cy_per_neigh = line.split(',')
|
|
atoms_per_unit_cell = int(atoms_per_unit_cell)
|
|
vol = float(neigh_vol)
|
|
cy_per_atom = float(cy_per_atom)
|
|
|
|
if plot_filter(atoms_per_unit_cell):
|
|
if atoms_per_unit_cell not in plot_data:
|
|
plot_data[atoms_per_unit_cell] = {}
|
|
|
|
cy_per_iter = cy_per_atom * vector_width / atoms_per_unit_cell
|
|
plot_data[atoms_per_unit_cell][vol] = cy_per_iter if vol not in plot_data[atoms_per_unit_cell] \
|
|
else min(cy_per_iter, plot_data[atoms_per_unit_cell][vol])
|
|
|
|
all_volumes = set()
|
|
all_configs = set()
|
|
for atoms_per_unit_cell in plot_data:
|
|
volumes = list(plot_data[atoms_per_unit_cell].keys())
|
|
volumes.sort()
|
|
cycles = [plot_data[atoms_per_unit_cell][vol] for vol in volumes]
|
|
ax.plot(volumes, cycles, marker='.', label=str(atoms_per_unit_cell))
|
|
|
|
for vol in volumes:
|
|
all_volumes.add(vol)
|
|
|
|
all_configs.add(atoms_per_unit_cell)
|
|
|
|
all_volumes = list(all_volumes)
|
|
all_volumes.sort()
|
|
all_configs = list(all_configs)
|
|
all_configs.sort()
|
|
with open(raw_output_file, 'w') as wp:
|
|
wp.write(" volume\config")
|
|
|
|
for conf in all_configs:
|
|
wp.write(",{0:14}".format(conf))
|
|
|
|
wp.write("\n")
|
|
|
|
for vol in all_volumes:
|
|
wp.write("{:14.6f}".format(vol))
|
|
for conf in all_configs:
|
|
try:
|
|
cycles = plot_data[conf][vol]
|
|
wp.write(",{:14.6f}".format(cycles))
|
|
except:
|
|
wp.write(',' + ' ' * 14)
|
|
|
|
|
|
wp.write("\n")
|
|
|
|
ax.vlines([32, 1000, 28000], 0, 1, transform=ax.get_xaxis_transform(), linestyles='dashed', color=['#444444', '#777777', '#aaaaaa'])
|
|
ax.set(xlabel='Neighbor data volume (kB)', ylabel='Cycles per iteration')
|
|
ax.set_xscale('log')
|
|
#ax.set_xticks([32, 1000, 28000])
|
|
#ax.set_xlim(0, 200000)
|
|
plt.legend(title="atoms/uc")
|
|
fig.savefig(plot_output_file, bbox_inches = 'tight', pad_inches = 0)
|