def test_get_qc_noisy_qpu_error(client_configuration: QCSClientConfiguration, dummy_compiler: DummyCompiler):
expected_message = (
"pyQuil currently does not support initializing a noisy QuantumComputer "
"based on a QCSQuantumProcessor. Change noisy to False or specify the name of a QVM."
)
with pytest.raises(ValueError, match=expected_message):
get_qc("Aspen-8", noisy=True)
def test_qc_error(client_configuration: QCSClientConfiguration):
# QVM is not a QPU
with pytest.raises(ValueError):
get_qc("9q-square-noisy-qvm", as_qvm=False, client_configuration=client_configuration)
with pytest.raises(ValueError):
get_qc("5q", as_qvm=False, client_configuration=client_configuration)
def test_qc_joint_calibration(client_configuration: QCSClientConfiguration):
# noise model with 95% symmetrized readout fidelity per qubit
noise_model = asymmetric_ro_model([0, 1], 0.945, 0.955)
qc = get_qc("2q-qvm", client_configuration=client_configuration)
qc.qam.noise_model = noise_model
# |01> state program
p = Program()
p += RESET()
p += X(0)
p.wrap_in_numshots_loop(10000)
# ZZ experiment
sz = ExperimentSetting(
in_state=_pauli_to_product_state(sZ(0) * sZ(1)), out_operator=sZ(0) * sZ(1), additional_expectations=[[0], [1]]
)
e = Experiment(settings=[sz], program=p)
results = qc.experiment(e)
# ZZ expectation value for state |01> with 95% RO fid on both qubits is about -0.81
assert np.isclose(results[0].expectation, -0.81, atol=0.01)
assert results[0].total_counts == 40000
# Z0 expectation value for state |01> with 95% RO fid on both qubits is about -0.9
assert np.isclose(results[0].additional_results[0].expectation, -0.9, atol=0.01)
assert results[0].additional_results[1].total_counts == 40000
# Z1 expectation value for state |01> with 95% RO fid on both qubits is about 0.9
assert np.isclose(results[0].additional_results[1].expectation, 0.9, atol=0.01)
assert results[0].additional_results[1].total_counts == 40000
def test_qc(client_configuration: QCSClientConfiguration):
qc = get_qc("9q-square-noisy-qvm", client_configuration=client_configuration)
assert isinstance(qc, QuantumComputer)
assert qc.qam.noise_model is not None
assert qc.qubit_topology().number_of_nodes() == 9
assert qc.qubit_topology().degree[0] == 2
assert qc.qubit_topology().degree[4] == 4
assert str(qc) == "9q-square-noisy-qvm"
def test_run_symmetrized_readout_error(client_configuration: QCSClientConfiguration):
# This test checks if the function runs for any possible input on a small number of qubits.
# Locally this test was run on all 8 qubits, but it was slow.
qc = get_qc("8q-qvm", client_configuration=client_configuration)
sym_type_vec = [-1, 0, 1, 2, 3]
prog_vec = [Program(I(x) for x in range(0, 3))[0:n] for n in range(0, 4)]
trials_vec = list(range(0, 5))
for prog, trials, sym_type in itertools.product(prog_vec, trials_vec, sym_type_vec):
print(qc.run_symmetrized_readout(prog, trials, sym_type))
def test_noisy(client_configuration: QCSClientConfiguration):
# https://github.com/rigetti/pyquil/issues/764
p = Program(
Declare("ro", "BIT", 1),
X(0),
MEASURE(0, ("ro", 0)),
).wrap_in_numshots_loop(10000)
qc = get_qc("1q-qvm", noisy=True, client_configuration=client_configuration)
result = qc.run(qc.compile(p))
assert result.mean() < 1.0
def test_qc_run(client_configuration: QCSClientConfiguration):
qc = get_qc("9q-square-noisy-qvm", client_configuration=client_configuration)
bs = qc.run(
qc.compile(
Program(
Declare("ro", "BIT", 1),
X(0),
MEASURE(0, ("ro", 0)),
).wrap_in_numshots_loop(3)
)
)
assert bs.shape == (3, 1)
def test_qc_calibration_2q(client_configuration: QCSClientConfiguration):
# noise model with 95% symmetrized readout fidelity per qubit
noise_model = asymmetric_ro_model([0, 1], 0.945, 0.955)
qc = get_qc("2q-qvm", client_configuration=client_configuration)
qc.qam.noise_model = noise_model
# bell state program (doesn't matter)
p = Program()
p += RESET()
p += H(0)
p += CNOT(0, 1)
p.wrap_in_numshots_loop(10000)
# ZZ experiment
sz = ExperimentSetting(in_state=_pauli_to_product_state(sZ(0) * sZ(1)), out_operator=sZ(0) * sZ(1))
e = Experiment(settings=[sz], program=p)
results = qc.calibrate(e)
# ZZ expectation should just be (1 - 2 * readout_error_q0) * (1 - 2 * readout_error_q1)
np.isclose(results[0].expectation, 0.81, atol=0.01)
assert results[0].total_counts == 40000
def test_qc_compile(dummy_compiler: DummyCompiler, client_configuration: QCSClientConfiguration):
qc = get_qc("5q", as_qvm=True, noisy=True, client_configuration=client_configuration)
qc.compiler = dummy_compiler
prog = Program()
prog += H(0)
assert qc.compile(prog) == prog
def test_qc_noisy(client_configuration: QCSClientConfiguration):
qc = get_qc("5q", as_qvm=True, noisy=True, client_configuration=client_configuration)
assert isinstance(qc, QuantumComputer)
def test_nq_qvm_qc(client_configuration: QCSClientConfiguration):
for n_qubits in [2, 4, 7, 19]:
qc = get_qc(f"{n_qubits}q-qvm", client_configuration=client_configuration)
for q1, q2 in itertools.permutations(range(n_qubits), r=2):
assert (q1, q2) in qc.qubit_topology().edges
assert qc.name == f"{n_qubits}q-qvm"