def test_spectral_decomposition(self, tol):
"""Test that the correct spectral decomposition is returned."""
a, P = spectral_decomposition(H)
# verify that H = \sum_k a_k P_k
assert np.allclose(H, np.einsum("i,ijk->jk", a, P), atol=tol, rtol=0)
def test_spectral_decomposition(self):
"""Test that the correct spectral decomposition is returned."""
self.logTestName()
a, P = spectral_decomposition(H)
# verify that H = \sum_k a_k P_k
self.assertAllAlmostEqual(H, np.einsum("i,ijk->jk", a, P), delta=self.tol)
def test_spectral_decomposition_qubit(self, tol):
"""Test that the correct spectral decomposition is returned."""
a, P = spectral_decomposition_qubit(H)
# verify that H = \sum_k a_k P_k
self.assertAllAlmostEqual(H, np.einsum('i,ijk->jk', a, P), delta=tol)
