def eliminacion_gaussiana_pivoteo_total(self, A, b):
"""
Reliza eliminacion gaussiana con pivoteo total,
dado A y b
"""
eg = eliminacion_gaussiana.EliminacionGaussiana()
eu = eliminacion_gaussiana.GaussUtilities()
(U, b_a) = eg.escalonar(A, b, pivot=2, u=True, b_a=True)
x = eu.solucionar_regr(U, b_a, marcs=eg.marcas())
return x
def eliminacion_gaussiana_simple(self, A, b):
"""
Reliza eliminacion gaussiana simple,
dado A y b
"""
eg = eliminacion_gaussiana.EliminacionGaussiana()
eu = eliminacion_gaussiana.GaussUtilities()
(U, b_a) = eg.escalonar(A, b, u=True, b_a=True)
x = eu.solucionar_regr(U, b_a)
return x
def factorizacion_LU_gauss_pivoteo_parcial(A, b):
"""
Reliza la factorizacion A = LU,
Lb = z, Ux = b y devuelve x,
dado A y b
"""
eg = eliminacion_gaussiana.EliminacionGaussiana()
eu = eliminacion_gaussiana.GaussUtilities()
(L, U) = eg.escalonar(A, b, pivot=1, l=True, u=True)
z = eu.solucionar_prog(L, b)
x = eu.solucionar_regr(U, z)
return x
def factorizacion_LU_gauss_simple(A, b):
"""
Reliza la factorizacion A = LU,
Lb = z, Ux = b y devuelve x,
dado A y b
"""
eg = eliminacion_gaussiana.EliminacionGaussiana()
eu = eliminacion_gaussiana.GaussUtilities()
(L, U) = eg.escalonar(A, b, pivot=0, l=True, u=True)
print("L= ", L)
print("U= ", U)
z = eu.solucionar_prog(L, b)
x = eu.solucionar_regr(U, z)
return x