def test_generate_cumulant_dqd(self):
Gamma_L = 1.; Gamma_R = 0.5; Tc = 3.; bias = 1.
expected_mean = utils.mean_dqd(Gamma_L, Gamma_R, Tc, bias)
expected_F2 = utils.zero_freq_F2_dqd(Gamma_L, Gamma_R, Tc, bias)
dqd_solver = utils.setup_dqd_solver_from_hilbert_space(Gamma_L, Gamma_R, Tc, bias)
self.assertAlmostEqual(expected_mean, dqd_solver.generate_cumulant(1)[0][0])
self.assertAlmostEqual(expected_F2, dqd_solver.generate_cumulant(2)[0][1])
def test_zero_freq_F2_dqd(self):
'''Functionality test for second order Fano factor of DQD model with single frequency.'''
Gamma_L = 1.
Gamma_R = 0.5
Tc = 3.
bias = 1.
expected_F2 = utils.zero_freq_F2_dqd(Gamma_L, Gamma_R, Tc, bias)
dqd_solver = utils.setup_dqd_solver_from_hilbert_space(
Gamma_L, Gamma_R, Tc, bias)
self.assertAlmostEqual(expected_F2,
dqd_solver.second_order_fano_factor(0))
def test_generate_cumulant_dqd(self):
Gamma_L = 1.
Gamma_R = 0.5
Tc = 3.
bias = 1.
expected_mean = utils.mean_dqd(Gamma_L, Gamma_R, Tc, bias)
expected_F2 = utils.zero_freq_F2_dqd(Gamma_L, Gamma_R, Tc, bias)
dqd_solver = utils.setup_dqd_solver_from_hilbert_space(
Gamma_L, Gamma_R, Tc, bias)
self.assertAlmostEqual(expected_mean,
dqd_solver.generate_cumulant(1)[0][0])
cumulants, states = dqd_solver.generate_cumulant(2)
self.assertAlmostEqual(expected_mean * expected_F2, cumulants[1])
def test_zero_freq_F2_dqd(self):
'''Functionality test for second order Fano factor of DQD model with single frequency.'''
Gamma_L = 1.; Gamma_R = 0.5; Tc = 3.; bias = 1.
expected_F2 = utils.zero_freq_F2_dqd(Gamma_L, Gamma_R, Tc, bias)
dqd_solver = utils.setup_dqd_solver_from_hilbert_space(Gamma_L, Gamma_R, Tc, bias)
self.assertAlmostEqual(expected_F2, dqd_solver.second_order_fano_factor(0))