Method eigh

ot/backend.py:1420–1421

(self, a)

        return np.einsum("...jk,...kl->...jl", Q, np.swapaxes(V, -1, -2))

    def eigh(self, a):
        return np.linalg.eigh(a)

    def kl_div(self, p, q, mass=False, eps=1e-16):
        value = np.sum(p * np.log(p / q + eps))