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) –

• triangular matrix (A) –

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

Returns

• alpha (numpy.ndarray)

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