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

  • 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.