def test_decomp_stdgates(): gate0 = qf.I(0, 1) gate1 = qf.canonical_decomposition(gate0).asgate() assert qf.gates_close(gate0, gate1) gate0 = qf.CNOT(0, 1) gate1 = qf.canonical_decomposition(gate0).asgate() assert qf.gates_close(gate0, gate1) gate0 = qf.SWAP(0, 1) gate1 = qf.canonical_decomposition(gate0).asgate() assert qf.gates_close(gate0, gate1) gate0 = qf.ISWAP(0, 1) gate1 = qf.canonical_decomposition(gate0).asgate() assert qf.gates_close(gate0, gate1) gate0 = qf.CNOT(0, 1)**0.5 gate1 = qf.canonical_decomposition(gate0).asgate() assert qf.gates_close(gate0, gate1) gate0 = qf.SWAP(0, 1)**0.5 gate1 = qf.canonical_decomposition(gate0).asgate() assert qf.gates_close(gate0, gate1) gate0 = qf.ISWAP(0, 1)**0.5 gate1 = qf.canonical_decomposition(gate0).asgate() assert qf.gates_close(gate0, gate1)
def test_gate_inverse(): inv = qf.S().H eye = qf.S() @ inv assert qf.gates_close(eye, qf.I()) inv = qf.ISWAP().H eye = qf.ISWAP() @ inv assert qf.gates_close(eye, qf.identity_gate(2))
def test_pswap(): for _ in range(REPS): theta = random.uniform(-4 * pi, +4 * pi) assert qf.almost_unitary(qf.PSWAP(theta)) assert qf.gates_close(qf.SWAP(), qf.PSWAP(0)) assert qf.gates_close(qf.ISWAP(), qf.PSWAP(pi / 2))
def test_gates_to_latex(): circ = qf.Circuit() circ += qf.I(7) circ += qf.X(0) circ += qf.Y(1) circ += qf.Z(2) circ += qf.H(3) circ += qf.S(4) circ += qf.T(5) circ += qf.S_H(6) circ += qf.T_H(7) circ += qf.RX(-0.5*pi, 0) circ += qf.RY(0.5*pi, 1) circ += qf.RZ((1/3)*pi, 1) circ += qf.RY(0.222, 1) circ += qf.TX(0.5, 0) circ += qf.TY(0.5, 1) circ += qf.TZ(0.4, 1) circ += qf.TZ(0.47276, 1) # Gate with cunning hack gate = qf.RZ(0.4, 1) gate.params['theta'] = qf.Parameter('\\theta') circ += gate circ += qf.CNOT(1, 2) circ += qf.CNOT(2, 1) circ += qf.CZ(1, 3) circ += qf.SWAP(1, 5) circ += qf.ISWAP(4, 2) # circ += qf.Barrier(0, 1, 2, 3, 4, 5, 6) # Not yet supported circ += qf.CCNOT(1, 2, 3) circ += qf.CSWAP(4, 5, 6) circ += qf.P0(0) circ += qf.P1(1) circ += qf.Reset(2) circ += qf.Reset(4, 5, 6) circ += qf.H(4) # circ += qf.Reset() # FIXME. Should fail with clear error message circ += qf.XX(0.25, 1, 3) circ += qf.YY(0.75, 1, 3) circ += qf.ZZ(1/3, 3, 1) circ += qf.Measure(0) latex = qf.circuit_to_latex(circ) print(latex)
def test_piswap(): for _ in range(REPS): theta = random.uniform(-4 * pi, +4 * pi) assert qf.almost_unitary(qf.PISWAP(theta)) for _ in range(REPS): theta = random.uniform(0, +pi) assert qf.gates_close(qf.PISWAP(0), qf.identity_gate(2)) assert qf.gates_close(qf.PISWAP(pi / 4), qf.ISWAP())
def prepare_w4(): """ Prepare a 4-qubit W state using sqrt(iswaps) and local gates """ circ = qf.Circuit() circ += qf.X(1) circ += qf.ISWAP(1, 2)**0.5 circ += qf.S(2) circ += qf.Z(2) circ += qf.ISWAP(2, 3)**0.5 circ += qf.S(3) circ += qf.Z(3) circ += qf.ISWAP(0, 1)**0.5 circ += qf.S(0) circ += qf.Z(0) ket = circ.run() return ket
def test_gatepow(): gates = [ qf.I(), qf.X(), qf.Y(), qf.Z(), qf.H(), qf.S(), qf.T(), qf.PHASE(0.1), qf.RX(0.2), qf.RY(0.3), qf.RZ(0.4), qf.CZ(), qf.CNOT(), qf.SWAP(), qf.ISWAP(), qf.CPHASE00(0.5), qf.CPHASE01(0.6), qf.CPHASE10(0.6), qf.CPHASE(0.7), qf.PSWAP(0.15), qf.CCNOT(), qf.CSWAP(), qf.TX(2.7), qf.TY(1.2), qf.TZ(0.3), qf.ZYZ(3.5, 0.9, 2.1), qf.CANONICAL(0.1, 0.2, 7.4), qf.XX(1.8), qf.YY(0.9), qf.ZZ(0.45), qf.PISWAP(0.2), qf.EXCH(0.1), qf.TH(0.3) ] for gate in gates: assert qf.gates_close(gate.H, gate**-1) for gate in gates: sqrt_gate = gate**(1 / 2) two_gate = sqrt_gate @ sqrt_gate assert qf.gates_close(gate, two_gate) for gate in gates: gate0 = gate**0.3 gate1 = gate**0.7 gate2 = gate0 @ gate1 assert qf.gates_close(gate, gate2) for K in range(1, 5): gate = qf.random_gate(K) # FIXME: Throw error on K=0 sqrt_gate = gate**0.5 two_gate = sqrt_gate @ sqrt_gate assert qf.gates_close(gate, two_gate) for gate in gates: rgate = qf.Gate((gate**0.5).tensor) tgate = rgate @ rgate assert qf.gates_close(gate, tgate)
def test_unitary_2qubit(): assert qf.almost_unitary(qf.CZ()) assert qf.almost_unitary(qf.CNOT()) assert qf.almost_unitary(qf.SWAP()) assert qf.almost_unitary(qf.ISWAP())
def test_almost_hermitian(): assert qf.almost_hermitian(qf.X()) assert not qf.almost_hermitian(qf.ISWAP())
def w16_circuit() -> qf.Circuit: """ Return a circuit that prepares the the 16-qubit W state using\ sqrt(iswaps) and local gates, respecting linear topology """ gates = [ qf.X(7), qf.ISWAP(7, 8)**0.5, qf.S(8), qf.Z(8), qf.SWAP(7, 6), qf.SWAP(6, 5), qf.SWAP(5, 4), qf.SWAP(8, 9), qf.SWAP(9, 10), qf.SWAP(10, 11), qf.ISWAP(4, 3)**0.5, qf.S(3), qf.Z(3), qf.ISWAP(11, 12)**0.5, qf.S(12), qf.Z(12), qf.SWAP(3, 2), qf.SWAP(4, 5), qf.SWAP(11, 10), qf.SWAP(12, 13), qf.ISWAP(2, 1)**0.5, qf.S(1), qf.Z(1), qf.ISWAP(5, 6)**0.5, qf.S(6), qf.Z(6), qf.ISWAP(10, 9)**0.5, qf.S(9), qf.Z(9), qf.ISWAP(13, 14)**0.5, qf.S(14), qf.Z(14), qf.ISWAP(1, 0)**0.5, qf.S(0), qf.Z(0), qf.ISWAP(2, 3)**0.5, qf.S(3), qf.Z(3), qf.ISWAP(5, 4)**0.5, qf.S(4), qf.Z(4), qf.ISWAP(6, 7)**0.5, qf.S(7), qf.Z(7), qf.ISWAP(9, 8)**0.5, qf.S(8), qf.Z(8), qf.ISWAP(10, 11)**0.5, qf.S(11), qf.Z(11), qf.ISWAP(13, 12)**0.5, qf.S(12), qf.Z(12), qf.ISWAP(14, 15)**0.5, qf.S(15), qf.Z(15), ] circ = qf.Circuit(gates) return circ