fix: type errors

Author: MoHermes

mpqp/core/instruction/measurement/expectation_value.py

@@ -12,6 +12,8 @@ class to define your observable, and a :class:`ExpectationMeasure` to perform
 
 import numpy as np
 import numpy.typing as npt
+from typeguard import typechecked
+
 from mpqp.core.instruction.gates.native_gates import SWAP
 from mpqp.core.instruction.measurement.measure import Measure
 from mpqp.core.instruction.measurement.pauli_string import (
@@ -23,7 +25,6 @@ class to define your observable, and a :class:`ExpectationMeasure` to perform
 from mpqp.tools.errors import NumberQubitsError
 from mpqp.tools.generics import Matrix, OneOrMany
 from mpqp.tools.maths import is_diagonal, is_hermitian, is_power_of_two
-from typeguard import typechecked
 
 if TYPE_CHECKING:
     from braket.circuits.observables import Hermitian
@@ -126,7 +127,7 @@ def matrix(self) -> Matrix:
         """The matrix representation of the observable."""
         if self._matrix is None:
             if self.is_diagonal and self._diag_elements is not None:
-                self._matrix = np.diag(self._diag_elements)
+                self._matrix = np.diag(np.array(self._diag_elements, dtype=np.float64))
             else:
                 self._matrix = self.pauli_string.to_matrix()
         matrix = copy.deepcopy(self._matrix).astype(np.complex64)
@@ -138,7 +139,7 @@ def pauli_string(self) -> PauliString:
         if self._pauli_string is None:
             if self.is_diagonal:
                 self._pauli_string = PauliString.from_diagonal_elements(
-                    self.diagonal_elements
+                    np.array(self.diagonal_elements, dtype=np.float64)
                 )
             else:
                 self._pauli_string = PauliString.from_matrix(self.matrix)
@@ -150,7 +151,7 @@ def diagonal_elements(self) -> npt.NDArray[np.float32]:
         """The diagonal elements of the matrix representing the observable (diagonal or not)."""
         if self._diag_elements is None:
             self._diag_elements = np.diagonal(self.matrix)
-        return copy.deepcopy(self._diag_elements).real
+        return copy.deepcopy(np.array(self._diag_elements, dtype=np.float32)).real
 
     @matrix.setter
     def matrix(self, matrix: Matrix):

mpqp/tools/obs_decomposition.py

@@ -404,7 +404,7 @@ def generate_and_explore_node_diagonal_case(
 
 @typechecked
 def decompose_diagonal_observable_ptdr(
-    diag_elements: list[Real] | npt.NDArray[np.float32], print_progression: bool = False
+    diag_elements: list[Real] | npt.NDArray[np.float64], print_progression: bool = False
 ) -> PauliString:
     """Decomposes a diagonal observable into a Pauli string representation.
 
@@ -505,7 +505,7 @@ def compute_coefficients_walsh(
 
 @typechecked
 def decompose_diagonal_observable_walsh_hadamard(
-    diag_elements: list[Real] | npt.NDArray[np.float32],
+    diag_elements: list[Real] | npt.NDArray[np.float64],
 ) -> PauliString:
     """Decomposes the observable represented by the diagonal elements into a
     Pauli string using the Walsh-Hadamard transform.