Estoy haciendo este código y cuando hago print()
de alguna de las variables: el, et o B
, me devuelve: none
. ¿Qué puede estar pasando? ¿Puede ser solo un error de indentación o hay algo más?
from numba import jit
import numpy as np
import random
random.seed(7)
D=2
N=32
Zc=0.25
kl=1e7
B=np.zeros((N+2,N+2))
@jit(nopython=True)
def M(B):
B=np.zeros((N+2,N+2))
C=np.zeros((N+2,N+2))
s=2*D+1
el=[]
et=[]
for i in range(kl):
e=0
for j in range(1,N+1):
for k in range(1,N+1):
Z=B[j,k]-((1/2*D))*(B[j+1,k]+B[j-1,k]+B[j,k+1]+B[j,k-1])
if Z > Zc:
C[j,k]= C[j,k]-((2*D)/s)*Zc
C[j+1,k]= C[j+1,k]+(1/s)*Zc
C[j-1,k]= C[j-1,k]+(1/s)*Zc
C[j,k+1]= C[j,k+1]+(1/s)*Zc
C[j,k-1]= C[j,k-1]+(1/s)*Zc
g=((2*D)/s)*((2*abs(Z)/Zc)-1)*(Zc**2)
e=e+g
else:
continue
if e>0:
for j in range(1,N+1):
for k in range(1,N+1):
B[j,k]= B[j,k] + C[j,k]
C[j,k]=0
else:
k0=[random.randint(1,N+1),random.randint(1,N+1)]
deltaB=random.random()-0.2
B[k0[0],k0[1]]= B[k0[0],k0[1]]+deltaB
B[0,:]=0
B[:,0]=0
B[N+1,:]=0
B[:,N+1]=0
el.append(e)
et.append(np.sum(np.square(B)))
return el ,et,B
B=M(B)
el=M(B)
et=M(B)
Puede que el problema esté en la indentación del return
, si no esta ahí el problema, no sé en que otro lugar puede ser el error. Acá esta el código completo y corregido el problema de las variables que no estaban definidas.