def run(self, circuit: Circuit, **kwargs): num_shots = kwargs.get('num_shots', 1024) initial_state = circuit.initial_state() gates = circuit.to_gates() measure_bits = circuit.measure_bits() counts = self.backend.run(initial_state, gates, measure_bits, num_shots) return QResult(counts, sig_bits=len(circuit.initial_state()))
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_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_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_entanglement(): circuit = Circuit() circuit.add(Qubits(2)) circuit.add(Hadamard(0)) circuit.add(CNOT(0, 1)) circuit.add(Measure(0, 1)) sess = QSession(backend=QuantumSimulator()) result = sess.run(circuit, num_shots=1024) assert result['01'] == 0 assert result['10'] == 0 assert result['00'] > 450 assert result['11'] > 450
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_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_Cz_int(): circuit = Circuit() circuit.add(Qubits(2)) circuit.add(Hadamard(0)) circuit.add(Hadamard(1)) circuit.add(Cz(0, 1)) circuit.add(Measure(0, 1)) sess = QSession(backend=QuantumSimulator()) result = sess.run(circuit, num_shots=1000) assert result['00'] > 210 assert result['01'] > 210 assert result['10'] > 210 assert result['11'] > 210
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)) sess = QSession(backend=QuantumSimulator()) result = sess.run(circuit, num_shots=1024) # 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_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_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_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_circuit_add_wrong_type(): class SomeClass(object): def __init__(self): pass with pytest.raises(TypeError) as type_error: circuit = Circuit().add(SomeClass())
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_circuit_add_circuit(): circuit = Circuit().add(_BaseLayer()).add(_BaseLayer()) circuit_to_add = Circuit() circuit_to_add.add(_BaseLayer()).add(_BaseLayer()).add(_BaseLayer()) assert len(circuit.layers) == 2 circuit.add(circuit_to_add) assert len(circuit.layers) == 5
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_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)) sess = QSession(backend=QuantumSimulator()) result_1 = sess.run(circuit_1, num_shots=1024) assert result_1['0'] > 450 assert result_1['1'] > 450
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 find_period(a, N): """ WIP: Quantum subroutine for shor's algorithm. Finds the period of a function of the form: f(x) = a^x % N This uses the quantum fourier transform. """ circuit = Circuit() circuit.add(Qubits(5)) circuit.add(QFT(0, 1, 2, 3)) circuit.add(quantum_amod_15(a)) circuit.add(QFT(3, 2, 1, 0)) # Inverse Quantum Fourier transform from shor.backends import QuantumSimulator, QSession sess = QSession(backend=QuantumSimulator()) result = sess.run(circuit, num_shots=1024) return result
def test_ccnot_integration(): circuit = Circuit() circuit.add(Qubits(3)) circuit.add(PauliX(0)) circuit.add(PauliX(1)) circuit.add(CCNOT(0, 1, 2)) circuit.add(Measure(0, 1, 2)) sess = QSession(backend=QuantumSimulator()) result = sess.run(circuit, num_shots=1024) assert result['000'] == 0 assert result['001'] == 0 assert result['010'] == 0 assert result['100'] == 0 assert result['110'] == 0 assert result['101'] == 0 assert result['011'] == 0 assert result['111'] == 1024
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_pauliz_integration(): circuit = Circuit() circuit.add(Qubits(1)) circuit.add(Hadamard(0)) circuit.add(PauliZ(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'] > 450 assert result['1'] > 450
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_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_swap_integration(): # circuit = Circuit() circuit.add(Qubits(2)) circuit.add(PauliX(0)) circuit.add(SWAP(0, 1)) circuit.add(Measure(0, 1)) sess = QSession(backend=QuantumSimulator()) result = sess.run(circuit, num_shots=1024) assert result['11'] == 0 assert result['00'] == 0 assert result['10'] == 0 assert result['01'] == 1024
def qft(qubits: List[int]) -> Circuit: qc = Circuit() for i in range(len(qubits)): for k in range(i): qc.add(Rx(qubits[i], qubits[k], angle=math.pi/float(2**(i-k)))) qc.add(H(qubits[i])) return qc