def test_gate_product(self):
filename = "qft.qasm"
filepath = Path(__file__).parent / 'qasm_files' / filename
qc = read_qasm(filepath)
U_list_expanded = qc.propagators()
U_list = qc.propagators(expand=False)
inds_list = []
for gate in qc.gates:
if isinstance(gate, Measurement):
continue
else:
inds_list.append(gate.get_inds(qc.N))
U_1, _ = gate_sequence_product(U_list,
inds_list=inds_list,
expand=True)
U_2 = gate_sequence_product(U_list_expanded,
left_to_right=True,
expand=False)
np.testing.assert_allclose(U_1, U_2)
def test_wstate(self):
filename = "w-state.qasm"
filepath = Path(__file__).parent / 'qasm_files' / filename
qc = read_qasm(filepath)
rand_state = rand_ket(2)
state = tensor(tensor(basis(2, 0), basis(2, 0), basis(2, 0)),
rand_state)
fourth = Measurement("test_rand", targets=[3])
_, probs_initial = fourth.measurement_comp_basis(state)
simulators = _simulators_sv(qc)
for simulator in simulators:
result = simulator.run_statistics(state)
final_states = result.get_final_states()
result_cbits = result.get_cbits()
for i, final_state in enumerate(final_states):
_, probs_final = fourth.measurement_comp_basis(final_state)
np.testing.assert_allclose(probs_initial, probs_final)
assert sum(result_cbits[i]) == 1
def test_export_import():
qc = QubitCircuit(3)
qc.add_gate("CRY", targets=1, controls=0, arg_value=np.pi)
qc.add_gate("CRX", targets=1, controls=0, arg_value=np.pi)
qc.add_gate("CRZ", targets=1, controls=0, arg_value=np.pi)
qc.add_gate("CNOT", targets=1, controls=0)
qc.add_gate("TOFFOLI", targets=2, controls=[0, 1])
# qc.add_gate("SQRTNOT", targets=0)
qc.add_gate("CS", targets=1, controls=0)
qc.add_gate("CT", targets=1, controls=0)
qc.add_gate("SWAP", targets=[0, 1])
qc.add_gate("QASMU", targets=[0], arg_value=[np.pi, np.pi, np.pi])
qc.add_gate("RX", targets=[0], arg_value=np.pi)
qc.add_gate("RY", targets=[0], arg_value=np.pi)
qc.add_gate("RZ", targets=[0], arg_value=np.pi)
qc.add_gate("SNOT", targets=[0])
qc.add_gate("X", targets=[0])
qc.add_gate("Y", targets=[0])
qc.add_gate("Z", targets=[0])
qc.add_gate("S", targets=[0])
qc.add_gate("T", targets=[0])
# qc.add_gate("CSIGN", targets=[0], controls=[1])
read_qc = read_qasm(circuit_to_qasm_str(qc), strmode=True)
props = qc.propagators()
read_props = read_qc.propagators()
for u0, u1 in zip(props, read_props):
assert (u0 - u1).norm() < 1e-12
def test_custom_gates():
filename = "test_custom_gates.qasm"
filepath = Path(__file__).parent / 'qasm_files' / filename
qc = read_qasm(filepath)
unitaries = qc.propagators()
assert (unitaries[0] - unitaries[1]).norm() < 1e-12
ry_cx = cnot() * tensor(identity(2), ry(np.pi/2))
assert (unitaries[2] - ry_cx).norm() < 1e-12
def test_qasm_addcircuit():
filename = "test_add.qasm"
filepath = Path(__file__).parent / 'qasm_files' / filename
qc = read_qasm(filepath)
assert qc.N == 2
assert qc.num_cbits == 2
check_gate_defn(qc.gates[0], "X", [1])
check_gate_defn(qc.gates[1], "SNOT", [0])
check_gate_defn(qc.gates[2], "SNOT", [1])
check_gate_defn(qc.gates[3], "CNOT", [1], [0])
check_gate_defn(qc.gates[4], "SNOT", [0])
check_gate_defn(qc.gates[5], "SNOT", [1])
check_gate_defn(qc.gates[6], "SNOT", [0], None, [0, 1], 0)
check_measurement_defn(qc.gates[7], "M", [0], 0)
check_measurement_defn(qc.gates[8], "M", [1], 1)
def test_qasm_addcircuit():
filename = "test_add.qasm"
filepath = Path(__file__).parent / 'qasm_files' / filename
with pytest.warns(UserWarning, match="not preserved in QubitCircuit"):
qc = read_qasm(filepath)
assert qc.N == 2
assert qc.num_cbits == 2
check_gate_defn(qc.gates[0], "X", [1])
check_gate_defn(qc.gates[1], "SNOT", [0])
check_gate_defn(qc.gates[2], "SNOT", [1])
check_gate_defn(qc.gates[3], "CNOT", [1], [0])
check_gate_defn(qc.gates[4], "SNOT", [0])
check_gate_defn(qc.gates[5], "SNOT", [1])
check_gate_defn(qc.gates[6], "SNOT", [0], None, [0, 1], 0)
check_measurement_defn(qc.gates[7], "M", [0], 0)
check_measurement_defn(qc.gates[8], "M", [1], 1)
def test_qasm_teleportation():
filename = "teleportation.qasm"
filepath = Path(__file__).parent / 'qasm_files' / filename
teleportation = read_qasm(filepath)
final_measurement = Measurement("start", targets=[2])
initial_measurement = Measurement("start", targets=[0])
state = tensor(rand_ket(2), basis(2, 0), basis(2, 0))
_, initial_probabilities = initial_measurement.measurement_comp_basis(state)
states, probabilites = teleportation.run_statistics(state)
for i, state in enumerate(states):
final = state
prob = probabilites[i]
_, final_probabilities = final_measurement.measurement_comp_basis(final)
np.testing.assert_allclose(initial_probabilities,
final_probabilities)
assert prob == pytest.approx(0.25, abs=1e-7)
def test_qasm_teleportation():
filename = "teleportation.qasm"
filepath = Path(__file__).parent / 'qasm_files' / filename
with pytest.warns(UserWarning, match="not preserved in QubitCircuit"):
teleportation = read_qasm(filepath)
final_measurement = Measurement("start", targets=[2])
initial_measurement = Measurement("start", targets=[0])
state = tensor(rand_ket(2), basis(2, 0), basis(2, 0))
_, initial_probabilities = initial_measurement.measurement_comp_basis(state)
teleportation_results = teleportation.run_statistics(state)
states = teleportation_results.get_final_states()
probabilities = teleportation_results.get_probabilities()
for i, state in enumerate(states):
final = state
prob = probabilities[i]
_, final_probabilities = final_measurement.measurement_comp_basis(final)
np.testing.assert_allclose(initial_probabilities,
final_probabilities)
assert prob == pytest.approx(0.25, abs=1e-7)
def test_qasm_errors(filename, error, error_message):
filepath = Path(__file__).parent / 'qasm_files' / filename
with pytest.raises(error) as exc_info:
read_qasm(filepath)
assert error_message in str(exc_info.value)
