qiskit 1.3.0__cp39-abi3-win32.whl

qiskit/circuit/library/generalized_gates/permutation.py

@@ -0,0 +1,194 @@
+# This code is part of Qiskit.
+#
+# (C) Copyright IBM 2017, 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.
+
+"""Permutation circuit (the old way to specify permutations, which is required for
+backward compatibility and which will be eventually deprecated) and the permutation
+gate (the new way to specify permutations, allowing a variety of synthesis algorithms).
+"""
+
+from __future__ import annotations
+
+import numpy as np
+
+from qiskit.circuit.quantumcircuit import QuantumCircuit
+from qiskit.circuit.quantumcircuit import Gate
+from qiskit.circuit.exceptions import CircuitError
+from qiskit.utils.deprecation import deprecate_func
+
+
+class Permutation(QuantumCircuit):
+    """An n_qubit circuit that permutes qubits."""
+
+    @deprecate_func(since="1.3", pending=True, additional_msg="Use PermutationGate instead.")
+    def __init__(
+        self,
+        num_qubits: int,
+        pattern: list[int] | np.ndarray | None = None,
+        seed: int | None = None,
+    ) -> None:
+        """Return an n_qubit permutation circuit implemented using SWAPs.
+
+        Args:
+            num_qubits: circuit width.
+            pattern: permutation pattern, describing which qubits occupy the
+                positions 0, 1, 2, etc. after applying the permutation, that
+                is ``pattern[k] = m`` when the permutation maps qubit ``m``
+                to position ``k``. As an example, the pattern ``[2, 4, 3, 0, 1]``
+                means that qubit ``2`` goes to position ``0``, qubit ``4``
+                goes to the position ``1``, etc. The pattern can also be ``None``,
+                in which case a random permutation over ``num_qubits`` is
+                created.
+            seed: random seed in case a random permutation is requested.
+
+        Raises:
+            CircuitError: if permutation pattern is malformed.
+
+        Reference Circuit:
+            .. plot::
+
+               from qiskit.circuit.library import Permutation
+               A = [2,4,3,0,1]
+               circuit = Permutation(5, A)
+               circuit.draw('mpl')
+
+        Expanded Circuit:
+            .. plot::
+
+               from qiskit.circuit.library import Permutation
+               from qiskit.visualization.library import _generate_circuit_library_visualization
+               A = [2,4,3,0,1]
+               circuit = Permutation(5, A)
+               _generate_circuit_library_visualization(circuit.decompose())
+        """
+        if pattern is not None:
+            if sorted(pattern) != list(range(num_qubits)):
+                raise CircuitError(
+                    "Permutation pattern must be some ordering of 0..num_qubits-1 in a list."
+                )
+            pattern = np.array(pattern)
+        else:
+            rng = np.random.default_rng(seed)
+            pattern = np.arange(num_qubits)
+            rng.shuffle(pattern)
+
+        name = "permutation_" + np.array_str(pattern).replace(" ", ",")
+
+        super().__init__(num_qubits, name=name)
+
+        # pylint: disable=cyclic-import
+        from qiskit.synthesis.permutation import synth_permutation_basic
+
+        circuit = synth_permutation_basic(pattern)
+        circuit.name = name
+
+        all_qubits = self.qubits
+        self.append(circuit.to_gate(), all_qubits)
+
+
+class PermutationGate(Gate):
+    """A gate that permutes qubits."""
+
+    def __init__(
+        self,
+        pattern: list[int],
+    ) -> None:
+        """Return a permutation gate.
+
+        Args:
+            pattern: permutation pattern, describing which qubits occupy the
+                positions 0, 1, 2, etc. after applying the permutation, that
+                is ``pattern[k] = m`` when the permutation maps qubit ``m``
+                to position ``k``. As an example, the pattern ``[2, 4, 3, 0, 1]``
+                means that qubit ``2`` goes to position ``0``, qubit ``4``
+                goes to the position ``1``, etc.
+
+        Raises:
+            CircuitError: if permutation pattern is malformed.
+
+        Reference Circuit:
+            .. plot::
+
+                from qiskit.circuit.quantumcircuit import QuantumCircuit
+                from qiskit.circuit.library import PermutationGate
+                A = [2, 4, 3, 0, 1]
+                permutation = PermutationGate(A)
+                circuit = QuantumCircuit(5)
+                circuit.append(permutation, [0, 1, 2, 3, 4])
+                circuit.draw("mpl")
+
+        Expanded Circuit:
+            .. plot::
+
+                from qiskit.circuit.quantumcircuit import QuantumCircuit
+                from qiskit.circuit.library import PermutationGate
+                from qiskit.visualization.library import _generate_circuit_library_visualization
+                A = [2, 4, 3, 0, 1]
+                permutation = PermutationGate(A)
+                circuit = QuantumCircuit(5)
+                circuit.append(permutation, [0, 1, 2, 3, 4])
+
+                _generate_circuit_library_visualization(circuit.decompose())
+        """
+        num_qubits = len(pattern)
+        if sorted(pattern) != list(range(num_qubits)):
+            raise CircuitError(
+                "Permutation pattern must be some ordering of 0..num_qubits-1 in a list."
+            )
+        pattern = np.array(pattern, dtype=np.int32)
+
+        super().__init__(name="permutation", num_qubits=num_qubits, params=[pattern])
+
+    def __array__(self, dtype=None, copy=None):
+        """Return a numpy.array for the Permutation gate."""
+        if copy is False:
+            raise ValueError("unable to avoid copy while creating an array as requested")
+
+        nq = len(self.pattern)
+        mat = np.zeros((2**nq, 2**nq), dtype=dtype)
+
+        for r in range(2**nq):
+            # convert row to bitstring, reverse, apply permutation pattern, reverse again,
+            # and convert to row
+            bit = bin(r)[2:].zfill(nq)[::-1]
+            permuted_bit = "".join([bit[j] for j in self.pattern])
+            pr = int(permuted_bit[::-1], 2)
+            mat[pr, r] = 1
+
+        return mat
+
+    def validate_parameter(self, parameter):
+        """Parameter validation."""
+        return parameter
+
+    @property
+    def pattern(self) -> np.ndarray[bool]:
+        """Returns the permutation pattern defining this permutation."""
+        return self.params[0]
+
+    def inverse(self, annotated: bool = False) -> PermutationGate:
+        """Returns the inverse of the permutation."""
+
+        # pylint: disable=cyclic-import
+        from qiskit.synthesis.permutation.permutation_utils import _inverse_pattern
+
+        return PermutationGate(pattern=_inverse_pattern(self.pattern))
+
+    def _qasm2_decomposition(self):
+        # pylint: disable=cyclic-import
+        from qiskit.synthesis.permutation import synth_permutation_basic
+
+        name = f"permutation__{'_'.join(str(n) for n in self.pattern)}_"
+
+        out = synth_permutation_basic(self.pattern)
+        out.name = name
+
+        return out.to_gate()

