projects/nnum/decomposition/solvers.py
def forward_elimination(l, y):
    """Perform forward elimination over a Lz = y linear system.

    :param l: a lower triangular matrix retrieved from a method such as LU Decomposition.
    :param y: the response column-vector of the linear system.
    :return: z, where y = Lz.
    """
    z = y.copy().astype(float)

    for k in range(y.shape[0]):
        z[k] -= np.sum(z[:k] * l[k, :k])
        z[k] /= l[k, k]

    return z
Perform forward elimination over a Lz = y linear system.

    :param l: a lower triangular matrix retrieved from a method such as LU Decomposition.
    :param y: the response column-vector of the linear system.
    :return: z, where y = Lz.
    
Perform forward elimination over a Lz = y linear system.
Perform forward elimination over a Lz = y linear system .
def forward_elimination(l, y):
    z = y.copy().astype(float)
    for k in range(y.shape[0]):
        z[k] -= np.sum(z[:k] * l[k, :k])
        z[k] /= l[k, k]
    return z
def forward_elimination ( l , y ) : z = y . copy ( ) . astype ( float ) for k in range ( y . shape [ 0 ] ) : z [ k ] -= np . sum ( z [ : k ] * l [ k , : k ] ) z [ k ] /= l [ k , k ] return z
Perform forward elimination over a Lz = y linear system.
Perform forward elimination over a Lz = y linear system .
"""Perform forward elimination over a Lz = y linear system.

    :param l: a lower triangular matrix retrieved from a method such as LU Decomposition.
    :param y: the response column-vector of the linear system.
    :return: z, where y = Lz.
    """
z, where y = Lz.
z where y = Lz .
l
a lower triangular matrix retrieved from a method such as LU Decomposition.
a lower triangular matrix retrieved from a method such as LU Decomposition .
y
the response column-vector of the linear system.
the response column - vector of the linear system .