def test_acceptance_with_qiskit():
"""Test gate operations with QASM interface"""
circuit = Circuit()
circuit += ops.RotateX(0, -np.pi)
circuit += ops.RotateY(0, -np.pi)
circuit += ops.RotateZ(0, -np.pi)
circuit += ops.CNOT(0, 1)
circuit += ops.Hadamard(0)
circuit += ops.PauliX(0)
circuit += ops.PauliY(0)
circuit += ops.PauliZ(0)
circuit += ops.SGate(0)
circuit += ops.TGate(0)
circuit += ops.SqrtPauliX(0)
circuit += ops.MolmerSorensenXX(0, 1)
circuit += ops.ControlledPauliY(0, 1)
circuit += ops.ControlledPauliZ(0, 1)
circuit += ops.SingleQubitGate(0, 1, 0, 1, 0, 1.0)
circuit += ops.PragmaRepeatedMeasurement('ro', 1, None)
circuit += ops.MeasureQubit(0, 'ro', 0)
circuit += ops.DefinitionFloat(name='rof', length=1, is_output=True)
circuit += ops.DefinitionBit(name='ro', length=2, is_output=True)
circuit += ops.DefinitionComplex(name='roc', length=1, is_output=True)
circuit += ops.InputSymbolic('other', 0)
circuit += ops.PragmaSetNumberOfMeasurements(20, 'ro')
backend = QasmBackend(number_qubits=2)
backend.run_circuit(circuit=circuit, overwrite=True)
q_circuit = QuantumCircuit.from_qasm_file(
"default_qasm_backend_output.qasm")
# Unless required by applicable law or agreed to in writing, software distributed under the License
# is distributed on an "AS IS" BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express
# or implied. See the License for the specific language governing permissions and limitations under
# the License.
import pytest
import sys
import numpy as np
import numpy.testing as npt
from qoqo import operations as ops
from qoqo import Circuit
from typing import List
@pytest.mark.parametrize('measurement', [
(ops.MeasureQubit(qubit=0, readout='ro', readout_index=0), int, 0),
(ops.PragmaRepeatedMeasurement(readout='ro',
number_measurements=10), int, 0),
(ops.PragmaGetPauliProduct(
qubit_paulis={
0: 1,
1: 2
}, readout='ro', circuit=Circuit()), float, 1),
(ops.PragmaGetOccupationProbability(readout='ro',
circuit=Circuit()), float, 1),
(ops.PragmaGetStateVector(readout='ro', circuit=Circuit()), complex, 2),
(ops.PragmaGetDensityMatrix(readout='ro', circuit=Circuit()), complex, 2),
])
def test_circuit(measurement):
"""Test building a circuit"""
circuit = Circuit()
circuit += ops.DefinitionFloat(name='ro', length=1, is_output=True)
circuit += ops.DefinitionComplex(name='ro', length=1, is_output=True)
ops.PragmaSetDensityMatrix(np.array([[0, 0], [0, 1]], dtype=complex)),
ops.PragmaDamping(0, 0.005, 0.02),
ops.PragmaDephasing(0, 0.005, 0.02),
ops.PragmaDepolarising(0, 0.005, 0.02),
ops.PragmaRandomNoise(0, 0.005, 0.02, 0.01),
ops.PragmaGeneralNoise(
0, 0.005, np.array([[1, 0, 0], [0, 1, 0], [0, 0, 1]], dtype=float)),
ops.PragmaConditional('ro', 0, Circuit()),
ops.PragmaRepeatGate(3),
ops.PragmaBoostNoise(0.004),
ops.PragmaStopParallelBlock([0, 1], 0.002),
ops.PragmaSleep([0, 1], 0.002),
ops.PragmaActiveReset(0),
ops.PragmaGlobalPhase(0.03),
ops.MeasureQubit(0, 'ro', 0),
ops.PragmaRepeatedMeasurement('ro', 20, {}),
ops.PragmaGetStateVector('ro', Circuit()),
ops.PragmaGetDensityMatrix('ro', Circuit()),
ops.PragmaGetOccupationProbability('ro', Circuit()),
ops.PragmaGetPauliProduct({}, 'ro', Circuit()),
ops.PragmaStartDecompositionBlock([0, 1], {}),
ops.PragmaStopDecompositionBlock([0, 1]),
ops.InputSymbolic('other', 0),
])
def test_all_operations(op: Any):
"""Test all operations with mocked interface"""
circuit = Circuit()
circuit += ops.DefinitionFloat(name='ro', length=1, is_output=False)
circuit += ops.DefinitionBit(name='ro', length=1, is_output=False)
circuit += ops.DefinitionComplex(name='ro', length=1, is_output=False)
circuit += ops.DefinitionUsize(name='ro', length=1, is_output=False)
(ops.RotateX(0, -np.pi), 'rx(-3.141592653589793) q[0]'),
(ops.RotateY(0, -np.pi), 'ry(-3.141592653589793) q[0]'),
(ops.RotateZ(0, -np.pi), 'rz(-3.141592653589793) q[0]'),
(ops.CNOT(0, 1), 'cx q[0],q[1]'),
(ops.Hadamard(0), 'h q[0]'),
(ops.PauliX(0), 'x q[0]'),
(ops.PauliY(0), 'y q[0]'),
(ops.PauliZ(0), 'z q[0]'),
(ops.SGate(0), 's q[0]'),
(ops.TGate(0), 't q[0]'),
(ops.SqrtPauliX(0), 'rx(1.5707963267948966) q[0]'),
(ops.MolmerSorensenXX(0, 1), 'rxx(pi/2) q[0],q[1]'),
(ops.ControlledPauliY(0, 1), 'cy q[0],q[1]'),
(ops.ControlledPauliZ(0, 1), 'cz q[0],q[1]'),
(ops.SingleQubitGate(0, 1, 0, 1, 0, 1.0), 'u3(0.0,0.0,-0.0) q[0]'),
(ops.PragmaRepeatedMeasurement('ro', 1, None), 'measure q -> ro;\n'),
(ops.PragmaRepeatedMeasurement('ro', 1, {0: 1, 1: 0}),
'measure q[1] -> ro[0];\nmeasure q[0] -> ro[1];\n'),
(ops.MeasureQubit(0, 'ro', 0), 'measure q[0] -> ro[0]'),
(ops.DefinitionFloat(name='ro', length=1, is_output=True), 'creg ro[1]'),
(ops.DefinitionUsize(name='ro', length=1, is_output=True), 'creg ro[1]'),
(ops.DefinitionBit(name='ro', length=1, is_output=True), 'creg ro[1]'),
(ops.DefinitionComplex(name='ro', length=1, is_output=True), 'creg ro[1]'),
(ops.InputSymbolic('other', 0), None),
(ops.PragmaSetNumberOfMeasurements(20, 'ro'), None)
])
def test_gate_translation(gate: Tuple[Any, str]):
"""Test gate operations with QASM interface"""
qasm_operation = qasm_call_operation(operation=gate[0],
number_qubits=2)