qiskit/circuit/library/generalized_gates/rv.py

@@ -0,0 +1,96 @@
+# This code is part of Qiskit.
+#
+# (C) Copyright IBM 2017, 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.
+
+"""Rotation around an arbitrary axis on the Bloch sphere."""
+
+import math
+import numpy
+from qiskit.circuit.gate import Gate
+from qiskit.circuit.exceptions import CircuitError
+
+
+class RVGate(Gate):
+    r"""Rotation around arbitrary rotation axis :math:`\vec{v}` where :math:`\|\vec{v}\|_2` is
+    angle of rotation in radians.
+
+    Can be applied to a :class:`~qiskit.circuit.QuantumCircuit`
+    with the :meth:`~qiskit.circuit.QuantumCircuit.rv` method.
+
+    **Circuit symbol:**
+
+    .. code-block:: text
+
+             ┌─────────────────┐
+        q_0: ┤ RV(v_x,v_y,v_z) ├
+             └─────────────────┘
+
+    **Matrix Representation:**
+
+    .. math::
+
+        \newcommand{\rotationangle}{\frac{\|\vec{v}\|_2}{2}}
+            R(\vec{v}) = e^{-i \vec{v}\cdot\vec{\sigma} / 2} =
+                \begin{pmatrix}
+                    \cos\left(\rotationangle\right)
+                    -i \frac{v_z}{\|\vec{v}\|_2} \sin\left(\rotationangle\right)
+                    & -(i \frac{v_x}{\|\vec{v}\|_2}
+                    + \frac{v_y}{\|\vec{v}\|_2}) \sin\left(\rotationangle\right) \\
+                    -(i \frac{v_x}{\|\vec{v}\|_2}
+                    - \frac{v_y}{\|\vec{v}\|_2}) \sin\left(\rotationangle\right)
+                    & \cos\left(\rotationangle\right)
+                    + i \frac{v_z}{\|\vec{v}\|_2} \sin\left(\rotationangle\right)
+                \end{pmatrix}
+    """
+
+    def __init__(self, v_x: float, v_y: float, v_z: float, basis: str = "U"):
+        """
+        Args:
+            v_x: x-component
+            v_y: y-component
+            v_z: z-component
+            basis: basis (see
+                :class:`~qiskit.synthesis.one_qubit.one_qubit_decompose.OneQubitEulerDecomposer`)
+        """
+        # pylint: disable=cyclic-import
+        from qiskit.synthesis.one_qubit.one_qubit_decompose import OneQubitEulerDecomposer
+
+        super().__init__("rv", 1, [v_x, v_y, v_z])
+        self._decomposer = OneQubitEulerDecomposer(basis=basis)
+
+    def _define(self):
+        try:
+            self.definition = self._decomposer(self.to_matrix())
+        except TypeError as ex:
+            raise CircuitError(
+                f"The {self.name} gate cannot be decomposed with unbound parameters"
+            ) from ex
+
+    def inverse(self, annotated: bool = False):
+        """Invert this gate."""
+        vx, vy, vz = self.params
+        return RVGate(-vx, -vy, -vz)
+
+    def to_matrix(self) -> numpy.ndarray:
+        """Return a numpy.array for the R(v) gate."""
+        v = numpy.asarray(self.params, dtype=float)
+        angle = math.sqrt(v.dot(v))
+        if angle == 0:
+            return numpy.array([[1, 0], [0, 1]])
+        nx, ny, nz = v / angle
+        sin = math.sin(angle / 2)
+        cos = math.cos(angle / 2)
+        return numpy.array(
+            [
+                [cos - 1j * nz * sin, (-ny - 1j * nx) * sin],
+                [(ny - 1j * nx) * sin, cos + 1j * nz * sin],
+            ]
+        )

