qiskit 1.0.0b1__cp38-abi3-manylinux_2_17_aarch64.manylinux2014_aarch64.whl

qiskit/pulse/instruction_schedule_map.py

@@ -0,0 +1,406 @@
+# This code is part of Qiskit.
+#
+# (C) Copyright IBM 2019.
+#
+# This code is licensed under the Apache License, Version 2.0. You may
+# obtain a copy of this license in the LICENSE.txt file in the root directory
+# of this source tree or at http://www.apache.org/licenses/LICENSE-2.0.
+#
+# Any modifications or derivative works of this code must retain this
+# copyright notice, and modified files need to carry a notice indicating
+# that they have been altered from the originals.
+
+# pylint: disable=unused-import
+
+"""
+A convenient way to track reusable subschedules by name and qubit.
+
+This can be used for scheduling circuits with custom definitions, for instance::
+
+    inst_map = InstructionScheduleMap()
+    inst_map.add('new_inst', 0, qubit_0_new_inst_schedule)
+
+    sched = schedule(quantum_circuit, backend, inst_map)
+
+An instance of this class is instantiated by Pulse-enabled backends and populated with defaults
+(if available)::
+
+    inst_map = backend.defaults().instruction_schedule_map
+
+"""
+import functools
+import warnings
+from collections import defaultdict
+from typing import Callable, Iterable, List, Tuple, Union, Optional
+
+from qiskit.circuit.instruction import Instruction
+from qiskit.circuit.parameterexpression import ParameterExpression
+from qiskit.pulse.calibration_entries import (
+    CalibrationEntry,
+    ScheduleDef,
+    CallableDef,
+    # for backward compatibility
+    PulseQobjDef,
+    CalibrationPublisher,
+)
+from qiskit.pulse.exceptions import PulseError
+from qiskit.pulse.schedule import Schedule, ScheduleBlock
+
+
+class InstructionScheduleMap:
+    """Mapping from :py:class:`~qiskit.circuit.QuantumCircuit`
+    :py:class:`qiskit.circuit.Instruction` names and qubits to
+    :py:class:`~qiskit.pulse.Schedule` s. In particular, the mapping is formatted as type::
+
+         Dict[str, Dict[Tuple[int], Schedule]]
+
+    where the first key is the name of a circuit instruction (e.g. ``'u1'``, ``'measure'``), the
+    second key is a tuple of qubit indices, and the final value is a Schedule implementing the
+    requested instruction.
+
+    These can usually be seen as gate calibrations.
+    """
+
+    def __init__(self):
+        """Initialize a circuit instruction to schedule mapper instance."""
+        # The processed and reformatted circuit instruction definitions
+
+        # Do not use lambda function for nested defaultdict, i.e. lambda: defaultdict(CalibrationEntry).
+        # This crashes qiskit parallel. Note that parallel framework passes args as
+        # pickled object, however lambda function cannot be pickled.
+        self._map = defaultdict(functools.partial(defaultdict, CalibrationEntry))
+
+        # A backwards mapping from qubit to supported instructions
+        self._qubit_instructions = defaultdict(set)
+
+    def has_custom_gate(self) -> bool:
+        """Return ``True`` if the map has user provided instruction."""
+        for qubit_inst in self._map.values():
+            for entry in qubit_inst.values():
+                if entry.user_provided:
+                    return True
+        return False
+
+    @property
+    def instructions(self) -> List[str]:
+        """Return all instructions which have definitions.
+
+        By default, these are typically the basis gates along with other instructions such as
+        measure and reset.
+
+        Returns:
+            The names of all the circuit instructions which have Schedule definitions in this.
+        """
+        return list(self._map.keys())
+
+    def qubits_with_instruction(
+        self, instruction: Union[str, Instruction]
+    ) -> List[Union[int, Tuple[int]]]:
+        """Return a list of the qubits for which the given instruction is defined. Single qubit
+        instructions return a flat list, and multiqubit instructions return a list of ordered
+        tuples.
+
+        Args:
+            instruction: The name of the circuit instruction.
+
+        Returns:
+            Qubit indices which have the given instruction defined. This is a list of tuples if the
+            instruction has an arity greater than 1, or a flat list of ints otherwise.
+
+        Raises:
+            PulseError: If the instruction is not found.
+        """
+        instruction = _get_instruction_string(instruction)
+        if instruction not in self._map:
+            return []
+        return [
+            qubits[0] if len(qubits) == 1 else qubits
+            for qubits in sorted(self._map[instruction].keys())
+        ]
+
+    def qubit_instructions(self, qubits: Union[int, Iterable[int]]) -> List[str]:
+        """Return a list of the instruction names that are defined by the backend for the given
+        qubit or qubits.
+
+        Args:
+            qubits: A qubit index, or a list or tuple of indices.
+
+        Returns:
+            All the instructions which are defined on the qubits.
+
+            For 1 qubit, all the 1Q instructions defined. For multiple qubits, all the instructions
+            which apply to that whole set of qubits (e.g. ``qubits=[0, 1]`` may return ``['cx']``).
+        """
+        if _to_tuple(qubits) in self._qubit_instructions:
+            return list(self._qubit_instructions[_to_tuple(qubits)])
+        return []
+
+    def has(self, instruction: Union[str, Instruction], qubits: Union[int, Iterable[int]]) -> bool:
+        """Is the instruction defined for the given qubits?
+
+        Args:
+            instruction: The instruction for which to look.
+            qubits: The specific qubits for the instruction.
+
+        Returns:
+            True iff the instruction is defined.
+        """
+        instruction = _get_instruction_string(instruction)
+        return instruction in self._map and _to_tuple(qubits) in self._map[instruction]
+
+    def assert_has(
+        self, instruction: Union[str, Instruction], qubits: Union[int, Iterable[int]]
+    ) -> None:
+        """Error if the given instruction is not defined.
+
+        Args:
+            instruction: The instruction for which to look.
+            qubits: The specific qubits for the instruction.
+
+        Raises:
+            PulseError: If the instruction is not defined on the qubits.
+        """
+        instruction = _get_instruction_string(instruction)
+        if not self.has(instruction, _to_tuple(qubits)):
+            if instruction in self._map:
+                raise PulseError(
+                    "Operation '{inst}' exists, but is only defined for qubits "
+                    "{qubits}.".format(
+                        inst=instruction, qubits=self.qubits_with_instruction(instruction)
+                    )
+                )
+            raise PulseError(f"Operation '{instruction}' is not defined for this system.")
+
+    def get(
+        self,
+        instruction: Union[str, Instruction],
+        qubits: Union[int, Iterable[int]],
+        *params: Union[complex, ParameterExpression],
+        **kwparams: Union[complex, ParameterExpression],
+    ) -> Union[Schedule, ScheduleBlock]:
+        """Return the defined :py:class:`~qiskit.pulse.Schedule` or
+        :py:class:`~qiskit.pulse.ScheduleBlock` for the given instruction on the given qubits.
+
+        If all keys are not specified this method returns schedule with unbound parameters.
+
+        Args:
+            instruction: Name of the instruction or the instruction itself.
+            qubits: The qubits for the instruction.
+            *params: Command parameters for generating the output schedule.
+            **kwparams: Keyworded command parameters for generating the schedule.
+
+        Returns:
+            The Schedule defined for the input.
+        """
+        return self._get_calibration_entry(instruction, qubits).get_schedule(*params, **kwparams)
+
+    def _get_calibration_entry(
+        self,
+        instruction: Union[str, Instruction],
+        qubits: Union[int, Iterable[int]],
+    ) -> CalibrationEntry:
+        """Return the :class:`.CalibrationEntry` without generating schedule.
+
+        When calibration entry is un-parsed Pulse Qobj, this returns calibration
+        without parsing it. :meth:`CalibrationEntry.get_schedule` method
+        must be manually called with assigned parameters to get corresponding pulse schedule.
+
+        This method is expected be directly used internally by the V2 backend converter
+        for faster loading of the backend calibrations.
+
+        Args:
+            instruction: Name of the instruction or the instruction itself.
+            qubits: The qubits for the instruction.
+
+        Returns:
+            The calibration entry.
+        """
+        instruction = _get_instruction_string(instruction)
+        self.assert_has(instruction, qubits)
+
+        return self._map[instruction][_to_tuple(qubits)]
+
+    def add(
+        self,
+        instruction: Union[str, Instruction],
+        qubits: Union[int, Iterable[int]],
+        schedule: Union[Schedule, ScheduleBlock, Callable[..., Union[Schedule, ScheduleBlock]]],
+        arguments: Optional[List[str]] = None,
+    ) -> None:
+        """Add a new known instruction for the given qubits and its mapping to a pulse schedule.
+
+        Args:
+            instruction: The name of the instruction to add.
+            qubits: The qubits which the instruction applies to.
+            schedule: The Schedule that implements the given instruction.
+            arguments: List of parameter names to create a parameter-bound schedule from the
+                associated gate instruction. If :py:meth:`get` is called with arguments rather
+                than keyword arguments, this parameter list is used to map the input arguments to
+                parameter objects stored in the target schedule.
+
+        Raises:
+            PulseError: If the qubits are provided as an empty iterable.
+        """
+        instruction = _get_instruction_string(instruction)
+
+        # validation of target qubit
+        qubits = _to_tuple(qubits)
+        if qubits == ():
+            raise PulseError(f"Cannot add definition {instruction} with no target qubits.")
+
+        # generate signature
+        if isinstance(schedule, (Schedule, ScheduleBlock)):
+            entry = ScheduleDef(arguments)
+        elif callable(schedule):
+            if arguments:
+                warnings.warn(
+                    "Arguments are overruled by the callback function signature. "
+                    "Input `arguments` are ignored.",
+                    UserWarning,
+                )
+            entry = CallableDef()
+        else:
+            raise PulseError(
+                "Supplied schedule must be one of the Schedule, ScheduleBlock or a "
+                "callable that outputs a schedule."
+            )
+        entry.define(schedule, user_provided=True)
+        self._add(instruction, qubits, entry)
+
+    def _add(
+        self,
+        instruction_name: str,
+        qubits: Tuple[int, ...],
+        entry: CalibrationEntry,
+    ):
+        """A method to resister calibration entry.
+
+        .. note::
+
+            This is internal fast-path function, and caller must ensure
+            the entry is properly formatted. This function may be used by other programs
+            that load backend calibrations to create Qiskit representation of it.
+
+        Args:
+            instruction_name: Name of instruction.
+            qubits: List of qubits that this calibration is applied.
+            entry: Calibration entry to register.
+
+        :meta public:
+        """
+        self._map[instruction_name][qubits] = entry
+        self._qubit_instructions[qubits].add(instruction_name)
+
+    def remove(
+        self, instruction: Union[str, Instruction], qubits: Union[int, Iterable[int]]
+    ) -> None:
+        """Remove the given instruction from the listing of instructions defined in self.
+
+        Args:
+            instruction: The name of the instruction to add.
+            qubits: The qubits which the instruction applies to.
+        """
+        instruction = _get_instruction_string(instruction)
+        qubits = _to_tuple(qubits)
+        self.assert_has(instruction, qubits)
+
+        del self._map[instruction][qubits]
+        if not self._map[instruction]:
+            del self._map[instruction]
+
+        self._qubit_instructions[qubits].remove(instruction)
+        if not self._qubit_instructions[qubits]:
+            del self._qubit_instructions[qubits]
+
+    def pop(
+        self,
+        instruction: Union[str, Instruction],
+        qubits: Union[int, Iterable[int]],
+        *params: Union[complex, ParameterExpression],
+        **kwparams: Union[complex, ParameterExpression],
+    ) -> Union[Schedule, ScheduleBlock]:
+        """Remove and return the defined schedule for the given instruction on the given
+        qubits.
+
+        Args:
+            instruction: Name of the instruction.
+            qubits: The qubits for the instruction.
+            *params: Command parameters for generating the output schedule.
+            **kwparams: Keyworded command parameters for generating the schedule.
+
+        Returns:
+            The Schedule defined for the input.
+        """
+        instruction = _get_instruction_string(instruction)
+        schedule = self.get(instruction, qubits, *params, **kwparams)
+        self.remove(instruction, qubits)
+        return schedule
+
+    def get_parameters(
+        self, instruction: Union[str, Instruction], qubits: Union[int, Iterable[int]]
+    ) -> Tuple[str]:
+        """Return the list of parameters taken by the given instruction on the given qubits.
+
+        Args:
+            instruction: Name of the instruction.
+            qubits: The qubits for the instruction.
+
+        Returns:
+            The names of the parameters required by the instruction.
+        """
+        instruction = _get_instruction_string(instruction)
+
+        self.assert_has(instruction, qubits)
+        signature = self._map[instruction][_to_tuple(qubits)].get_signature()
+        return tuple(signature.parameters.keys())
+
+    def __str__(self):
+        single_q_insts = "1Q instructions:\n"
+        multi_q_insts = "Multi qubit instructions:\n"
+        for qubits, insts in self._qubit_instructions.items():
+            if len(qubits) == 1:
+                single_q_insts += f"  q{qubits[0]}: {insts}\n"
+            else:
+                multi_q_insts += f"  {qubits}: {insts}\n"
+        instructions = single_q_insts + multi_q_insts
+        return f"<{self.__class__.__name__}({instructions})>"
+
+    def __eq__(self, other):
+        if not isinstance(other, InstructionScheduleMap):
+            return False
+
+        for inst in self.instructions:
+            for qinds in self.qubits_with_instruction(inst):
+                try:
+                    if self._map[inst][_to_tuple(qinds)] != other._map[inst][_to_tuple(qinds)]:
+                        return False
+                except KeyError:
+                    return False
+        return True
+
+
+def _to_tuple(values: Union[int, Iterable[int]]) -> Tuple[int, ...]:
+    """Return the input as a tuple.
+
+    Args:
+        values: An integer, or iterable of integers.
+
+    Returns:
+        The input values as a sorted tuple.
+    """
+    try:
+        return tuple(values)
+    except TypeError:
+        return (values,)
+
+
+def _get_instruction_string(inst: Union[str, Instruction]):
+    if isinstance(inst, str):
+        return inst
+    else:
+        try:
+            return inst.name
+        except AttributeError as ex:
+            raise PulseError(
+                'Input "inst" has no attribute "name". This should be a circuit "Instruction".'
+            ) from ex

