def test_teleportation():
circuit = Circuit()
circuit.add(Qubits(3))
circuit.add(CNOT(1, 2))
circuit.add(CNOT(0, 1))
circuit.add(Hadamard(0))
circuit.add(Measure([0, 1]))
result = circuit.run(1024).result
# All 16 states should be relatively equal probability
if result["11"] == 1024:
circuit.add(PauliX(2))
circuit.add(PauliZ(2))
elif result["10"] == 1024:
circuit.add(PauliZ(2))
elif result["01"] == 1024:
circuit.add(PauliX(2))
circuit.add(Measure([2]))
result2 = circuit.run(1024).result
# TODO: Fix this test.
assert result2
def test_cx_int():
circuit_1 = Circuit()
circuit_1.add(Qubits(2))
circuit_1.add(Hadamard(0))
circuit_1.add(Cx(0, 1))
circuit_1.add(Measure(0, 1))
circuit_2 = Circuit()
circuit_2.add(Qubits(2))
circuit_2.add(Hadamard(1))
circuit_2.add(Cx(0, 1))
circuit_2.add(Measure(0, 1))
sess = QSession(backend=QuantumSimulator())
result_1 = sess.run(circuit_1, num_shots=1024)
result_2 = sess.run(circuit_2, num_shots=1024)
assert result_1['01'] == 0
assert result_1['10'] > 450
assert result_1['00'] > 450
assert result_1['11'] == 0
assert result_2['01'] == 0
assert result_2['10'] == 0
assert result_2['00'] > 450
assert result_2['11'] > 450
def simple_circuit():
qc = QC()
qc.add(Qbits(3))
qc.add(H(1))
qc.add(CNOT(1, 0))
qc.add(Measure([0, 1]))
return qc
def test_apis():
qbits = Qbits(4)
# Qbit
# id
# name
# Implements iterable, so we get all of the nice
# python indexing functionality
# even = qbits[::2] # Get the even indexed qbits
# odd = qbits[1::2] # Get the odd indexed qbits
# TODO: Support multiple qbits, and this test.
# TODO: Add empty circuit test
# TODO: Infer qbits used from ciruit rather than from Qbits layer
# H(even) does not work.
qc = QuantumCircuit()
qc += Qbits(1)
# qc += H(even)
# qc += X(odd)
qc += H(0)
qc.add(Measure([qbits]))
job = qc.run(100)
result = job.result
assert result["1111"] == 0
# Any single qbit should accept
# A single qbit
H(qbits[1])
def test_u3_int():
circuit_1 = Circuit()
circuit_1.add(Qubits(1))
circuit_1.add(U3(0, theta=np.pi / 2, phi=-np.pi / 2, lam=np.pi / 2))
circuit_1.add(Measure([0]))
circuit_2 = Circuit()
circuit_2.add(Qubits(1))
circuit_2.add(Rx(theta=np.pi))
circuit_2.add(Measure([0]))
result_1 = circuit_1.run(1024).result
result_2 = circuit_2.run(1024).result
assert result_1["0"] > 450
assert result_1["1"] > 450
assert result_2["0"] > 450
assert result_2["1"] > 450
def test_u2_int():
circuit_1 = Circuit()
circuit_1.add(Qubits(1))
circuit_1.add(U2(0, phi=-np.pi / 2, alpha=np.pi / 2))
circuit_1.add(Measure([0]))
result_1 = circuit_1.run(1024).result
assert result_1["0"] > 450
assert result_1["1"] > 450
def test_t_integration():
circuit = Circuit()
circuit.add(Qubits(1))
circuit.add(T(0)) # Can also use H()
circuit.add(Measure([0]))
job = circuit.run(1024)
result = job.result
assert result["0"] == 1024
assert result["1"] == 0
def test_ry_int():
circuit = Circuit()
circuit.add(Qubits(1))
circuit.add(Ry(0, angle=np.pi / 2))
circuit.add(Measure([0]))
job = circuit.run(1000)
result = job.result
assert result["0"] > 450
assert result["1"] > 450
def test_inity_int():
circuit_1 = Circuit()
circuit_1.add(Qubits(1))
circuit_1.add(Init_y(0))
circuit_1.add(Measure([0]))
result_1 = circuit_1.run(1024).result
assert result_1["0"] > 450
assert result_1["1"] > 450
def test_u1_integration():
circuit = Circuit()
circuit.add(Qubits(1))
circuit.add(PauliX(0))
circuit.add(U1(0))
circuit.add(Measure([0]))
job = circuit.run(1024)
result = job.result
assert result["0"] == 0
assert result["1"] == 1024
def test_id_qubit():
circuit = Circuit()
circuit.add(Qubits(1))
circuit.add(ID(0))
circuit.add(Measure([0]))
job = circuit.run(1024)
result = job.result
# Accounting for random noise, results won't be exact
assert result["1"] == 0
assert result["0"] == 1024
def test_pauliy_integration():
circuit = Circuit()
circuit.add(Qubits(1))
circuit.add(PauliY(0)) # Can also use H()
circuit.add(Measure([0]))
job = circuit.run(1024)
result = job.result
# Accounting for random noise, results won't be exact
assert result["0"] == 0
assert result["1"] == 1024
def test_modulus_circuit():
circuit = Circuit()
circuit.add(Qbits(5))
circuit.add(quantum_amod_15(4))
circuit.add(Measure([0, 1]))
job = circuit.run(1024)
result = job.result
print(result)
def test_U3_int():
circuit_1 = Circuit()
circuit_1.add(Qubits(1))
circuit_1.add(U3(0, theta=np.pi / 2, phi=-np.pi / 2, alpha=np.pi / 2))
circuit_1.add(Measure(0))
circuit_2 = Circuit()
circuit_2.add(Qubits(1))
circuit_2.add(Rx(theta=np.pi))
circuit_2.add(Measure(0))
sess = QSession(backend=QuantumSimulator())
result_1 = sess.run(circuit_1, num_shots=1024)
result_2 = sess.run(circuit_2, num_shots=1024)
assert result_1['0'] > 450
assert result_1['1'] > 450
assert result_2['0'] > 450
assert result_2['1'] > 450
def test_s_integration():
circuit = Circuit()
circuit.add(Qubits(1))
circuit.add(S(0)) # Can also use H()
circuit.add(Measure())
sess = QSession(backend=QuantumSimulator())
result = sess.run(circuit, num_shots=1024)
assert result['0'] == 1024
assert result['1'] == 0
def test_crz_integration():
circuit = Circuit()
circuit.add(Qubits(2))
circuit.add(CRZ(0, 1, angle=math.pi / 3))
circuit.add(Measure(0, 1))
job = circuit.run(1024)
result = job.result
assert result["11"] == 0
assert result["00"] == 1024
assert result["10"] == 0
assert result["01"] == 0
def test_rz_int():
circuit = Circuit()
circuit.add(Qubits(1))
circuit.add(Rz(0, angle=np.pi / 2))
circuit.add(Measure(0))
sess = QSession(backend=QuantumSimulator())
result = sess.run(circuit, num_shots=1000)
assert result['0'] == 1000
assert result['1'] == 0
def test_single_qubit():
circuit = Circuit()
circuit.add(Qubits(1))
circuit.add(Hadamard(0))
circuit.add(Measure([0]))
job = circuit.run(1024, provider=QiskitProvider())
result = job.result
# Accounting for random noise, results won't be exact
assert result[bin(0)] > 450
assert result[bin(1)] > 450
def test_crk_int():
circuit = Circuit()
circuit.add(Qubits(2))
circuit.add(CRk(0, 1, k=2))
circuit.add(Measure(0, 1))
result = circuit.run(1000).result
assert result["00"] == 1000
assert result["01"] == 0
assert result["10"] == 0
assert result["11"] == 0
def test_u1_integration():
circuit = Circuit()
circuit.add(Qubits(1))
circuit.add(PauliX(0))
circuit.add(U1(0))
circuit.add(Measure())
sess = QSession(backend=QuantumSimulator())
result = sess.run(circuit, num_shots=1024)
assert result['0'] == 0
assert result['1'] == 1024
def test_ID_qubit():
circuit = Circuit()
circuit.add(Qubits(1))
circuit.add(ID(0))
circuit.add(Measure())
sess = QSession(backend=QuantumSimulator())
result = sess.run(circuit, num_shots=1024)
# Accounting for random noise, results won't be exact
assert result['1'] == 0
assert result['0'] == 1024
def test_pauliy_integration():
circuit = Circuit()
circuit.add(Qubits(1))
circuit.add(PauliY(0)) # Can also use H()
circuit.add(Measure())
sess = QSession(backend=QuantumSimulator())
result = sess.run(circuit, num_shots=1024)
# Accounting for random noise, results won't be exact
assert result['0'] == 0
assert result['1'] == 1024
def test_entanglement():
circuit = Circuit()
circuit.add(Qubits(2))
circuit.add(Hadamard(0))
circuit.add(CNOT(0, 1))
circuit.add(Measure([0, 1]))
job = circuit.run(1024)
result = job.result
assert result["01"] == 0
assert result["10"] == 0
assert result["00"] > 450
assert result["11"] > 450
def test_swap_integration():
circuit = Circuit()
circuit.add(Qubits(2))
circuit.add(PauliX(0))
circuit.add(SWAP(0, 1))
circuit.add(Measure([0, 1]))
job = circuit.run(1024)
result = job.result
assert result["11"] == 0
assert result["00"] == 0
assert result["01"] == 0
assert result["10"] == 1024
def test_mult_gate_inputs1():
circuit_1 = Circuit()
circuit_1.add(Qubits(2))
circuit_1.add(H([0, 1]))
circuit_1.add(Measure(0, 1))
result_1 = circuit_1.run(1024).result
assert result_1["00"] > 215
assert result_1["11"] > 215
assert result_1["10"] > 215
assert result_1["01"] > 215
def test_circuit_on_simulator(self):
qc = QC()
qc.add(Qbits(3))
qc.add(H(1))
qc.add(CNOT(1, 0))
qc.add(Measure([0, 1]))
ibm_provider = IBMQ()
job = qc.run(1024, ibm_provider)
result = job.result
counts = result.counts
assert counts[0] > 450 < counts[3]
def test_cr_int():
circuit = Circuit()
circuit.add(Qubits(2))
circuit.add(Cr(0, 1, angle=np.pi / 2))
circuit.add(Measure(0, 1))
result = circuit.run(1000).result
assert result["00"] == 1000
assert result["01"] == 0
assert result["10"] == 0
assert result["11"] == 0
def test_multi_hadamard():
circuit = Circuit()
circuit.add(Qubits(4))
circuit.add(Hadamard(0))
circuit.add(Hadamard(1))
circuit.add(Hadamard(2))
circuit.add(Hadamard(3))
circuit.add(Measure([0, 1, 2, 3]))
job = circuit.run(1024)
result = job.result
# All 16 states should be relatively equal probability
assert len(result.counts) == 16
assert max(result.counts.values()) - min(result.counts.values()) < 50
def test_cz2_int():
circuit_1 = Circuit()
circuit_1.add(Qubits(2))
circuit_1.add(PauliX(0))
circuit_1.add(Cz(0, 1))
circuit_1.add(Measure(0, 1))
result_1 = circuit_1.run(1024).result
assert result_1["00"] == 0
assert result_1["11"] == 0
assert result_1["01"] == 1024
assert result_1["10"] == 0
def test_unitary_symmetry_does_nothing():
symmetric_circuit_1 = Circuit()
symmetric_circuit_1.add(Qubits(2))
symmetric_circuit_1.add(Hadamard(0))
symmetric_circuit_1.add(Hadamard(0))
symmetric_circuit_1.add(CNOT(0, 1))
symmetric_circuit_1.add(CNOT(0, 1))
symmetric_circuit_1.add(Measure([0, 1]))
symmetric_circuit_2 = Circuit()
symmetric_circuit_2.add(Qubits(2, state=1))
symmetric_circuit_2.add(Hadamard(0))
symmetric_circuit_2.add(CNOT(0, 1))
symmetric_circuit_2.add(CNOT(0, 1))
symmetric_circuit_2.add(Hadamard(0))
symmetric_circuit_2.add(Measure([0, 1]))
result_1 = symmetric_circuit_1.run(1024).result
result_2 = symmetric_circuit_2.run(1024).result
assert result_1.counts.get(0) == 1024
assert result_2.counts.get(0) == 1024