qiskit/circuit/library/generalized_gates/uc.py

@@ -0,0 +1,213 @@
+# This code is part of Qiskit.
+#
+# (C) Copyright IBM 2020, 2024.
+#
+# 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 structure of the code is based on Emanuel Malvetti's semester thesis at
+# ETH in 2018, which was supervised by Raban Iten and Prof. Renato Renner.
+
+# pylint: disable=invalid-name
+# pylint: disable=missing-param-doc
+# pylint: disable=missing-type-doc
+
+"""Uniformly controlled gates (also called multiplexed gates)."""
+
+from __future__ import annotations
+
+import math
+
+import numpy as np
+
+from qiskit.circuit.gate import Gate
+from qiskit.circuit.library.standard_gates.h import HGate
+from qiskit.quantum_info.operators.predicates import is_unitary_matrix
+from qiskit.circuit.quantumregister import QuantumRegister
+from qiskit.circuit.quantumcircuit import QuantumCircuit
+from qiskit.circuit.exceptions import CircuitError
+from qiskit.exceptions import QiskitError
+from qiskit._accelerate import uc_gate
+
+from .diagonal import Diagonal
+
+_EPS = 1e-10  # global variable used to chop very small numbers to zero
+
+
+class UCGate(Gate):
+    r"""Uniformly controlled gate (also called multiplexed gate).
+
+    These gates can have several control qubits and a single target qubit.
+    If the k control qubits are in the state :math:`|i\rangle` (in the computational basis),
+    a single-qubit unitary :math:`U_i` is applied to the target qubit.
+
+    This gate is represented by a block-diagonal matrix, where each block is a
+    :math:`2\times 2` unitary, that is
+
+    .. math::
+
+        \begin{pmatrix}
+            U_0 & 0 & \cdots & 0 \\
+            0 & U_1 & \cdots & 0 \\
+            \vdots  &     & \ddots & \vdots \\
+            0 & 0   &  \cdots & U_{2^{k-1}}
+        \end{pmatrix}.
+
+    The decomposition is based on Ref. [1].
+
+    **References:**
+
+    [1] Bergholm et al., Quantum circuits with uniformly controlled one-qubit gates (2005).
+        `Phys. Rev. A 71, 052330 <https://journals.aps.org/pra/abstract/10.1103/PhysRevA.71.052330>`__.
+
+    """
+
+    def __init__(self, gate_list: list[np.ndarray], up_to_diagonal: bool = False):
+        r"""
+        Args:
+            gate_list: List of two qubit unitaries :math:`[U_0, ..., U_{2^{k-1}}]`, where each
+                single-qubit unitary :math:`U_i` is given as a :math:`2 \times 2` numpy array.
+            up_to_diagonal: Determines if the gate is implemented up to a diagonal.
+                or if it is decomposed completely (default: False).
+                If the ``UCGate`` :math:`U` is decomposed up to a diagonal :math:`D`, this means
+                that the circuit implements a unitary :math:`U'` such that :math:`D U' = U`.
+
+        Raises:
+            QiskitError: in case of bad input to the constructor
+        """
+        # check input format
+        if not isinstance(gate_list, list):
+            raise QiskitError("The single-qubit unitaries are not provided in a list.")
+        for gate in gate_list:
+            if not gate.shape == (2, 2):
+                raise QiskitError("The dimension of a controlled gate is not equal to (2,2).")
+        if not gate_list:
+            raise QiskitError("The gate list cannot be empty.")
+
+        # Check if number of gates in gate_list is a positive power of two
+        num_contr = math.log2(len(gate_list))
+        if num_contr < 0 or not num_contr.is_integer():
+            raise QiskitError(
+                "The number of controlled single-qubit gates is not a non-negative power of 2."
+            )
+
+        # Check if the single-qubit gates are unitaries
+        for gate in gate_list:
+            if not is_unitary_matrix(gate, _EPS):
+                raise QiskitError("A controlled gate is not unitary.")
+
+        # Create new gate.
+        super().__init__("multiplexer", int(num_contr) + 1, gate_list)
+        self.up_to_diagonal = up_to_diagonal
+
+    def inverse(self, annotated: bool = False) -> Gate:
+        """Return the inverse.
+
+        This does not re-compute the decomposition for the multiplexer with the inverse of the
+        gates but simply inverts the existing decomposition.
+        """
+        if not annotated:
+            inverse_gate = Gate(
+                name=self.name + "_dg", num_qubits=self.num_qubits, params=[]
+            )  # removing the params because arrays are deprecated
+
+            definition = QuantumCircuit(*self.definition.qregs)
+            for inst in reversed(self._definition):
+                definition._append(
+                    inst.replace(operation=inst.operation.inverse(annotated=annotated))
+                )
+
+            definition.global_phase = -self.definition.global_phase
+
+            inverse_gate.definition = definition
+        else:
+            inverse_gate = super().inverse(annotated=annotated)
+        return inverse_gate
+
+    def _get_diagonal(self):
+        # Important: for a control list q_controls = [q[0],...,q_[k-1]] the
+        # diagonal gate is provided in the computational basis of the qubits
+        # q[k-1],...,q[0],q_target, decreasingly ordered with respect to the
+        # significance of the qubit in the computational basis
+        _, diag = self._dec_ucg()
+        return diag
+
+    def _define(self):
+        ucg_circuit, _ = self._dec_ucg()
+        self.definition = ucg_circuit
+
+    def _dec_ucg(self):
+        """
+        Call to create a circuit that implements the uniformly controlled gate. If
+        up_to_diagonal=True, the circuit implements the gate up to a diagonal gate and
+        the diagonal gate is also returned.
+        """
+        diag = np.ones(2**self.num_qubits).tolist()
+        q = QuantumRegister(self.num_qubits, "q")
+        q_controls = q[1:]
+        q_target = q[0]
+        circuit = QuantumCircuit(q, name="uc")
+        # If there is no control, we use the ZYZ decomposition
+        if not q_controls:
+            circuit.unitary(self.params[0], [q])
+            return circuit, diag
+        # If there is at least one control, first,
+        # we find the single qubit gates of the decomposition.
+        (single_qubit_gates, diag) = self._dec_ucg_help()
+        # Now, it is easy to place the C-NOT gates and some Hadamards and Rz(pi/2) gates
+        # (which are absorbed into the single-qubit unitaries) to get back the full decomposition.
+        for i, gate in enumerate(single_qubit_gates):
+            # Absorb Hadamards and Rz(pi/2) gates
+            if i == 0:
+                squ = HGate().to_matrix().dot(gate)
+            elif i == len(single_qubit_gates) - 1:
+                squ = gate.dot(UCGate._rz(np.pi / 2)).dot(HGate().to_matrix())
+            else:
+                squ = (
+                    HGate()
+                    .to_matrix()
+                    .dot(gate.dot(UCGate._rz(np.pi / 2)))
+                    .dot(HGate().to_matrix())
+                )
+            # Add single-qubit gate
+            circuit.unitary(squ, [q_target])
+            # The number of the control qubit is given by the number of zeros at the end
+            # of the binary representation of (i+1)
+            binary_rep = np.binary_repr(i + 1)
+            num_trailing_zeros = len(binary_rep) - len(binary_rep.rstrip("0"))
+            q_contr_index = num_trailing_zeros
+            # Add C-NOT gate
+            if not i == len(single_qubit_gates) - 1:
+                circuit.cx(q_controls[q_contr_index], q_target)
+                circuit.global_phase -= 0.25 * np.pi
+        if not self.up_to_diagonal:
+            # Important: the diagonal gate is given in the computational basis of the qubits
+            # q[k-1],...,q[0],q_target (ordered with decreasing significance),
+            # where q[i] are the control qubits and t denotes the target qubit.
+            diagonal = Diagonal(diag)
+            circuit.append(diagonal, q)
+        return circuit, diag
+
+    def _dec_ucg_help(self):
+        """
+        This method finds the single qubit gate arising in the decomposition of UCGates given in
+        https://arxiv.org/pdf/quant-ph/0410066.pdf.
+        """
+        single_qubit_gates = [gate.astype(complex) for gate in self.params]
+        return uc_gate.dec_ucg_help(single_qubit_gates, self.num_qubits)
+
+    @staticmethod
+    def _rz(alpha):
+        return np.array([[np.exp(1j * alpha / 2), 0], [0, np.exp(-1j * alpha / 2)]])
+
+    def validate_parameter(self, parameter):
+        """Uniformly controlled gate parameter has to be an ndarray."""
+        if isinstance(parameter, np.ndarray):
+            return parameter
+        else:
+            raise CircuitError(f"invalid param type {type(parameter)} in gate {self.name}")

