paul-corbalan/Projet_Modelisation_mathematique_de_la_foudre
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Switch MAXDELTA and N

Browse Source 
Add dossier variable
This commit is contained in:
Paul-Corbalan 2020-04-03 13:45:43 +02:00
parent 1a69c5deb7
commit 955dd2475b
1 changed files with 18 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

27
Simulation_DBM/Simulation_DBM_2.py
View File

@ -238,18 +238,23 @@ def Print_plt(Mat, N, eta, nbPart, fichier, afficher=False):
plt.title("N="+str(N)+" eta="+str(eta)+" nbpart="+str(nbPart)+" MAXDELTA="+str(MAXDELTA))
# plt.legend(loc='upper left')
plt.axis([0,N,0,N])
plt.savefig("..\\Image\\"+fichier+".png")
plt.savefig(fichier+".png")
if afficher:
plt.show()
plt.close()
# [0.001, 0.0005, 0.0001]
for MAXDELTA in [0.001, 0.0005, 0.0001]:
eta = 6
# [64, 128, 256, 512]
for N in [64, 128, 256, 512]:
# MAXDELTA = MAXDELTA/1000.0
# [64, 128, 256, 512]
for N in [64, 128, 256, 512]:
# [1, 0.1, 0.05, 0.017, 0.01, 0.005, 0.001, 0.0005, 0.0001]
for MAXDELTA in [0.001, 0.0005, 0.0001]:
for i in range(50):
nbPart=1
nbPart = 1
dossier = "..\\Image\\"+"eta="+str(eta)+" N="+str(N)+" MAXDELTA="+str(MAXDELTA)+"\\"
# la grille est de NxN mais en réalite la mesure d'une unité h est 1/N
pixelsBitmap = [[0 for j in range(N)] for i in range(N)]
# que pour le dessin pour éviter la transposition
@ -262,7 +267,11 @@ for MAXDELTA in [0.001, 0.0005, 0.0001]:
nbPart=main(pixelsBitmap, pixelsCroissance, potentielsCroissance, nbPart)
fichier = "foudre-"+str(N)+"-"+str(eta)+"-"+str(MAXDELTA)+"-"+str(nbPart)+"-"+datetime.now().strftime("%Y%m%d%H%M%S")
np.save("..\\Image\\"+fichier,np.array(dessinBitmap))
Save_txt(np.array(dessinBitmap), N, "..\\Image\\"+fichier)
# fichier = "foudre-"+str(N)+"-"+str(eta)+"-"+str(MAXDELTA)+"-"+str(nbPart)+"-"+datetime.now().strftime("%Y%m%d%H%M%S")
# eta=6 N=128 MAXDELTA=0,0001
fichier = dossier+"foudre-"+str(N)+"-"+str(eta)+"-"+str(MAXDELTA)+"-"+str(nbPart)+"-"+datetime.now().strftime("%Y%m%d%H%M%S")
np.save(fichier,np.array(dessinBitmap))
Save_txt(np.array(dessinBitmap), N, fichier)
Print_plt(np.array(dessinBitmap), N, eta, nbPart, fichier)