def test_three_qubit_matrix_to_operations_errors():
a, b, c = cirq.LineQubit.range(3)
with pytest.raises(ValueError, match="(8,8)"):
cirq.three_qubit_matrix_to_operations(a, b, c, np.eye(2))
with pytest.raises(ValueError, match="not unitary"):
cirq.three_qubit_matrix_to_operations(a, b, c,
cirq.unitary(cirq.CCX) * 2)
def test_three_qubit_matrix_to_operations(u):
a, b, c = cirq.LineQubit.range(3)
operations = cirq.three_qubit_matrix_to_operations(a, b, c, u)
final_circuit = cirq.Circuit(operations)
final_unitary = final_circuit.unitary(
qubits_that_should_be_present=[a, b, c])
cirq.testing.assert_allclose_up_to_global_phase(u,
final_unitary,
atol=1e-9)
num_two_qubit_gates = len([
op for op in list(final_circuit.all_operations())
if isinstance(op.gate, (cirq.CZPowGate, cirq.CNotPowGate))
])
assert num_two_qubit_gates <= 20, f"expected at most 20 CZ/CNOTs got {num_two_qubit_gates}"
def matrix_to_sycamore_operations(
target_qubits: List[cirq.GridQubit],
matrix: np.ndarray) -> Tuple[cirq.OP_TREE, List[cirq.GridQubit]]:
"""A method to convert a unitary matrix to a list of Sycamore operations.
This method will return a list of `cirq.Operation`s using the qubits and (optionally) ancilla
qubits to implement the unitary matrix `matrix` on the target qubits `qubits`.
The operations are also supported by `cirq.google.gate_sets.SYC_GATESET`.
Args:
target_qubits: list of qubits the returned operations will act on. The qubit order defined by the list
is assumed to be used by the operations to implement `matrix`.
matrix: a matrix that is guaranteed to be unitary and of size (2**len(qs), 2**len(qs)).
Returns:
A tuple of operations and ancilla qubits allocated.
Operations: In case the matrix is supported, a list of operations `ops` is returned.
`ops` acts on `qs` qubits and for which `cirq.unitary(ops)` is equal to `matrix` up
to certain tolerance. In case the matrix is not supported, it might return NotImplemented to
reduce the noise in the judge output.
Ancilla qubits: In case ancilla qubits are allocated a list of ancilla qubits. Otherwise
an empty list.
.
"""
num_qubit = len(target_qubits)
ops = []
if num_qubit == 1:
ops = cirq.optimizers.single_qubit_matrix_to_gates(matrix)
for i in range(len(ops)):
ops[i] = ops[i].on(target_qubits[0])
elif num_qubit == 2:
if np.all(matrix == cirq.unitary(cirq.google.SYC)):
ops = [cirq.google.SYC(target_qubits[0], target_qubits[1])]
else:
ops = cirq.two_qubit_matrix_to_operations(target_qubits[0],
target_qubits[1],
matrix,
allow_partial_czs=False)
elif num_qubit == 3:
if (matrix == np.diag([1, 1, 1, 1, 1, 1, 1, 1])).all():
ops = []
else:
ops = cirq.three_qubit_matrix_to_operations(
target_qubits[0], target_qubits[1], target_qubits[2], matrix)
# print('target_qubit:', target_qubits)
# print('matrix:', matrix, end='\n\n')
# print('ops:')
# for i in ops:
# print(i)
return ops, []
def matrix_to_sycamore_operations(
target_qubits: List[cirq.GridQubit],
matrix: np.ndarray) -> Tuple[cirq.OP_TREE, List[cirq.GridQubit]]:
"""A method to convert a unitary matrix to a list of Sycamore operations.
This method will return a list of `cirq.Operation`s using the qubits and (optionally) ancilla
qubits to implement the unitary matrix `matrix` on the target qubits `qubits`.
The operations are also supported by `cirq.google.gate_sets.SYC_GATESET`.
Args:
target_qubits: list of qubits the returned operations will act on. The qubit order defined by the list
is assumed to be used by the operations to implement `matrix`.
matrix: a matrix that is guaranteed to be unitary and of size (2**len(qs), 2**len(qs)).
Returns:
A tuple of operations and ancilla qubits allocated.
Operations: In case the matrix is supported, a list of operations `ops` is returned.
`ops` acts on `qs` qubits and for which `cirq.unitary(ops)` is equal to `matrix` up
to certain tolerance. In case the matrix is not supported, it might return NotImplemented to
reduce the noise in the judge output.
Ancilla qubits: In case ancilla qubits are allocated a list of ancilla qubits. Otherwise
an empty list.
.
"""
ops = NotImplemented
if len(target_qubits) == 1:
ops = cirq.single_qubit_matrix_to_gates(matrix)
for i, x in enumerate(ops):
ops[i] = cirq.GateOperation(x, [target_qubits[0]])
elif len(target_qubits) == 2:
ops = cirq.two_qubit_matrix_to_operations(target_qubits[0],
target_qubits[1],
matrix,
allow_partial_czs=False,
atol=1e-4)
elif len(target_qubits) == 3:
ops = cirq.three_qubit_matrix_to_operations(target_qubits[0],
target_qubits[1],
target_qubits[2],
matrix,
atol=1e-4)
return ops, []
def test_three_qubit_matrix_to_operations_scipy_error():
a, b, c = cirq.LineQubit.range(3)
with pytest.raises(ImportError, match="three_qubit.*1.5.0+"):
cirq.three_qubit_matrix_to_operations(a, b, c, np.eye(8))
