Buenas, ando liado con una matriz que se crea con un bucle for múltiple que a su vez dentro del bucle se integra a través de la función quad el tamaño de la matriz depende de unos valores de m y n. Para unos valores de m y n de 8 tarda aproximadamente unos 180 segundos (3 minutos) en crear la matriz. Algún consejo para optimizar el llenado de la matriz.
Un saludo.
M=N=8
def I00(x,i,m,a):
I00=X(x,i,a)*X(x,m,a)
return I00
K=np.zeros(shape=(m*n,m*n))
M=np.zeros(shape=(m*n,m*n))
mmm=np.linspace(0,m,num=m,endpoint=False,dtype='int')
ii=np.linspace(0,m,num=m,endpoint=False,dtype='int')
nnn=np.linspace(0,n,num=n,endpoint=False,dtype='int')
jj=np.linspace(0,n,num=n,endpoint=False,dtype='int')
mii=np.linspace(0,n*m,num=m,endpoint=False,dtype='int')
njj=np.linspace(0,n*m,num=n,endpoint=False,dtype='int')
for i,mi in zip(ii,mii):
for j in jj:
for mm,nj in zip(mmm,njj):
for nn in nnn:
K[mi+j][nj+nn]=(D[0][0]*integrate.quad(I22,0,a,args=(mm+1,i+1,a))[0]
*integrate.quad(I00,0,b,args=(nn+1,j+1,b))[0]) +
(D[1][1]*integrate.quad(I00,0,a,args=(mm+1,i+1,a))[0]
*integrate.quad(I22,0,b,args=(nn+1,j+1,b))[0]) +
(D[0][1]*((integrate.quad(I20,0,a,args=(mm+1,i+1,a))[0]
*integrate.quad(I02,0,b,args=(nn+1,j+1,b))[0]) +
(integrate.quad(I02,0,a,args=(mm+1,i+1,a))[0]
*integrate.quad(I20,0,b,args=(nn+1,j+1,b))[0]))) +
(4*D[2][2]*integrate.quad(I11,0,a,args=(mm+1,i+1,a))[0]
*integrate.quad(I11,0,b,args=(nn+1,j+1,b))[0]) +
(2*D[0][2]*((integrate.quad(I21,0,a,args=(mm+1,i+1,a))[0]
*integrate.quad(I01,0,b,args=(nn+1,j+1,b))[0]) +
(integrate.quad(I12,0,a,args=(mm+1,i+1,a))[0]
*integrate.quad(I10,0,b,args=(nn+1,j+1,b))[0]))) +
(2*D[1][2]*((integrate.quad(I10,0,a,args=(mm+1,i+1,a))[0]
*integrate.quad(I12,0,b,args=(nn+1,j+1,b))[0]) +
(integrate.quad(I01,0,a,args=(mm+1,i+1,a))[0]
*integrate.quad(I21,0,b,args=(nn+1,j+1,b))[0])))
M[mi+j][nj+nn]=rho*h*integrate.quad(I00, 0, a, args=(mm+1,i+1,a))[0]
*integrate.quad(I00, 0, b, args=(nn+1,j+1,b))[0]