qiskit/pulse/instructions/__init__.py

@@ -0,0 +1,69 @@
+# This code is part of Qiskit.
+#
+# (C) Copyright IBM 2020.
+#
+# This code is licensed under the Apache License, Version 2.0. You may
+# obtain a copy of this license in the LICENSE.txt file in the root directory
+# of this source tree or at http://www.apache.org/licenses/LICENSE-2.0.
+#
+# Any modifications or derivative works of this code must retain this
+# copyright notice, and modified files need to carry a notice indicating
+# that they have been altered from the originals.
+
+r"""
+.. _pulse-insts:
+
+===============================================
+Instructions (:mod:`qiskit.pulse.instructions`)
+===============================================
+
+The ``instructions`` module holds the various :obj:`Instruction`\ s which are supported by
+Qiskit Pulse. Instructions have operands, which typically include at least one
+:py:class:`~qiskit.pulse.channels.Channel` specifying where the instruction will be applied.
+
+Every instruction has a duration, whether explicitly included as an operand or implicitly defined.
+For instance, a :py:class:`~qiskit.pulse.instructions.ShiftPhase` instruction can be instantiated
+with operands *phase* and *channel*, for some float ``phase`` and a
+:py:class:`~qiskit.pulse.channels.Channel` ``channel``::
+
+    ShiftPhase(phase, channel)
+
+The duration of this instruction is implicitly zero. On the other hand, the
+:py:class:`~qiskit.pulse.instructions.Delay` instruction takes an explicit duration::
+
+    Delay(duration, channel)
+
+An instruction can be added to a :py:class:`~qiskit.pulse.Schedule`, which is a
+sequence of scheduled Pulse ``Instruction`` s over many channels. ``Instruction`` s and
+``Schedule`` s implement the same interface.
+
+.. autosummary::
+   :toctree: ../stubs/
+
+   Acquire
+   Call
+   Reference
+   Delay
+   Play
+   RelativeBarrier
+   SetFrequency
+   ShiftFrequency
+   SetPhase
+   ShiftPhase
+   Snapshot
+   TimeBlockade
+
+These are all instances of the same base class:
+
+.. autoclass:: Instruction
+"""
+from .acquire import Acquire
+from .delay import Delay
+from .directives import Directive, RelativeBarrier, TimeBlockade
+from .call import Call
+from .instruction import Instruction
+from .frequency import SetFrequency, ShiftFrequency
+from .phase import ShiftPhase, SetPhase
+from .play import Play
+from .snapshot import Snapshot
+from .reference import Reference

