def printFrame(d):
"""
Función que imprime el error calculado
"""
# Calculo el error porcentual y agrego al DataFrame
# una columna con esos datos llamada 'Error %'
d['Error %'] = d['Error'] / d['Analytic']
print(DataFrame(d))
print('.'+ '-'*70 + '.')
def calcError(phiA, phiN):
"""
Función que calcula el error entre la aproximación calculada y el valor real de la solución.
@return: error
"""
return np.absolute(phiA - phiN)
if __name__ == '__main__':
Coefficients.alloc(5)
m = Mesh(nodes = 5)
d = Diffusion1D(m.volumes())
ma = Matrix(m.volumes())
a = Advection1D(m.volumes())
print(m.delta(), d.aP(), a.aP(), ma.mat(), sep='\n')
printData(nvx =5, nx = 6, longitud = 1.3)
