Ejemplo n.º 1
0
 def test_thermal_state_precomp(self):
     full = rand_herm(2**4)
     beta = 0.624
     rhoth1 = thermal_state(full, beta)
     func = thermal_state(full, None, precomp_func=True)
     rhoth2 = func(beta)
     assert_allclose(rhoth1, rhoth2)
Ejemplo n.º 2
0
 def test_thermal_state_tuple(self):
     full = rand_herm(2**4)
     l, v = eigsys(full)
     for beta in (0, 0.5, 1, 10):
         rhoth1 = thermal_state(full, beta)
         rhoth2 = thermal_state((l, v), beta)
         assert_allclose(rhoth1, rhoth2)
Ejemplo n.º 3
0
 def test_thermal_state_normalization(self):
     full = rand_herm(2**4)
     for beta in (0, 0.5, 1, 10):
         rhoth = thermal_state(full, beta)
         assert_allclose(tr(rhoth), 1)
Ejemplo n.º 4
0
 def test_thermal_state_cold(self):
     full = ham_j1j2(4, j2=0.1253)
     rhoth = thermal_state(full, 100)
     gs = groundstate(full)
     assert_allclose(tr(gs.H @ rhoth @ gs), 1.0, rtol=1e-4)
Ejemplo n.º 5
0
 def test_thermal_state_hot(self):
     full = rand_herm(2**4)
     rhoth = chop(thermal_state(full, 0.0))
     assert_allclose(rhoth, eye(2**4) / 2**4)