def test_cache_refresh_D_rates(self):
'''Functionality test for refresh cache behaviour.'''
Gamma_L = 1.
Gamma_R = 1.
srl_solver = utils.setup_srl_solver_from_hilbert_space(Gamma_L, Gamma_R)
expected_F2 = utils.zero_freq_F2_srl(Gamma_L, Gamma_R)
self.assertAlmostEqual(expected_F2, srl_solver.second_order_fano_factor(0))
Gamma_R = 0.5
srl_solver.D_rates = np.array([Gamma_R, Gamma_L])
expected_F2 = utils.zero_freq_F2_srl(Gamma_L, Gamma_R)
self.assertAlmostEqual(expected_F2, srl_solver.second_order_fano_factor(0))
def test_cache_refresh_D_rates(self):
'''Functionality test for refresh cache behaviour.'''
Gamma_L = 1.
Gamma_R = 1.
srl_solver = utils.setup_srl_solver_from_hilbert_space(
Gamma_L, Gamma_R)
expected_F2 = utils.zero_freq_F2_srl(Gamma_L, Gamma_R)
self.assertAlmostEqual(expected_F2,
srl_solver.second_order_fano_factor(0))
Gamma_R = 0.5
srl_solver.D_rates = np.array([Gamma_R, Gamma_L])
expected_F2 = utils.zero_freq_F2_srl(Gamma_L, Gamma_R)
self.assertAlmostEqual(expected_F2,
srl_solver.second_order_fano_factor(0))
def test_generate_cumulant_srl(self):
Gamma_R = 0.5; Gamma_L = 1.
expected_mean = utils.mean_srl(Gamma_L, Gamma_R)
expected_F2 = utils.zero_freq_F2_srl(Gamma_L, Gamma_R)
expected_F3 = utils.skewness_srl(0, 0, Gamma_L, Gamma_R)
srl_solver = utils.setup_srl_solver_from_hilbert_space(Gamma_L, Gamma_R)
self.assertAlmostEqual(expected_mean, srl_solver.generate_cumulant(1)[0][0])
self.assertAlmostEqual(expected_F2, srl_solver.generate_cumulant(2)[0][1])
self.assertAlmostEqual(expected_F3, srl_solver.generate_cumulant(3)[0][2])
def test_zero_freq_F2_srl(self):
'''Functionality test for second order Fano factor of SRL model with single frequency.'''
Gamma_R = 0.5
Gamma_L = 1.
expected_F2 = utils.zero_freq_F2_srl(Gamma_L, Gamma_R)
srl_solver = utils.setup_srl_solver_from_hilbert_space(
Gamma_L, Gamma_R)
self.assertAlmostEqual(expected_F2,
srl_solver.second_order_fano_factor(0))
def test_generate_cumulant_srl(self):
Gamma_R = 0.5
Gamma_L = 1.
expected_mean = utils.mean_srl(Gamma_L, Gamma_R)
expected_F2 = utils.zero_freq_F2_srl(Gamma_L, Gamma_R)
expected_F3 = utils.skewness_srl(0, 0, Gamma_L, Gamma_R)
srl_solver = utils.setup_srl_solver_from_hilbert_space(
Gamma_L, Gamma_R)
self.assertAlmostEqual(expected_mean,
srl_solver.generate_cumulant(1)[0][0])
self.assertAlmostEqual(expected_F2,
srl_solver.generate_cumulant(2)[0][1])
self.assertAlmostEqual(expected_F3,
srl_solver.generate_cumulant(3)[0][2])
def test_zero_freq_F2_srl(self):
'''Functionality test for second order Fano factor of SRL model with single frequency.'''
Gamma_R = 0.5; Gamma_L = 1.
expected_F2 = utils.zero_freq_F2_srl(Gamma_L, Gamma_R)
srl_solver = utils.setup_srl_solver_from_hilbert_space(Gamma_L, Gamma_R)
self.assertAlmostEqual(expected_F2, srl_solver.second_order_fano_factor(0))
