coef.alloc(nvx)
#fvm.Coefficients.alloc(nvx)
#
# Calculamos los coeficientes de FVM de la Difusión
#
dif = fvm.Diffusion1D(nvx, Gamma = Gamma, dx = delta)
DE, DW, DP = dif.calcCoef()
#print('aW = {}'.format(dif.aW()),
# 'aE = {}'.format(dif.aE()),
# 'Su = {}'.format(dif.Su()),
# 'aP = {}'.format(dif.aP()), sep='\n')
#print('.'+'-'*70+'.')
#
# Calculamos los coeficientes de FVM de la Advección
#
adv = fvm.Advection1D(nvx, rho = rho, dx = delta)
adv.setU(u)
adv.calcCoef('Quick', phiA, phiB)
print('aW = {}'.format(dif.aW()),
'aE = {}'.format(dif.aE()),
'Su = {}'.format(dif.Su()),
'aP = {}'.format(dif.aP()), sep='\n')
#print('u = {}'.format(adv.u()))
#print('.'+'-'*70+'.')
#
# Se construye el arreglo donde se guardará la solución
#
phi = np.zeros(nvx) # El arreglo contiene ceros
phi[0] = phiA # Condición de frontera izquierda
phi[-1] = phiB # Condición de frontera derecha
#
#
# Se construye el arreglo donde se guardará la solución
#
Phi = np.zeros(nvx) # El arreglo contiene ceros
Phi[0] = PhiA # Condición de frontera izquierda
Phi[-1] = PhiB # Condición de frontera derecha
#
# Creamos los coeficientes de FVM
#
coef = fvm.Coefficients()
coef.alloc(nvx)
df1 = fvm.Diffusion1D(nvx, Gamma=gamma, dx=dx)
adv1 = fvm.Advection1D(nvx, rho=rho, dx=dx)
adv1.setU(u)
D = gamma / dx
tem = fvm.Temporal1D(nvx, rho=rho, dx=dx, dt=dt)
# Calculamos la solución analítica
#
Phiaf = analyticSol(x, u, Tf - Ti, gamma)
vfl.grafica(x, Phiaf, kind='-', label='solución analítica')
for i in range(1, int((Tf - Ti) / dt) + 1):
#
# Iniciamos con coeficientes cero
#