def test_two_qubit_linear_inv(two_q_tomo_fixture):
qubits = [0, 1]
results, rho_true = two_q_tomo_fixture
rho_est = linear_inv_state_estimate(results, qubits)
np.testing.assert_allclose(rho_true, rho_est, atol=0.01)
def test_single_qubit_linear_inv(single_q_tomo_fixture):
qubits = [0]
results, rho_true = single_q_tomo_fixture
rho_est = linear_inv_state_estimate(results, qubits)
np.testing.assert_allclose(rho_true, rho_est, atol=1e-2)
