Ejemplo n.º 1
0
ncolor = 0
nmarker = 0

Eta = []
Eps = []
Com = []
Efr = []


max_com = 0
max_efr = 0
min_com = 1e30
min_efr = 1e30

for filename in [filename1, filename2]:
    (eta, eps) = classique.lecture(filename, 0, 1)
    (com, efr) = classique.lecture(filename, 2, 4)

    compt = 0
    offset = 6
    for i in range(0, 3):
        Eta.append(eta[compt + 0])
        Eps.append(eps[compt + 0 : compt + offset])
        Com.append(com[compt + 0 : compt + offset])
        Efr.append(efr[compt + 0 : compt + offset])

        max_com = max(max_com, max(com[compt + 0 : compt + offset]))
        max_efr = max(max_efr, max(efr[compt + 0 : compt + offset]))
        min_com = min(min_com, min(com[compt + 0 : compt + offset]))
        min_efr = min(min_efr, min(efr[compt + 0 : compt + offset]))
Ejemplo n.º 2
0
	outputname="/".join(filename.split("/")[0:-1])+"/graphe_tmv_"+(filename.split("/")[-1]).split(".")[0]
else:
	outputname="graphe_tmv_"+(filename.split("/")[-1]).split(".")[0]

print(filename+" will be printed in "+outputname)

########################################################################################
#-------------------------------          Figure           -----------------------------
########################################################################################

colors=["m","b","c","r","g","y","k","firebrick","purple"]
markers=["v","o","*","^","."]
ncolor=0
nmarker=0

(eta,eps)   = classique.lecture(filename,0,1)
(com,tmv)   = classique.lecture(filename,2,5)

Eta  = []
Eps  = []
Com  = []
Tmv  = []

compt=0
max_com = 0
max_tmv = 0
min_com = 1e30
min_tmv = 1e30


offset = 6
########################################################################################
#-------------------------------          Figure           -----------------------------
########################################################################################

colors=["m","b","c","r"]
markers=["^","o",".","v"]
ncolor=0
nmarker=0



nc = 450*3
nr = 450*3
surface = np.zeros((nr,nc))
classique.lecturesurface(filename,surface)



# 				for j in range(0,len(ordonnee)):
# 					surface[int(round(abscisse_1[j]))+(int(1./(2*lc))-1)/2,int(round(abscisse_2[j]))+(int(1./(2*lc))-1)/2]=ordonnee[j]
# 					mask[int(round(abscisse_1[j]))+(int(1./(2*lc))-1)/2,int(round(abscisse_2[j]))+(int(1./(2*lc))-1)/2]=0
					#surface[int(round(abscisse_1[j])),int(round(abscisse_2[j]))]=ordonnee[j]
					#mask[int(round(abscisse_1[j])),int(round(abscisse_2[j]))]=0
				#print(abscisse_1[j],abscisse_2[j],int(round(abscisse_1[j])),int(round(abscisse_2[j])),ordonnee[j])
				#print(mask)
				#print(surface)
				#if maillage==3:
# 				xextent=[lc*min(abscisse_1)-lc*(int(1./(2*lc))-1)/2,lc*max(abscisse_1)+lc*(int(1./(2*lc))-1)/2]
# 				yextent=[lc*min(abscisse_2)-lc*(int(1./(2*lc))-1)/2,lc*max(abscisse_2)+lc*(int(1./(2*lc))-1)/2]
				#else:
filename_refb[1] = "sparseMatrix"
filename_ref[-1] = "_".join(filename_refb)
filename_ref = "/".join(filename_ref)

print(filename + " will be printed in " + outputname)

########################################################################################
# -------------------------------          Figure           -----------------------------
########################################################################################

colors = ["m", "b", "c", "r", "g", "y", "k", "firebrick", "purple"]
markers = ["v", "o", "*", "^", "."]
ncolor = 0
nmarker = 0

(eta, eps) = classique.lecture(filename, 0, 1)
(com, efr) = classique.lecture(filename, 2, 4)

(com_ref, efr_ref) = classique.lecture(filename_ref, 0, 2)

Eta = []
Eps = []
Com = []
Efr = []

compt = 0
max_com = 0
max_efr = 0
min_com = 1e30
min_efr = 1e30
markers=["^","o",".","v"]

nlinestyle=0
ymax_err=0
ymin_err=1e30

Courbes_dist = []
Courbes_rank = []
Courbes_err1 = []
Courbes_err2 = []

for filename in filenames:
	test="_".join((filename.split("/")[-1]).split(".")[0].split("_")[3:])


	(dist,rank)   = classique.lecture(filename,0,1)
	(err1,err2)   = classique.lecture(filename,2,3)


	Dist  = []
	Rank  = []
	Err1  = []
	Err2  = []



	compt=0


	offset = 49
	for i in range(1,10):