qiskit/circuit/library/generalized_gates/uc_pauli_rot.py

@@ -0,0 +1,164 @@
+# 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 structure of the code is based on Emanuel Malvetti's semester thesis at ETH in 2018,
+# which was supervised by Raban Iten and Prof. Renato Renner.
+
+"""Uniformly controlled Pauli rotations."""
+
+from __future__ import annotations
+
+import math
+
+import numpy as np
+
+from qiskit.circuit.gate import Gate
+from qiskit.circuit.quantumcircuit import QuantumCircuit
+from qiskit.circuit.quantumregister import QuantumRegister
+from qiskit.exceptions import QiskitError
+
+_EPS = 1e-10  # global variable used to chop very small numbers to zero
+
+
+class UCPauliRotGate(Gate):
+    r"""Uniformly controlled Pauli rotations.
+
+    Implements the :class:`.UCGate` for the special case that all unitaries are Pauli rotations,
+    :math:`U_i = R_P(a_i)` where :math:`P \in \{X, Y, Z\}` and :math:`a_i \in \mathbb{R}` is
+    the rotation angle.
+    """
+
+    def __init__(self, angle_list: list[float], rot_axis: str) -> None:
+        r"""
+        Args:
+            angle_list: List of rotation angles :math:`[a_0, ..., a_{2^{k-1}}]`.
+            rot_axis: Rotation axis. Must be either of ``"X"``, ``"Y"`` or ``"Z"``.
+        """
+        self.rot_axes = rot_axis
+        # Check if angle_list has type "list"
+        if not isinstance(angle_list, list):
+            raise QiskitError("The angles are not provided in a list.")
+        # Check if the angles in angle_list are real numbers
+        for angle in angle_list:
+            try:
+                float(angle)
+            except TypeError as ex:
+                raise QiskitError(
+                    "An angle cannot be converted to type float (real angles are expected)."
+                ) from ex
+        num_contr = math.log2(len(angle_list))
+        if num_contr < 0 or not num_contr.is_integer():
+            raise QiskitError(
+                "The number of controlled rotation gates is not a non-negative power of 2."
+            )
+        if rot_axis not in ("X", "Y", "Z"):
+            raise QiskitError("Rotation axis is not supported.")
+        # Create new gate.
+        num_qubits = int(num_contr) + 1
+        super().__init__("ucr" + rot_axis.lower(), num_qubits, angle_list)
+
+    def _define(self):
+        ucr_circuit = self._dec_ucrot()
+        gate = ucr_circuit.to_instruction()
+        q = QuantumRegister(self.num_qubits, "q")
+        ucr_circuit = QuantumCircuit(q)
+        ucr_circuit.append(gate, q[:])
+        self.definition = ucr_circuit
+
+    def _dec_ucrot(self):
+        """
+        Finds a decomposition of a UC rotation gate into elementary gates
+        (C-NOTs and single-qubit rotations).
+        """
+        q = QuantumRegister(self.num_qubits, "q")
+        circuit = QuantumCircuit(q)
+        q_target = q[0]
+        q_controls = q[1:]
+        if not q_controls:  # equivalent to: if len(q_controls) == 0
+            if self.rot_axes == "X":
+                if np.abs(self.params[0]) > _EPS:
+                    circuit.rx(self.params[0], q_target)
+            if self.rot_axes == "Y":
+                if np.abs(self.params[0]) > _EPS:
+                    circuit.ry(self.params[0], q_target)
+            if self.rot_axes == "Z":
+                if np.abs(self.params[0]) > _EPS:
+                    circuit.rz(self.params[0], q_target)
+        else:
+            # First, we find the rotation angles of the single-qubit rotations acting
+            #  on the target qubit
+            angles = self.params.copy()
+            UCPauliRotGate._dec_uc_rotations(angles, 0, len(angles), False)
+            # Now, it is easy to place the C-NOT gates to get back the full decomposition.
+            for i, angle in enumerate(angles):
+                if self.rot_axes == "X":
+                    if np.abs(angle) > _EPS:
+                        circuit.rx(angle, q_target)
+                if self.rot_axes == "Y":
+                    if np.abs(angle) > _EPS:
+                        circuit.ry(angle, q_target)
+                if self.rot_axes == "Z":
+                    if np.abs(angle) > _EPS:
+                        circuit.rz(angle, q_target)
+                # Determine the index of the qubit we want to control the C-NOT gate.
+                # Note that it corresponds
+                # to the number of trailing zeros in the binary representation of i+1
+                if not i == len(angles) - 1:
+                    binary_rep = np.binary_repr(i + 1)
+                    q_contr_index = len(binary_rep) - len(binary_rep.rstrip("0"))
+                else:
+                    # Handle special case:
+                    q_contr_index = len(q_controls) - 1
+                # For X rotations, we have to additionally place some Ry gates around the
+                # C-NOT gates. They change the basis of the NOT operation, such that the
+                # decomposition of for uniformly controlled X rotations works correctly by symmetry
+                # with the decomposition of uniformly controlled Z or Y rotations
+                if self.rot_axes == "X":
+                    circuit.ry(np.pi / 2, q_target)
+                circuit.cx(q_controls[q_contr_index], q_target)
+                if self.rot_axes == "X":
+                    circuit.ry(-np.pi / 2, q_target)
+        return circuit
+
+    @staticmethod
+    def _dec_uc_rotations(angles, start_index, end_index, reversed_dec):
+        """
+        Calculates rotation angles for a uniformly controlled R_t gate with a C-NOT gate at
+        the end of the circuit. The rotation angles of the gate R_t are stored in
+        angles[start_index:end_index]. If reversed_dec == True, it decomposes the gate such that
+        there is a C-NOT gate at the start of the circuit (in fact, the circuit topology for
+        the reversed decomposition is the reversed one of the original decomposition)
+        """
+        interval_len_half = (end_index - start_index) // 2
+        for i in range(start_index, start_index + interval_len_half):
+            if not reversed_dec:
+                angles[i], angles[i + interval_len_half] = UCPauliRotGate._update_angles(
+                    angles[i], angles[i + interval_len_half]
+                )
+            else:
+                angles[i + interval_len_half], angles[i] = UCPauliRotGate._update_angles(
+                    angles[i], angles[i + interval_len_half]
+                )
+        if interval_len_half <= 1:
+            return
+        else:
+            UCPauliRotGate._dec_uc_rotations(
+                angles, start_index, start_index + interval_len_half, False
+            )
+            UCPauliRotGate._dec_uc_rotations(
+                angles, start_index + interval_len_half, end_index, True
+            )
+
+    @staticmethod
+    def _update_angles(angle1, angle2):
+        """Calculate the new rotation angles according to Shende's decomposition."""
+        return (angle1 + angle2) / 2.0, (angle1 - angle2) / 2.0