Feat multi-observable for IBM and diagonal observables

examples/notebooks/3_Expectation_value_of_observables.ipynb

@@ -554,6 +554,26 @@
     ")\n",
     "print(results)"
    ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "ad9fed49",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "#TODO: show examples with multi observable expectation value\n",
+    "\n",
+    "#TODO: maybe show also the set of commuting pauli string when those observables are decomposed."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "a6f5f27c",
+   "metadata": {},
+   "outputs": [],
+   "source": []
   }
  ],
  "metadata": {

examples/scripts/aws_execution_trials.py

@@ -42,7 +42,8 @@
 circuit.add(SWAP(2, 0))
 circuit.add(CNOT(0, 2))
 circuit.add(Ry(3.14 / 2, 2))
-circuit.add(S(1))
+circuit.add(S(0))
+circuit.add(S_dagger(1))
 circuit.add(H(3))
 circuit.add(CNOT(1, 2))
 circuit.add(Rx(3.14, 1))

examples/scripts/circuit_to_qasm.py

@@ -15,7 +15,8 @@
 circuit.add(SWAP(2, 0))
 circuit.add(CNOT(0, 2))
 circuit.add(Ry(3.14 / 2, 2))
-circuit.add(S(1))
+circuit.add(S(0))
+circuit.add(S_dagger(1))
 circuit.add(H(3))
 circuit.add(CNOT(1, 2))
 circuit.add(Rx(3.14, 1))

examples/scripts/native_gates_examples.py

@@ -9,7 +9,7 @@
 # Declaration of the circuit with the right size
 circuit = QCircuit(3, label="Test native gates")
 # Constructing the circuit by adding gates and measurements
-circuit.add([H(0), X(1), Y(2), Z(0), S(1), T(0)])
+circuit.add([H(0), X(1), Y(2), Z(0), S(0), S_dagger(1), T(0)])
 circuit.add([Rx(1.2324, 2), Ry(-2.43, 0), Rz(1.04, 1), Rk(-1, 1), P(-323, 2)])
 circuit.add(U(1.2, 2.3, 3.4, 2))
 circuit.add(SWAP(2, 0))

examples/scripts/observable_job.py

@@ -71,3 +71,5 @@
     ],
 )
 print(results)
+
+# TODO: show examples of multi observable run

mpqp/all.py

@@ -49,6 +49,7 @@
     Ry,
     Rz,
     S,
+    S_dagger,
     T,
     U,
     UnitaryMatrix,

mpqp/core/circuit.py

@@ -760,29 +760,29 @@ def inverse(self) -> QCircuit:
             The inverse circuit.
 
         Examples:
-            >>> c1 = QCircuit([S(0), CZ(0,1), H(1), Ry(4.56, 1)])
+            >>> c1 = QCircuit([T(0), CZ(0,1), H(1), Ry(4.56, 1)])
             >>> print(c1)  # doctest: +NORMALIZE_WHITESPACE
                  ┌───┐
-            q_0: ┤ S ├─■──────────────────
+            q_0: ┤ T ├─■──────────────────
                  └───┘ │ ┌───┐┌──────────┐
             q_1: ──────■─┤ H ├┤ Ry(4.56) ├
                          └───┘└──────────┘
             >>> print(c1.inverse())  # doctest: +NORMALIZE_WHITESPACE
                                       ┌────┐
-            q_0: ───────────────────■─┤ S† ├
+            q_0: ───────────────────■─┤ T† ├
                  ┌───────────┐┌───┐ │ └────┘
             q_1: ┤ Ry(-4.56) ├┤ H ├─■───────
                  └───────────┘└───┘
-             >>> c2 = QCircuit([S(0), CRk(2, 0, 1), Barrier(), H(1), Ry(4.56, 1)])
+             >>> c2 = QCircuit([T(0), CRk(2, 0, 1), Barrier(), H(1), Ry(4.56, 1)])
             >>> print(c2)  # doctest: +NORMALIZE_WHITESPACE
                  ┌───┐          ░
