Commit 1ff66a6a authored by Antoine MARTIN's avatar Antoine MARTIN

modifs fractions intermédiaires 2

parent 1c51da8e
......@@ -38,7 +38,7 @@ for c in [100,0.01]:
print("échantillon "+str(i)+'\n')
fname = "c="+str(c)+"/c="+str(c)+"ech"+str(i)+".hdf5"
fi = h5py.File(fname, "w")
im,centers_diam,f0=image_init(N,D,0.8)
im,centers_diam,f0=image_init(N,D,0.9)
l=len(centers_diam)
fi['centers']=np.zeros((8,l,3))
fi['images']=np.zeros((8,N,N))
......@@ -51,10 +51,7 @@ for c in [100,0.01]:
for j,f in enumerate([0.8,0.6,0.5,0.4,0.3,0.2,0.1,0.05]):
os.system("mkdir c="+str(c)+"/ech="+str(i)+"/frac="+str(f))
print("fraction volumique "+str(f)+'\n')
if f!=0.8:
f1=image_diminue(im,centers_diam,f)
else:
f1=f0
f1=image_diminue(im,centers_diam,f)
# plt.figure()
# plt.imshow(im[:,:],origin='lower')
# plt.show()
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment