Add ZX-Calculus functionality with pyzx and OpenQASM import/export ability #2

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open

clausia wants to merge 24 commits into devel from add-pyzx-and-qasm

requirements.txt

-Original file line number
+Diff line change
@@ Expand Up @@
     #pyquil # you need to install the forest-sdk
     #qulacs-gpu
     #qibo <= 0.1.1
+    #optional third party libraries
+    #pyzx

src/tequila/__init__.py

-Original file line number
+Diff line change
@@ Expand Up / @@ -14,6 +14,8 @@ @@
         INSTALLED_SAMPLERS, \
         INSTALLED_SIMULATORS, SUPPORTED_BACKENDS, INSTALLED_BACKENDS, show_available_simulators
     from tequila.wavefunction import QubitWaveFunction
+    from tequila.circuit.qasm import export_open_qasm, import_open_qasm, import_open_qasm_from_file
+    from tequila.circuit.pyzx import convert_to_pyzx, convert_from_pyzx
     import tequila.quantumchemistry as chemistry # shortcut
     from tequila.quantumchemistry import Molecule, MoleculeFromOpenFermion
@@ Expand Down @@

src/tequila/circuit/_gates_impl.py

-Original file line number
+Diff line change
@@ Expand Up / @@ -30,11 +30,16 @@ def control(self): @@
         @property
         def qubits(self):
-            return self._qubits
+            # Set the active qubits
+            if self.control:
+                qubits = self.target + self.control
+            else:
+                qubits = self.target
+            return sorted(tuple(set(qubits)))
         @property
         def max_qubit(self):
-            return self._max_qubit
+            return self.compute_max_qubit()
         def extract_variables(self):
             return []
@@ Expand All @@
             self._target = tuple(list_assignment(target))
             self._control = tuple(list_assignment(control))
             self.finalize()
-            # Set the active qubits
-            if self.control:
-                self._qubits = self.target + self.control
-            else:
-                self._qubits = self.target
-            self._qubits = sorted(tuple(set(self._qubits)))
-            self._max_qubit = self.compute_max_qubit()
         def copy(self):
             return copy.deepcopy(self)
@@ Expand Down Expand Up / @@ -134,12 +132,8 @@ def __eq__(self, other): @@
         def map_qubits(self, qubit_map: dict):
             mapped = copy.deepcopy(self)
             mapped._target = tuple([qubit_map[i] for i in self.target])
-            qubits = mapped._target
             if self.control is not None:
                 mapped._control = tuple([qubit_map[i] for i in self.control])
-                qubits += mapped._control
-            mapped._qubits = sorted(tuple(set(qubits)))
-            mapped._max_qubit = mapped.compute_max_qubit()
             if hasattr(self, "generator"):
                 mapped.generator = self.generator.map_qubits(qubit_map=qubit_map)
             if hasattr(self, "generators"):
@@ Expand Down @@

src/tequila/circuit/compiler.py

-Original file line number
+Diff line change
@@ -1,6 +1,6 @@
     from tequila import TequilaException
     from tequila.circuit.circuit import QCircuit
-    from tequila.circuit.gates import Rx, Ry, H, X, Rz, ExpPauli, CNOT, Phase, T, Z, Y
+    from tequila.circuit.gates import Rx, Ry, H, X, Rz, ExpPauli, CNOT, Phase, T, Z, Y, S, CX
     from tequila.circuit._gates_impl import RotationGateImpl, PhaseGateImpl, QGateImpl, \
         ExponentialPauliGateImpl, TrotterizedGateImpl, PowerGateImpl
     from tequila.utils import to_float
@@ Expand Down Expand Up / @@ -54,7 +54,10 @@ def __init__(self, @@
                      controlled_rotation=False,
                      swap=False,
                      cc_max=False,
-                     gradient_mode=False
+                     gradient_mode=False,
+                     ry_gate=False,
+                     y_gate=False,
+                     ch_gate=False
                      ):
             """
@@ Expand Down Expand Up / @@ -93,6 +96,12 @@ def __init__(self, @@
                 whether or not to break down swap gates into CNOT gates.
             cc_max:
                 whether or not to break down all controlled gates with 2 or more controls.
+            ry_gate:
+                whether or not to break down all rotational y gates
+            y_gate:
+                whether or not to break down all y gates
+            ch_gate:
+                whether or not to break down all controlled-H gates
             """
             self.multitarget = multitarget
             self.multicontrol = multicontrol
@@ Expand All / @@ -111,6 +120,9 @@ def __init__(self, @@
             self.swap = swap
             self.cc_max = cc_max
             self.gradient_mode = gradient_mode
+            self.ry_gate = ry_gate
+            self.y_gate = y_gate
+            self.ch_gate = ch_gate
         def __call__(self, objective: typing.Union[Objective, QCircuit, ExpectationValueImpl], variables=None, *args,
                      **kwargs):
@@ Expand Down Expand Up @@
                     cg = compile_power_gate(gate=cg)
                 if self.phase:
                     cg = compile_phase(gate=cg)
+                if self.ch_gate:
+                    cg = compile_ch(gate=cg)
+                if self.y_gate:
+                    cg = compile_y(gate=cg)
+                if self.ry_gate:
+                    cg = compile_ry(gate=cg, controlled_rotation=self.controlled_rotation)
                 if controlled:
                     if self.cc_max:
                         cg = compile_to_cc(gate=cg)
@@ Expand Down Expand Up @@
             return compile_exponential_pauli_gate(result)
         else:
             return result
+    @compiler
+    def compile_ry(gate: RotationGateImpl, controlled_rotation: bool = False) -> QCircuit:
+        """
+        Compile Ry gates into Rx and Rz.
+        Parameters
+        ----------
+        gate:
+            the gate.
+        controlled_rotation:
+            determines if the decomposition of the controlled-Ry gate will be performed in compile_controlled_rotation,
+            if not, decomposition will be performed here
+        Returns
+        -------
+        QCircuit, the result of compilation.
+        """
+        if gate.name.lower() == "ry":
+            if not (gate.is_controlled() and controlled_rotation):
+                return Rz(target=gate.target, control=None, angle=-numpy.pi / 2) \
+                       + Rx(target=gate.target, control=gate.control, angle=gate.parameter) \
+                       + Rz(target=gate.target, control=None, angle=numpy.pi / 2)
+        return QCircuit.wrap_gate(gate)
+    @compiler
+    def compile_y(gate) -> QCircuit:
+        """
+        Compile Y gates into X and Rz.
+        Parameters
+        ----------
+        gate:
+            the gate.
+        Returns
+        -------
+        QCircuit, the result of compilation.
+        """
+        if gate.name.lower() == "y":
+            return Rz(target=gate.target, control=None, angle=-numpy.pi / 2) \
+                   + X(target=gate.target, control=gate.control, power=gate.parameter if gate.is_parametrized() else None) \
+                   + Rz(target=gate.target, control=None, angle=numpy.pi / 2)
+        else:
+            return QCircuit.wrap_gate(gate)
+    @compiler
+    def compile_ch(gate: QGateImpl) -> QCircuit:
+        """
+        Compile CH gates into its equivalent:
+            CH = Ry(0.25pi) CZ Ry(-0.25pi)
+        Parameters
+        ----------
+        gate:
+            the gate.
+        Returns
+        -------
+        QCircuit, the result of compilation.
+        """
+        if gate.name.lower() == "h" and gate.is_controlled():
+            return Ry(target=gate.target, control=None, angle=-numpy.pi / 4) \
+                   + Z(target=gate.target, control=gate.control, power=gate.parameter if gate.is_parametrized() else None) \
+                   + Ry(target=gate.target, control=None, angle=numpy.pi / 4)
+        else:
+            return QCircuit.wrap_gate(gate)

src/tequila/circuit/gates.py

-Original file line number
+Diff line change
@@ Expand Up @@
         return Rz(target=target, control=control, angle=angle)
+    def U(theta, phi, lambd, target: typing.Union[list, int], control: typing.Union[list, int] = None) -> QCircuit:
+        """
+        Notes
+        ----------
+        Convenient gate, one of the abstract gates defined by OpenQASM.
+        .. math::
+            U(\\theta, \\phi, \\lambda) = R_z(\\phi)R_x(-\\pi/2)R_z(\\theta)R_x(\\pi/2)R_z(\\lambda)
+            U(\\theta, \\phi, \\lambda) = \\begin{pmatrix}
+                                                e^{-i \\frac{\\phi}{2}} & 0 \\\\
+& e^{i \\frac{\\phi}{2}}
+                                            \\end{pmatrix}
+                                            \\begin{pmatrix}
+                                                \\cos{-\\frac{\\pi}{4}} & -i \\sin{-\\frac{\\pi}{4}} \\\\
+                                                -i \\sin{-\\frac{\\pi}{4}} & \\cos{-\\frac{\\pi}{4}}
+                                            \\end{pmatrix}
+                                            \\begin{pmatrix}
+                                                e^{-i \\frac{\\theta}{2}} & 0 \\\\
+& e^{i \\frac{\\theta}{2}}
+                                            \\end{pmatrix}
+                                            \\begin{pmatrix}
+                                                \\cos{\\frac{\\pi}{4}} & -i \\sin{\\frac{\\pi}{4}} \\\\
+                                                -i \\sin{\\frac{\\pi}{4}} & \\cos{\\frac{\\pi}{4}}
+                                            \\end{pmatrix}
+                                            \\begin{pmatrix}
+                                                e^{-i \\frac{\\lambda}{2}} & 0 \\\\
+& e^{i \\frac{\\lambda}{2}}
+                                            \\end{pmatrix}
+            U(\\theta, \\phi, \\lambda) = \\begin{pmatrix}
+                                            \\cos{\\frac{\\theta}{2}} &
+                                            -e^{i \\lambda} \\sin{\\frac{\\theta}{2}} \\\\
+                                            e^{i \\phi} \\sin{\\frac{\\theta}{2}} &
+                                            e^{i (\\phi+\\lambda)} \\cos{\\frac{\\theta}{2}}
+                                          \\end{pmatrix}
+        Parameters
+        ----------
+        theta
+            first parameter angle
+        phi
+            second parameter angle
+        lamnd
+            third parameter angle
+        target
+            int or list of int
+        control
+            int or list of int
+        Returns
+        -------
+        QCircuit object
+        """
+        theta = assign_variable(theta)
+        phi = assign_variable(phi)
+        lambd = assign_variable(lambd)
+        pi_half = assign_variable(np.pi / 2)
+        return Rz(angle=lambd, target=target, control=control) + \
+               Rx(angle=pi_half, target=target, control=control) + \
+               Rz(angle=theta, target=target, control=control) + \
+               Rx(angle=-pi_half, target=target, control=control) + \
+               Rz(angle=phi, target=target, control=control)
+    def u1(lambd, target: typing.Union[list, int], control: typing.Union[list, int] = None) -> QCircuit:
+        """
+        Notes
+        ----------
+        Convenient gate, one of the abstract gates defined by Quantum Experience Standard Header.
+        Changes the phase of a carrier without applying any pulses.
+        .. math::
+            from OpenQASM 2.0 specification:
+                u1(\\lambda) \\sim U(0, 0, \\lambda) = R_z(\\lambda) = e^{-i\\frac{\\lambda}{2} \\sigma_{z}}
+            also is equal to:
+                u1(\\lambda) = \\begin{pmatrix} 1 & 0 \\\\ 0 & e^{i\\lambda} \\end{pmatrix}
+            which is the Tequila Phase gate:
+                u1(\\lambda) = Phase(\\lambda)
+        Parameters
+        ----------
+        lambd
+            parameter angle
+        target
+            int or list of int
+        control
+            int or list of int
+        Returns
+        -------
+        QCircuit object
+        """
+        return Phase(phi=lambd, target=target, control=control)
+    def u2(phi, lambd, target: typing.Union[list, int], control: typing.Union[list, int] = None) -> QCircuit:
+        """
+        Notes
+        ----------
+        Convenient gate, one of the abstract gates defined by Quantum Experience Standard Header.
+        Uses a single \\pi/2-pulse.
+        .. math::
+            u2(\\phi, \\lambda) = U(\\pi/2, \\phi, \\lambda) = R_z(\\phi + \\pi/2)R_x(\\pi/2)R_z(\\lambda - \\pi/2)
+            u2(\\phi, \\lambda) = \\frac{1}{\\sqrt{2}}
+                                  \\begin{pmatrix}
+& -e^{i\\lambda} \\\\
+                                        e^{i\\phi} & e^{i(\\phi+\\lambda)}
+                                  \\end{pmatrix}
+        Parameters
+        ----------
+        phi
+            first parameter angle
+        lambd
+            second parameter angle
+        target
+            int or list of int
+        control
+            int or list of int
+        Returns
+        -------
+        QCircuit object
+        """
+        return U(theta=np.pi/2, phi=phi, lambd=lambd, target=target, control=control)
+    def u3(theta, phi, lambd, target: typing.Union[list, int], control: typing.Union[list, int] = None) -> QCircuit:
+        """
+        Notes
+        ----------
+        Convenient gate, one of the abstract gates defined by Quantum Experience Standard Header
+        The most general single-qubit gate.
+        Uses a pair of \\pi/2-pulses.
+        .. math::
+            u3(\\theta, \\phi, \\lambda) = U(\\theta, \\phi, \\lambda)
+                                         = \\begin{pmatrix}
+                                                \\cos{\\frac{\\5theta}{2}} &
+                                                -e^{i \\lambda} \\sin{\\frac{\\theta}{2}} \\\\
+                                                e^{i \\phi} \\sin{\\frac{\\theta}{2}} &
+                                                e^{i (\\phi+\\lambda)} \\cos{\\frac{\\theta}{2}}
+                                           \\end{pmatrix}
+        Parameters
+        ----------
+        theta
+            first parameter angle
+        phi
+            second parameter angle
+        lambd
+            third parameter angle
+        target
+            int or list of int
+        control
+            int or list of int
+        Returns
+        -------
+        QCircuit object
+        """
+        return U(theta=theta, phi=phi, lambd=lambd, target=target, control=control)
     if __name__ == "__main__":
         G = CRx(1, 0, 2.0)
@@ Expand Down @@

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Add ZX-Calculus functionality with pyzx and OpenQASM import/export ability #2

Uh oh!

Diff view

Diff view

There are no files selected for viewing

Uh oh!

Uh oh!

Uh oh!

Add ZX-Calculus functionality with pyzx and OpenQASM import/export ability #2

Are you sure you want to change the base?

Uh oh!

Add ZX-Calculus functionality with pyzx and OpenQASM import/export ability #2

Uh oh!

Uh oh!

Diff view

Diff view

There are no files selected for viewing

Uh oh!

Uh oh!

Uh oh!