def setup_dqd_solver(Gamma_L, Gamma_R, Tc, bias):
L = reduced_dqd_liouvillian(bias, Tc, Gamma_L, Gamma_R)
jump_op = np.zeros((5, 5))
jump_op[0, 4] = 1
return FCSSolver(L, jump_op, np.array([1, 1, 1, 0, 0]))
def setup_srl_solver(Gamma_L, Gamma_R):
return FCSSolver(reduced_srl_liouvillian(Gamma_L, Gamma_R),
np.array([[0, 1], [0, 0]]), np.array([1, 1]))
