size = np.zeros(len(bonds))
s = np.zeros(len(bonds))
pvec = np.zeros(len(bonds))
j = 0
for i in range(nn):
for bond in bonds[i * step:(i + 1) * step]:
sites.activate([bond])
P[j] = sites.giantComponent
s[j] = sites.averageSquaredSize
size[j] = sites.sizeOfLargestCluster
pvec[j] = j / len(bonds)
j += 1
print("start: {}, stop: {}, index: {}".format(i * step, (i + 1) * step, i))
image = sites.makeImage()
ax.imshow(image,
aspect="equal",
origin="upper",
vmin=0,
vmax=1,
cmap="Purples")
axes[i].imshow(image,
aspect="equal",
origin="upper",
vmin=0,
vmax=1,
cmap="Purples")
# p = np.sum(np.sum(image))/(L*L)
# print(p)
