Пример #1
0
 def factorizacion_LU_cholesky(self, A, b):
     directo = lu_directo.LU_Directo()
     eu = eliminacion_gaussiana.GaussUtilities()
     (L, U) = directo.cholesky(A)
     z = eu.solucionar_prog(L, b)
     x = eu.solucionar_regr(U, z)
     return x
Пример #2
0
 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
Пример #3
0
 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 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
Пример #5
0
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
Пример #6
0
    def factorizacion_LU_cholesky(self, A, b):
        """
        Realiza la factorizacion A = LU,
        Lb = z, Ux = b y devuelve x,
        dado A y b.
        """

        directo = lu_directo.LU_Directo()
        eu = eliminacion_gaussiana.GaussUtilities()
        (L, U) = directo.cholesky(A)
        print(L)
        print(U)
        z = eu.solucionar_prog(L, b)
        x = eu.solucionar_regr(U, z)
        return x