physbo.misc.gauss_elim module¶
-
physbo.misc.gauss_elim.
gauss_elim
(U, t)[source]¶ 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.