-            q_0: ┤ S ├─■────────░──────────────────
+            q_0: ┤ T ├─■────────░──────────────────
                  └───┘ │P(π/2)  ░ ┌───┐┌──────────┐
             q_1: ──────■────────░─┤ H ├┤ Ry(4.56) ├
                                 ░ └───┘└──────────┘
             >>> print(c2.inverse())  # doctest: +NORMALIZE_WHITESPACE
                                     ░           ┌────┐
-            q_0: ───────────────────░──■────────┤ S† ├
+            q_0: ───────────────────░──■────────┤ T† ├
                  ┌───────────┐┌───┐ ░  │P(-π/2) └────┘
             q_1: ┤ Ry(-4.56) ├┤ H ├─░──■──────────────
                  └───────────┘└───┘ ░
@@ -906,7 +906,7 @@ def measurements(self) -> list[Measure]:
             ... ])
             >>> circuit.measurements  # doctest: +NORMALIZE_WHITESPACE
             [BasisMeasure(shots=1000),
-            ExpectationMeasure(Observable(array([[1.+0.j, 0.+0.j], [0.+0.j, 1.+0.j]], dtype=complex64)), [1], shots=1000)]
+            ExpectationMeasure([Observable(array([[1.+0.j, 0.+0.j], [0.+0.j, 1.+0.j]], dtype=complex64))], [1], shots=1000)]
 
         """
         return [inst for inst in self.instructions if isinstance(inst, Measure)]

mpqp/core/instruction/gates/__init__.py

@@ -20,6 +20,7 @@
     Ry,
     Rz,
     S,
+    S_dagger,
     T,
     U,
     X,

mpqp/core/instruction/gates/native_gates.py

@@ -378,8 +378,10 @@ class OneQubitNoParamGate(SingleQubitGate, NoParameterGate, SimpleClassReprABC):
         target: Index referring to the qubits on which the gate will be applied.
     """
 
-    def __init__(self, target: int):
-        SingleQubitGate.__init__(self, target, type(self).__name__)
+    def __init__(self, target: int, label: Optional[str] = None):
+        SingleQubitGate.__init__(
+            self, target, type(self).__name__ if label is None else label
+        )
 
 
 class Id(OneQubitNoParamGate, InvolutionGate):
@@ -767,6 +769,56 @@ def __init__(self, target: int):
         super().__init__(target)
         self.matrix = np.array([[1, 0], [0, 1j]])
 
+    def inverse(self) -> Gate:
+        return S_dagger(self.targets[0])
+
+
+class S_dagger(OneQubitNoParamGate):
+    r"""One qubit S adjoint gate. It's equivalent to ``P(-pi/2)``.
+
+    `\begin{pmatrix}1&0\\0&-i\end{pmatrix}`
+
+    Args:
+        target: Index referring to the qubit on which the gate will be applied.
+
+    Example:
+        >>> pprint(S_dagger(0).to_matrix())
+        [[1, 0  ],
+         [0, -1j]]
+
+    """
+
+    @classproperty
+    def braket_gate(cls):
+        from braket.circuits import gates
+
+        return gates.Si
+
+    @classproperty
+    def qiskit_gate(cls):
+        from qiskit.circuit.library import SdgGate
+
+        return SdgGate
+
+    @classproperty
+    def cirq_gate(cls):
+        from cirq.ops.common_gates import ZPowGate
+
+        return ZPowGate(exponent=-1 / 2)
+
+    qlm_aqasm_keyword = "DAG(S)"
+    qiskit_string = "sdg"
+
+    def __init__(self, target: int):
+        super().__init__(target, "S†")
+        self.matrix = np.array([[1, 0], [0, -1j]])
+
+    def __repr__(self):
+        return f"{type(self).__name__}({self.targets[0]})"
+
+    def inverse(self) -> Gate:
+        return S(self.targets[0])
+
 
 class T(OneQubitNoParamGate):
     r"""One qubit T gate. It is also referred to as the `\pi/4` gate because it