Commit 34ee12d2 authored by Antoine MARTIN's avatar Antoine MARTIN

ajout de fractions volumiques et de cartes erreur de déformations

parent 6d0e87bf
......@@ -27,6 +27,26 @@ def reference(tabC,res1,c,ec,fr):
k,ll,mm,nn=tabC.shape
mesh = UnitSquareMesh(res1, res1, "left/right")
def pixel_0(x,on_boundary):
return x[0]<=0.5 and x[1]<=0.5
pixels=np.array([[pixel_0 for i in range(res1)] for j in range(res1)])
for i in range(res1):
for j in range(res1):
def pixel(x,on_boundary):
return x[0]<=(j+1)/res1 and x[0]>=j/res1 and x[1]<=(i+1)/res1 and x[1]>=i/res1
pixels[i,j]=pixel()
domains = MeshFunction("size_t", mesh,2)
domains.set_all(0)
for i in range(res1):
for j in range(res1):
AutoSubDomain(pixels[i,j]).mark(domains, i*res1+j)
dx = Measure("dx", domain=mesh, subdomain_data=domains)
class LocalC(UserExpression):
def __init__(self, Cij, mesh):
self.Cij = Cij
......@@ -123,13 +143,16 @@ def reference(tabC,res1,c,ec,fr):
a += dot(lamb_,dv)*dx + dot(dlamb,v_)*dx
Chom = np.zeros((3, 3))
tA=np.array([[np.array([[1., 0., 0.],
[0., 1., 0.],
[0., 0., 1.]]) for i in range(res1)] for j in range(res1)])
for (j, case) in enumerate(["Exx", "Eyy", "Exy"]):
ffile = XDMFFile('c='+str(c)+'/ech='+str(ec)+'/frac='+str(fr)+'/Strain_EF_{}.xdmf'.format(case))
ee = np.zeros((3,))
ee[j] = 1
Eps.assign(Constant(ee))
solve(a == L, w, solver_parameters={"linear_solver": "mumps"})#, solver_parameters={"linear_solver": "bicgstab", "preconditioner":"ilu"})
solve(a == L, w)#, solver_parameters={"linear_solver": "mumps"})#, solver_parameters={"linear_solver": "bicgstab", "preconditioner":"ilu"})
(v, lamb) = split(w)
vf = Function(V, name='Strain-'+case)
vf.assign(local_project(eps(v), V))
......@@ -154,15 +177,20 @@ def reference(tabC,res1,c,ec,fr):
# #plt.show()
# plt.close()
for ii in range(res1):
for jj in range(res1):
for kk in range(3):
tA[res1-ii,jj][kk,j]=assemble(eps(v)[kk]*dx(ii*res1+jj))+Eps[kk]
Sigma = np.zeros((3,))
for kk in range(3):
Sigma[kk] = assemble(sigma(v, Eps)[kk]*dx)
Chom[j, :] = Sigma
ffile.close()
return Chom
return Chom,tA
Chom = reference_problem(Cvar, mesh,k)
Chom,tA = reference_problem(Cvar, mesh,k)
CVoigt = np.array([[0., 0., 0.],
[0., 0., 0.],
[0., 0., 0.]])
......@@ -173,7 +201,7 @@ def reference(tabC,res1,c,ec,fr):
CVoigt=CVoigt/k/k
print(Chom)
return Chom,CVoigt
return Chom,CVoigt,tA
......
......@@ -3,6 +3,7 @@ import h5py
import matplotlib.pyplot as plt
from matplotlib import cm
from matplotlib.ticker import LinearLocator, FormatStrFormatter
import os
......@@ -65,7 +66,7 @@ def verif_im(c,i):
def verif_eps(c,i,j,niter):
fi = h5py.File("c="+str(c)+"/c="+str(c)+"ech"+str(i)+".hdf5","r")
EpsHH=fi['EpsHH'][j,niter-1]
frac=[0.05,0.1,0.2,0.3,0.4,0.5,0.6,0.8]
frac=[0.05,0.1,0.2,0.3,0.4,0.5,0.6,0.67,0.8,0.9]
f=frac[j]
plt.figure()
pp=plt.imshow(EpsHH[:,:,0,0],origin='lower')
......@@ -86,7 +87,7 @@ def verif_eps(c,i,j,niter):
#plt.show()
plt.close()
plt.figure()
pp=plt.imshow(EpsHH[:,:,1,1],origin='lower')
pp=plt.imshow(EpsHH[:,:,2,2],origin='lower')
plt.colorbar(pp)
plt.savefig('c='+str(c)+'/ech='+str(i)+'/frac='+str(f)+'/Aloc HH 33 iter='+str(niter)+'.pdf')
#plt.show()
......@@ -97,45 +98,69 @@ def verif_eps(c,i,j,niter):
def courbes(c,i):
os.system('mkdir courbes')
fi = h5py.File("c="+str(c)+"/c="+str(c)+"ech"+str(i)+".hdf5","r")
fractions=fi['fractions']
fr=fi['fractions']
ChomEF=fi['ChomEF']
ChomHH=fi['ChomHH']
C_11=np.zeros((8,))
C_11_EF=np.zeros((8,))
C_22=np.zeros((8,))
C_22_EF=np.zeros((8,))
C_12=np.zeros((8,))
C_12_EF=np.zeros((8,))
C_33=np.zeros((8,))
C_33_EF=np.zeros((8,))
C_11_0=np.zeros((8,))
C_22_0=np.zeros((8,))
C_12_0=np.zeros((8,))
C_33_0=np.zeros((8,))
for j,f in enumerate([0.05,0.1,0.2,0.3,0.4,0.5,0.6,0.8]):
C_11[j]=ChomHH[j,-1,1,0,0]/1.11111*np.sqrt(c)/10
C_11_0[j]=ChomHH[j,-1,0,0,0]/1.11111*np.sqrt(c)/10
C_11_EF[j]=ChomEF[j,0,0]/1.11111*np.sqrt(c)/10
C_22[j]=ChomHH[j,-1,1,1,1]/1.11111*np.sqrt(c)/10
C_22_0[j]=ChomHH[j,-1,0,1,1]/1.11111*np.sqrt(c)/10
C_22_EF[j]=ChomEF[j,1,1]/1.11111*np.sqrt(c)/10
C_12[j]=ChomHH[j,-1,1,0,1]/0.277777*np.sqrt(c)/10
C_12_0[j]=ChomHH[j,-1,0,0,1]/0.277777*np.sqrt(c)/10
C_12_EF[j]=ChomEF[j,0,1]/0.277777*np.sqrt(c)/10
C_33[j]=ChomHH[j,-1,1,2,2]/2/0.416666*np.sqrt(c)/10
C_33_0[j]=ChomHH[j,-1,0,2,2]/2/0.416666*np.sqrt(c)/10
C_33_EF[j]=ChomEF[j,2,2]/2/0.416666*np.sqrt(c)/10
C_11=np.ones((12,))
C_11_EF=np.ones((12,))
C_22=np.ones((12,))
C_22_EF=np.ones((12,))
C_12=np.ones((12,))
C_12_EF=np.ones((12,))
C_33=np.ones((12,))
C_33_EF=np.ones((12,))
C_11_0=np.ones((12,))
C_22_0=np.ones((12,))
C_12_0=np.ones((12,))
C_33_0=np.ones((12,))
fractions=np.ones((12,))
for j,f in enumerate([0.,0.05,0.1,0.2,0.3,0.4,0.5,0.6,0.67,0.8,0.9,1.]):
if c==0.01:
jj=11-j
else:
jj=j
if f==0. or f==1.:
fractions[jj]=f
if jj==11:
C_11[jj]=0.01
C_11_0[jj]=0.01
C_11_EF[jj]=0.01
C_22[jj]=0.01
C_22_0[jj]=0.01
C_22_EF[jj]=0.01
C_12[jj]=0.01
C_12_0[jj]=0.01
C_12_EF[jj]=0.01
C_33[jj]=0.01
C_33_0[jj]=0.01
C_33_EF[jj]=0.01
if f!=0. and f!=1.:
fractions[jj]=fr[j-1]
C_11[jj]=ChomHH[j-1,-1,1,0,0]/1.11111*np.sqrt(c)/10
C_11_0[jj]=ChomHH[j-1,-1,0,0,0]/1.11111*np.sqrt(c)/10
C_11_EF[jj]=ChomEF[j-1,0,0]/1.11111*np.sqrt(c)/10
C_22[jj]=ChomHH[j-1,-1,1,1,1]/1.11111*np.sqrt(c)/10
C_22_0[jj]=ChomHH[j-1,-1,0,1,1]/1.11111*np.sqrt(c)/10
C_22_EF[jj]=ChomEF[j-1,1,1]/1.11111*np.sqrt(c)/10
C_12[jj]=ChomHH[j-1,-1,1,0,1]/0.277777*np.sqrt(c)/10
C_12_0[jj]=ChomHH[j-1,-1,0,0,1]/0.277777*np.sqrt(c)/10
C_12_EF[jj]=ChomEF[j-1,0,1]/0.277777*np.sqrt(c)/10
C_33[jj]=ChomHH[j-1,-1,1,2,2]/2/0.416666*np.sqrt(c)/10
C_33_0[jj]=ChomHH[j-1,-1,0,2,2]/2/0.416666*np.sqrt(c)/10
C_33_EF[jj]=ChomEF[j-1,2,2]/2/0.416666*np.sqrt(c)/10
plt.figure()
g1,=plt.plot(fractions,C_11)
g2,=plt.plot(fractions,C_11_0)
g11,=plt.plot(fractions,C_11_EF)
plt.legend([g1,g2,g11],['HH','HH C0','EF'])
plt.title('C_11 contraste '+str(c)+' échantillon '+str(i))
plt.xlabel('f')
plt.xlabel('f_souple')
plt.ylabel('Cij/Cij_rigide')
plt.axis([0,1,0,1])
plt.show()
plt.savefig('courbes/C_11 contraste '+str(c)+' échantillon '+str(i)+'.pdf')
#plt.show()
plt.close()
plt.figure()
......@@ -144,10 +169,11 @@ def courbes(c,i):
g11,=plt.plot(fractions,C_22_EF)
plt.legend([g1,g2,g11],['HH','HH C0','EF'])
plt.title('C_22 contraste '+str(c)+' échantillon '+str(i))
plt.xlabel('f')
plt.xlabel('f_souple')
plt.ylabel('Cij/Cij_rigide')
plt.axis([0,1,0,1])
plt.show()
plt.savefig('courbes/C_22 contraste '+str(c)+' échantillon '+str(i)+'.pdf')
#plt.show()
plt.close()
plt.figure()
......@@ -156,10 +182,11 @@ def courbes(c,i):
g11,=plt.plot(fractions,C_12_EF)
plt.legend([g1,g2,g11],['HH','HH C0','EF'])
plt.title('C_12 contraste '+str(c)+' échantillon '+str(i))
plt.xlabel('f')
plt.xlabel('f_souple')
plt.ylabel('Cij/Cij_rigide')
plt.axis([0,1,0,1])
plt.show()
plt.savefig('courbes/C_12 contraste '+str(c)+' échantillon '+str(i)+'.pdf')
#plt.show()
plt.close()
plt.figure()
......@@ -168,10 +195,11 @@ def courbes(c,i):
g11,=plt.plot(fractions,C_33_EF)
plt.legend([g1,g2,g11],['HH','HH C0','EF'])
plt.title('C_33 contraste '+str(c)+' échantillon '+str(i))
plt.xlabel('f')
plt.xlabel('f_souple')
plt.ylabel('Cij/Cij_rigide')
plt.axis([0,1,0,1])
plt.show()
plt.savefig('courbes/C_33 contraste '+str(c)+' échantillon '+str(i)+'.pdf')
#plt.show()
plt.close()
......@@ -181,34 +209,54 @@ def courbe_moyenne(c):
fi = h5py.File("c="+str(c)+"/c="+str(c)+"ech"+str(i)+".hdf5","r")
fichiers.append(fi)
C_11=np.zeros((8,))
C_11_EF=np.zeros((8,))
C_22=np.zeros((8,))
C_22_EF=np.zeros((8,))
C_12=np.zeros((8,))
C_12_EF=np.zeros((8,))
C_33=np.zeros((8,))
C_33_EF=np.zeros((8,))
C_11_0=np.zeros((8,))
C_22_0=np.zeros((8,))
C_12_0=np.zeros((8,))
C_33_0=np.zeros((8,))
C_11=np.ones((12,))
C_11_EF=np.ones((12,))
C_22=np.ones((12,))
C_22_EF=np.ones((12,))
C_12=np.ones((12,))
C_12_EF=np.ones((12,))
C_33=np.ones((12,))
C_33_EF=np.ones((12,))
C_11_0=np.ones((12,))
C_22_0=np.ones((12,))
C_12_0=np.ones((12,))
C_33_0=np.ones((12,))
fractions=np.ones((12,))
for i in range(10):
ChomHH=fichiers[i]['ChomHH']
ChomEF=fichiers[i]['ChomEF']
for j,f in enumerate([0.05,0.1,0.2,0.3,0.4,0.5,0.6,0.8]):
C_11[j]+=ChomHH[j,-1,1,0,0]/1.11111*np.sqrt(c)/10/10
C_11_0[j]+=ChomHH[j,-1,0,0,0]/1.11111*np.sqrt(c)/10/10
C_11_EF[j]+=ChomEF[j,0,0]/1.11111*np.sqrt(c)/10/10
C_22[j]+=ChomHH[j,-1,1,1,1]/1.11111*np.sqrt(c)/10/10
C_22_0[j]+=ChomHH[j,-1,0,1,1]/1.11111*np.sqrt(c)/10/10
C_22_EF[j]+=ChomEF[j,1,1]/1.11111*np.sqrt(c)/10/10
C_12[j]+=ChomHH[j,-1,1,0,1]/0.277777*np.sqrt(c)/10/10
C_12_0[j]+=ChomHH[j,-1,0,0,1]/0.277777*np.sqrt(c)/10/10
C_12_EF[j]+=ChomEF[j,0,1]/0.277777*np.sqrt(c)/10/10
C_33[j]+=ChomHH[j,-1,1,2,2]/2/0.416666*np.sqrt(c)/10/10
C_33_0[j]+=ChomHH[j,-1,0,2,2]/2/0.416666*np.sqrt(c)/10/10
C_33_EF[j]+=ChomEF[j,2,2]/2/0.416666*np.sqrt(c)/10/10
for j,f in enumerate([0.,0.05,0.1,0.2,0.3,0.4,0.5,0.6,0.67,0.8,0.9,1.]):
if c==0.01:
jj=11-j
else:
jj=j
fractions[jj]=f
if jj==11:
C_11[jj]=0.01
C_11_0[jj]=0.01
C_11_EF[jj]=0.01
C_22[jj]=0.01
C_22_0[jj]=0.01
C_22_EF[jj]=0.01
C_12[jj]=0.01
C_12_0[jj]=0.01
C_12_EF[jj]=0.01
C_33[jj]=0.01
C_33_0[jj]=0.01
C_33_EF[jj]=0.01
if f!=0. and f!=1.:
C_11[jj]+=ChomHH[j-1,-1,1,0,0]/1.11111*np.sqrt(c)/10/10
C_11_0[jj]+=ChomHH[j-1,-1,0,0,0]/1.11111*np.sqrt(c)/10/10
C_11_EF[jj]+=ChomEF[j-1,0,0]/1.11111*np.sqrt(c)/10/10
C_22[jj]+=ChomHH[j-1,-1,1,1,1]/1.11111*np.sqrt(c)/10/10
C_22_0[jj]+=ChomHH[j-1,-1,0,1,1]/1.11111*np.sqrt(c)/10/10
C_22_EF[jj]+=ChomEF[j-1,1,1]/1.11111*np.sqrt(c)/10/10
C_12[jj]+=ChomHH[j-1,-1,1,0,1]/0.277777*np.sqrt(c)/10/10
C_12_0[jj]+=ChomHH[j-1,-1,0,0,1]/0.277777*np.sqrt(c)/10/10
C_12_EF[jj]+=ChomEF[j-1,0,1]/0.277777*np.sqrt(c)/10/10
C_33[jj]+=ChomHH[j-1,-1,1,2,2]/2/0.416666*np.sqrt(c)/10/10
C_33_0[jj]+=ChomHH[j-1,-1,0,2,2]/2/0.416666*np.sqrt(c)/10/10
C_33_EF[jj]+=ChomEF[j-1,2,2]/2/0.416666*np.sqrt(c)/10/10
plt.figure()
g1,=plt.plot(fractions,C_11)
......@@ -216,10 +264,11 @@ def courbe_moyenne(c):
g11,=plt.plot(fractions,C_11_EF)
plt.legend([g1,g2,g11],['HH','HH C0','EF'])
plt.title('C_11 contraste '+str(c))
plt.xlabel('f')
plt.xlabel('f_souple')
plt.ylabel('Cij/Cij_rigide')
plt.axis([0,1,0,1])
plt.show()
plt.savefig('courbes/C_11 moyen contraste '+str(c)+'.pdf')
#plt.show()
plt.close()
plt.figure()
......@@ -228,10 +277,11 @@ def courbe_moyenne(c):
g11,=plt.plot(fractions,C_22_EF)
plt.legend([g1,g2,g11],['HH','HH C0','EF'])
plt.title('C_22 contraste '+str(c))
plt.xlabel('f')
plt.xlabel('f_souple')
plt.ylabel('Cij/Cij_rigide')
plt.axis([0,1,0,1])
plt.show()
plt.savefig('courbes/C_22 moyen contraste '+str(c)+'.pdf')
#plt.show()
plt.close()
plt.figure()
......@@ -240,10 +290,11 @@ def courbe_moyenne(c):
g11,=plt.plot(fractions,C_12_EF)
plt.legend([g1,g2,g11],['HH','HH C0','EF'])
plt.title('C_12 contraste '+str(c))
plt.xlabel('f')
plt.xlabel('f_souple')
plt.ylabel('Cij/Cij_rigide')
plt.axis([0,1,0,1])
plt.show()
plt.savefig('courbes/C_12 moyen contraste '+str(c)+'.pdf')
#plt.show()
plt.close()
plt.figure()
......@@ -252,13 +303,79 @@ def courbe_moyenne(c):
g11,=plt.plot(fractions,C_33_EF)
plt.legend([g1,g2,g11],['HH','HH C0','EF'])
plt.title('C_33 contraste '+str(c))
plt.xlabel('f')
plt.xlabel('f_souple')
plt.ylabel('Cij/Cij_rigide')
plt.axis([0,1,0,1])
plt.show()
plt.savefig('courbes/C_33 moyen contraste '+str(c)+'.pdf')
#plt.show()
plt.close()
def erreur_abs_def(c,i,j):
fi = h5py.File("c="+str(c)+"/c="+str(c)+"ech"+str(i)+".hdf5","r")
tAHH=fi['EpsHH'][j,9]
tAEF=fi['EpsEF'][j]
frac=[0.05,0.1,0.2,0.3,0.4,0.5,0.6,0.67,0.8,0.9]
f=frac[j]
plt.figure()
pp=plt.imshow(tAHH[:,:,0,0]-tAEF[:,:,0,0],origin='lower')
plt.colorbar(pp)
plt.savefig('c='+str(c)+'/ech='+str(i)+'/frac='+str(f)+'/erreur abs eps 11.pdf')
#plt.show()
plt.close()
plt.figure()
pp=plt.imshow(tAHH[:,:,0,1]-tAEF[:,:,0,1],origin='lower')
plt.colorbar(pp)
plt.savefig('c='+str(c)+'/ech='+str(i)+'/frac='+str(f)+'/erreur abs eps 12.pdf')
#plt.show()
plt.close()
plt.figure()
pp=plt.imshow(tAHH[:,:,1,1]-tAEF[:,:,1,1],origin='lower')
plt.colorbar(pp)
plt.savefig('c='+str(c)+'/ech='+str(i)+'/frac='+str(f)+'/erreur abs eps 22.pdf')
#plt.show()
plt.close()
plt.figure()
pp=plt.imshow(tAHH[:,:,2,2]-tAEF[:,:,2,2],origin='lower')
plt.colorbar(pp)
plt.savefig('c='+str(c)+'/ech='+str(i)+'/frac='+str(f)+'/erreur abs eps 33.pdf')
#plt.show()
plt.close()
def erreur_rel_def(c,i,j):
fi = h5py.File("c="+str(c)+"/c="+str(c)+"ech"+str(i)+".hdf5","r")
tAHH=fi['EpsHH'][j,9]
tAEF=fi['EpsEF'][j]
ii,jj,kk,ll=tAHH.shape
err=np.zeros((ii,jj,kk,ll))
for k in range(ii):
for l in range(jj):
for m in range(kk):
for n in range(ll):
err[k,l,m,n]=tAHH[k,l,m,n]/tAEF[k,l,m,n]-1
frac=[0.05,0.1,0.2,0.3,0.4,0.5,0.6,0.67,0.8,0.9]
f=frac[j]
plt.figure()
pp=plt.imshow(err[:,:,0,0],origin='lower')
plt.colorbar(pp)
plt.savefig('c='+str(c)+'/ech='+str(i)+'/frac='+str(f)+'/erreur rel eps 11.pdf')
#plt.show()
plt.close()
plt.figure()
pp=plt.imshow(err[:,:,0,1],origin='lower')
plt.colorbar(pp)
plt.savefig('c='+str(c)+'/ech='+str(i)+'/frac='+str(f)+'/erreur rel eps 12.pdf')
#plt.show()
plt.close()
plt.figure()
pp=plt.imshow(err[:,:,1,1],origin='lower')
plt.colorbar(pp)
plt.savefig('c='+str(c)+'/ech='+str(i)+'/frac='+str(f)+'/erreur rel eps 22.pdf')
#plt.show()
plt.close()
plt.figure()
pp=plt.imshow(err[:,:,2,2],origin='lower')
plt.colorbar(pp)
plt.savefig('c='+str(c)+'/ech='+str(i)+'/frac='+str(f)+'/erreur rel eps 33.pdf')
#plt.show()
plt.close()
\ No newline at end of file
......@@ -40,29 +40,31 @@ 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.9)
im,centers_diam,f1=image_init(N,D,0.9)
l=len(centers_diam)
fi['centers']=np.zeros((8,l,3))
fi['images']=np.zeros((8,N,N))
fi['fractions']=np.zeros((8,))
fi['ChomEF']=np.zeros((8,3,3))
fi['CVoigt']=np.zeros((8,3,3))
fi['ChomHH']=np.zeros((8,10,2,3,3))
fi['EpsHH']=np.zeros((8,10,N,N,3,3))
for j,f in enumerate([0.8,0.6,0.5,0.4,0.3,0.2,0.1,0.05]):
fi['centers']=np.zeros((10,l,3))
fi['images']=np.zeros((10,N,N))
fi['fractions']=np.zeros((10,))
fi['ChomEF']=np.zeros((10,3,3))
fi['CVoigt']=np.zeros((10,3,3))
fi['ChomHH']=np.zeros((10,10,2,3,3))
fi['EpsHH']=np.zeros((10,10,N,N,3,3))
fi['EpsEF']=np.zeros((10,N,N,3,3))
for j,f in enumerate([0.9,0.8,0.67,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')
f1=image_diminue(im,centers_diam,f)
if f!=0.9:
f1=image_diminue(im,centers_diam,f)
# plt.figure()
# plt.imshow(im[:,:],origin='lower')
# plt.show()
# plt.close()
fi['images'][7-j]=im
fi['images'][9-j]=im
l=len(centers_diam)
for ii in range(l):
fi['centers'][7-j][ii]=centers_diam[ii]
fi['fractions'][7-j]=f1
fi['centers'][9-j][ii]=centers_diam[ii]
fi['fractions'][9-j]=f1
tC=np.array([[C[0] for ii in range(N)] for jj in range(N)])
for ii in range(N):
......@@ -70,23 +72,24 @@ for c in [100,0.01]:
if im[ii,jj]==1:
tC[ii,jj]=C[1]
ChomEF,CVoigt = reference(tC,256,c,i,f)
fi['ChomEF'][7-j]=ChomEF
fi['CVoigt'][7-j]=CVoigt
ChomEF,CVoigt,tA = reference(tC,256,c,i,f)
fi['ChomEF'][9-j]=ChomEF
fi['CVoigt'][9-j]=CVoigt
fi['EpsEF'][9-j]=tA
ChomHH_list,tA_list=homog_hierar(tC,CVoigt)
l=len(ChomHH_list)
if l<10:
for k in range(l):
fi['ChomHH'][7-j,k]=ChomHH_list[k]
fi['EpsHH'][7-j,k]=tA_list[k]
fi['ChomHH'][9-j,k]=ChomHH_list[k]
fi['EpsHH'][9-j,k]=tA_list[k]
for k in range(l,10):
fi['ChomHH'][7-j,k]=ChomHH_list[l-1]
fi['EpsHH'][7-j,k]=tA_list[l-1]
fi['ChomHH'][9-j,k]=ChomHH_list[l-1]
fi['EpsHH'][9-j,k]=tA_list[l-1]
else:
for k in range(10):
fi['ChomHH'][7-j,k]=ChomHH_list[k]
fi['EpsHH'][7-j,k]=tA_list[k]
fi['ChomHH'][9-j,k]=ChomHH_list[k]
fi['EpsHH'][9-j,k]=tA_list[k]
......
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