physbo.misc.gauss_elim module

physbo.misc.gauss_elim.gauss_elim(U, t)

Calculate alpha using scipy.linalg.solve_triangular. alpha = (U^T U)^-1 t = U^-1 [(U^T)-1 t]

Parameters:

• U ((M, M) array_like)

• matrix (A triangular)

• t ((M,) or (M, N) array_like)

Returns:

• alpha (numpy.ndarray)

• Solution to the system L^T alpha = t. Shape of return matches t.