qiskit/pulse/instructions/acquire.py

@@ -0,0 +1,150 @@
+# This code is part of Qiskit.
+#
+# (C) Copyright IBM 2020.
+#
+# This code is licensed under the Apache License, Version 2.0. You may
+# obtain a copy of this license in the LICENSE.txt file in the root directory
+# of this source tree or at http://www.apache.org/licenses/LICENSE-2.0.
+#
+# Any modifications or derivative works of this code must retain this
+# copyright notice, and modified files need to carry a notice indicating
+# that they have been altered from the originals.
+
+"""The Acquire instruction is used to trigger the qubit measurement unit and provide
+some metadata for the acquisition process, for example, where to store classified readout data.
+"""
+from typing import Optional, Union, Tuple
+
+from qiskit.circuit import ParameterExpression
+from qiskit.pulse.channels import MemorySlot, RegisterSlot, AcquireChannel
+from qiskit.pulse.configuration import Kernel, Discriminator
+from qiskit.pulse.exceptions import PulseError
+from qiskit.pulse.instructions.instruction import Instruction
+
+
+class Acquire(Instruction):
+    """The Acquire instruction is used to trigger the ADC associated with a particular qubit;
+    e.g. instantiated with AcquireChannel(0), the Acquire command will trigger data collection
+    for the channel associated with qubit 0 readout. This instruction also provides acquisition
+    metadata:
+
+     * the number of cycles during which to acquire (in terms of dt),
+
+     * the register slot to store classified, intermediary readout results,
+
+     * the memory slot to return classified results,
+
+     * the kernel to integrate raw data for each shot, and
+
+     * the discriminator to classify kerneled IQ points.
+    """
+
+    def __init__(
+        self,
+        duration: Union[int, ParameterExpression],
+        channel: AcquireChannel,
+        mem_slot: Optional[MemorySlot] = None,
+        reg_slot: Optional[RegisterSlot] = None,
+        kernel: Optional[Kernel] = None,
+        discriminator: Optional[Discriminator] = None,
+        name: Optional[str] = None,
+    ):
+        """Create a new Acquire instruction.
+
+        Args:
+            duration: Length of time to acquire data in terms of dt.
+            channel: The channel that will acquire data.
+            mem_slot: The classical memory slot in which to store the classified readout result.
+            reg_slot: The fast-access register slot in which to store the classified readout
+                      result for fast feedback.
+            kernel: A ``Kernel`` for integrating raw data.
+            discriminator: A ``Discriminator`` for discriminating kerneled IQ data into 0/1
+                           results.
+            name: Name of the instruction for display purposes.
+        """
+        super().__init__(
+            operands=(duration, channel, mem_slot, reg_slot, kernel, discriminator),
+            name=name,
+        )
+
+    def _validate(self):
+        """Called after initialization to validate instruction data.
+
+        Raises:
+            PulseError: If the input ``channel`` is not type :class:`AcquireChannel`.
+            PulseError: If the input ``mem_slot`` is not type :class:`MemorySlot`.
+            PulseError: If the input ``reg_slot`` is not type :class:`RegisterSlot`.
+            PulseError: When memory slot and register slot are both empty.
+        """
+        if not isinstance(self.channel, AcquireChannel):
+            raise PulseError(f"Expected an acquire channel, got {self.channel} instead.")
+
+        if self.mem_slot and not isinstance(self.mem_slot, MemorySlot):
+            raise PulseError(f"Expected a memory slot, got {self.mem_slot} instead.")
+
+        if self.reg_slot and not isinstance(self.reg_slot, RegisterSlot):
+            raise PulseError(f"Expected a register slot, got {self.reg_slot} instead.")
+
+        if self.mem_slot is None and self.reg_slot is None:
+            raise PulseError("Neither MemorySlots nor RegisterSlots were supplied.")
+
+    @property
+    def channel(self) -> AcquireChannel:
+        """Return the :py:class:`~qiskit.pulse.channels.Channel` that this instruction is
+        scheduled on.
+        """
+        return self.operands[1]
+
+    @property
+    def channels(self) -> Tuple[Union[AcquireChannel, MemorySlot, RegisterSlot]]:
+        """Returns the channels that this schedule uses."""
+        return tuple(self.operands[ind] for ind in (1, 2, 3) if self.operands[ind] is not None)
+
+    @property
+    def duration(self) -> Union[int, ParameterExpression]:
+        """Duration of this instruction."""
+        return self.operands[0]
+
+    @property
+    def kernel(self) -> Kernel:
+        """Return kernel settings."""
+        return self._operands[4]
+
+    @property
+    def discriminator(self) -> Discriminator:
+        """Return discrimination settings."""
+        return self._operands[5]
+
+    @property
+    def acquire(self) -> AcquireChannel:
+        """Acquire channel to acquire data. The ``AcquireChannel`` index maps trivially to
+        qubit index.
+        """
+        return self.channel
+
+    @property
+    def mem_slot(self) -> MemorySlot:
+        """The classical memory slot which will store the classified readout result."""
+        return self.operands[2]
+
+    @property
+    def reg_slot(self) -> RegisterSlot:
+        """The fast-access register slot which will store the classified readout result for
+        fast-feedback computation.
+        """
+        return self.operands[3]
+
+    def is_parameterized(self) -> bool:
+        """Return True iff the instruction is parameterized."""
+        return isinstance(self.duration, ParameterExpression) or super().is_parameterized()
+
+    def __repr__(self) -> str:
+        return "{}({}{}{}{}{}{})".format(
+            self.__class__.__name__,
+            self.duration,
+            ", " + str(self.channel),
+            ", " + str(self.mem_slot) if self.mem_slot else "",
+            ", " + str(self.reg_slot) if self.reg_slot else "",
+            ", " + str(self.kernel) if self.kernel else "",
+            ", " + str(self.discriminator) if self.discriminator else "",
+        )