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]

パラメータ:
  • U ((M, M) array_like)

  • matrix (A triangular)

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

戻り値:

  • alpha (numpy.ndarray)

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