qadence 1.7.8__py3-none-any.whl → 1.9.0__py3-none-any.whl

qadence/operations/parametric.py

@@ -15,6 +15,7 @@ from qadence.blocks.utils import (
     add,  # noqa
     chain,
 )
+from qadence.noise import NoiseHandler
 from qadence.parameters import (
     Parameter,
     ParamMap,
@@ -32,10 +33,15 @@ class PHASE(ParametricBlock):
 
     name = OpName.PHASE
 
-    def __init__(self, target: int, parameter: Parameter | TNumber | sympy.Expr | str):
+    def __init__(
+        self,
+        target: int,
+        parameter: Parameter | TNumber | sympy.Expr | str,
+        noise: NoiseHandler | None = None,
+    ) -> None:
         self.parameters = ParamMap(parameter=parameter)
         self.generator = I(target) - Z(target)
-        super().__init__((target,))
+        super().__init__((target,), noise=noise)
 
     @classmethod
     def num_parameters(cls) -> int:
@@ -56,13 +62,18 @@ class RX(ParametricBlock):
 
     name = OpName.RX
 
-    def __init__(self, target: int, parameter: Parameter | TParameter | ParamMap):
+    def __init__(
+        self,
+        target: int,
+        parameter: Parameter | TParameter | ParamMap,
+        noise: NoiseHandler | None = None,
+    ) -> None:
         # TODO: should we give them more meaningful names? like 'angle'?
         self.parameters = (
             parameter if isinstance(parameter, ParamMap) else ParamMap(parameter=parameter)
         )
         self.generator = X(target)
-        super().__init__((target,))
+        super().__init__((target,), noise=noise)
 
     @classmethod
     def num_parameters(cls) -> int:
@@ -84,12 +95,17 @@ class RY(ParametricBlock):
 
     name = OpName.RY
 
-    def __init__(self, target: int, parameter: Parameter | TParameter | ParamMap):
+    def __init__(
+        self,
+        target: int,
+        parameter: Parameter | TParameter | ParamMap,
+        noise: NoiseHandler | None = None,
+    ) -> None:
         self.parameters = (
             parameter if isinstance(parameter, ParamMap) else ParamMap(parameter=parameter)
         )
         self.generator = Y(target)
-        super().__init__((target,))
+        super().__init__((target,), noise=noise)
 
     @classmethod
     def num_parameters(cls) -> int:
@@ -111,12 +127,17 @@ class RZ(ParametricBlock):
 
     name = OpName.RZ
 
-    def __init__(self, target: int, parameter: Parameter | TParameter | ParamMap):
+    def __init__(
+        self,
+        target: int,
+        parameter: Parameter | TParameter | ParamMap,
+        noise: NoiseHandler | None = None,
+    ) -> None:
         self.parameters = (
             parameter if isinstance(parameter, ParamMap) else ParamMap(parameter=parameter)
         )
         self.generator = Z(target)
-        super().__init__((target,))
+        super().__init__((target,), noise=noise)
 
     @classmethod
     def num_parameters(cls) -> int:
@@ -147,10 +168,11 @@ class U(ParametricBlock):
         phi: Parameter | TParameter,
         theta: Parameter | TParameter,
         omega: Parameter | TParameter,
-    ):
+        noise: NoiseHandler | None = None,
+    ) -> None:
         self.parameters = ParamMap(phi=phi, theta=theta, omega=omega)
         self.generator = chain(Z(target), Y(target), Z(target))
-        super().__init__((target,))
+        super().__init__((target,), noise=noise)
 
     @classmethod
     def num_parameters(cls) -> int:

qadence/operations/primitive.py

@@ -18,6 +18,7 @@ from qadence.blocks.utils import (
     chain,
     kron,
 )
+from qadence.noise import NoiseHandler
 from qadence.parameters import (
     Parameter,
 )
@@ -31,8 +32,8 @@ class X(PrimitiveBlock):
 
     name = OpName.X
 
-    def __init__(self, target: int):
-        super().__init__((target,))
+    def __init__(self, target: int, noise: NoiseHandler | None = None) -> None:
+        super().__init__((target,), noise=noise)
 
     @property
     def generator(self) -> AbstractBlock:
@@ -52,8 +53,8 @@ class Y(PrimitiveBlock):
 
     name = OpName.Y
 
-    def __init__(self, target: int):
-        super().__init__((target,))
+    def __init__(self, target: int, noise: NoiseHandler | None = None) -> None:
+        super().__init__((target,), noise=noise)
 
     @property
     def generator(self) -> AbstractBlock:
@@ -73,8 +74,8 @@ class Z(PrimitiveBlock):
 
     name = OpName.Z
 
-    def __init__(self, target: int):
-        super().__init__((target,))
+    def __init__(self, target: int, noise: NoiseHandler | None = None) -> None:
+        super().__init__((target,), noise=noise)
 
     @property
     def generator(self) -> AbstractBlock:
@@ -94,8 +95,8 @@ class I(PrimitiveBlock):
 
     name = OpName.I
 
-    def __init__(self, target: int):
-        super().__init__((target,))
+    def __init__(self, target: int, noise: NoiseHandler | None = None) -> None:
+        super().__init__((target,), noise=noise)
 
     def __ixor__(self, other: AbstractBlock | int) -> AbstractBlock:
         if not isinstance(other, AbstractBlock):
@@ -113,7 +114,7 @@ class I(PrimitiveBlock):
 
     @property
     def generator(self) -> AbstractBlock:
-        return I(*self.qubit_support)
+        return I(self.qubit_support[0])
 
     @property
     def eigenvalues_generator(self) -> Tensor:
@@ -137,8 +138,9 @@ class Projector(ProjectorBlock):
         ket: str,
         bra: str,
         qubit_support: int | tuple[int, ...],
-    ):
-        super().__init__(ket=ket, bra=bra, qubit_support=qubit_support)
+        noise: NoiseHandler | None = None,
+    ) -> None:
+        super().__init__(ket=ket, bra=bra, qubit_support=qubit_support, noise=noise)
 
     @property
     def generator(self) -> None:
@@ -154,8 +156,13 @@ class N(Projector):
 
     name = OpName.N
 
-    def __init__(self, target: int, state: str = "1"):
-        super().__init__(ket=state, bra=state, qubit_support=(target,))
+    def __init__(
+        self,
+        target: int,
+        state: str = "1",
+        noise: NoiseHandler | None = None,
+    ) -> None:
+        super().__init__(ket=state, bra=state, qubit_support=(target,), noise=noise)
 
     @property
     def generator(self) -> None:
@@ -175,9 +182,13 @@ class S(PrimitiveBlock):
 
     name = OpName.S
 
-    def __init__(self, target: int):
+    def __init__(
+        self,
+        target: int,
+        noise: NoiseHandler | None = None,
+    ) -> None:
         self.generator = I(target) - Z(target)
-        super().__init__((target,))
+        super().__init__((target,), noise=noise)
 
     @property
     def eigenvalues_generator(self) -> Tensor:
@@ -188,7 +199,7 @@ class S(PrimitiveBlock):
         return tensor([1, 1j], dtype=cdouble)
 
     def dagger(self) -> SDagger:
-        return SDagger(*self.qubit_support)
+        return SDagger(self.qubit_support[0])
 
 
 class SDagger(PrimitiveBlock):
@@ -196,9 +207,13 @@ class SDagger(PrimitiveBlock):
 
     name = OpName.SDAGGER
 
-    def __init__(self, target: int):
+    def __init__(
+        self,
+        target: int,
+        noise: NoiseHandler | None = None,
+    ) -> None:
         self.generator = I(target) - Z(target)
-        super().__init__((target,))
+        super().__init__((target,), noise=noise)
 
     @property
     def eigenvalues_generator(self) -> Tensor:
@@ -209,7 +224,7 @@ class SDagger(PrimitiveBlock):
         return tensor([1, -1j], dtype=cdouble)
 
     def dagger(self) -> S:
-        return S(*self.qubit_support)
+        return S(self.qubit_support[0])
 
 
 class H(PrimitiveBlock):
@@ -217,9 +232,13 @@ class H(PrimitiveBlock):
 
     name = OpName.H
 
-    def __init__(self, target: int):
+    def __init__(
+        self,
+        target: int,
+        noise: NoiseHandler | None = None,
+    ) -> None:
         self.generator = (1 / np.sqrt(2)) * (X(target) + Z(target) - np.sqrt(2) * I(target))
-        super().__init__((target,))
+        super().__init__((target,), noise=noise)
 
     @property
     def eigenvalues_generator(self) -> Tensor:
@@ -275,9 +294,13 @@ class T(PrimitiveBlock):
 
     name = OpName.T
 
-    def __init__(self, target: int):
+    def __init__(
+        self,
+        target: int,
+        noise: NoiseHandler | None = None,
+    ) -> None:
         self.generator = I(target) - Z(target)
-        super().__init__((target,))
+        super().__init__((target,), noise)
 
     @property
     def eigenvalues_generator(self) -> Tensor:
@@ -292,7 +315,7 @@ class T(PrimitiveBlock):
         return 1
 
     def dagger(self) -> TDagger:
-        return TDagger(*self.qubit_support)
+        return TDagger(self.qubit_support[0])
 
 
 class TDagger(PrimitiveBlock):
@@ -301,9 +324,13 @@ class TDagger(PrimitiveBlock):
     # FIXME: this gate is not support by any backend
     name = "T_dagger"
 
-    def __init__(self, target: int):
+    def __init__(
+        self,
+        target: int,
+        noise: NoiseHandler | None = None,
+    ) -> None:
         self.generator = I(target) - Z(target)
-        super().__init__((target,))
+        super().__init__((target,), noise)
 
     @property
     def eigenvalues_generator(self) -> Tensor:
@@ -318,7 +345,7 @@ class TDagger(PrimitiveBlock):
         return 1
 
     def dagger(self) -> T:
-        return T(*self.qubit_support)
+        return T(self.qubit_support[0])
 
 
 class SWAP(PrimitiveBlock):
@@ -326,7 +353,12 @@ class SWAP(PrimitiveBlock):
 
     name = OpName.SWAP
 
-    def __init__(self, control: int, target: int) -> None:
+    def __init__(
+        self,
+        control: int,
+        target: int,
+        noise: NoiseHandler | None = None,
+    ) -> None:
         a11 = 0.5 * (Z(control) - I(control))
         a22 = -0.5 * (Z(target) + I(target))
         a12 = 0.5 * (chain(X(control), Z(control)) + X(control))
@@ -334,7 +366,7 @@ class SWAP(PrimitiveBlock):
         self.generator = (
             kron(-1.0 * a22, a11) + kron(-1.0 * a11, a22) + kron(a12, a21) + kron(a21, a12)
         )
-        super().__init__((control, target))
+        super().__init__((control, target), noise=noise)
 
     @property
     def eigenvalues_generator(self) -> Tensor:

qadence/overlap.py

@@ -102,9 +102,6 @@ def _select_overlap_method(
             return overlap
 
     elif method == OverlapMethod.SWAP_TEST:
-        if backend == BackendName.BRAKET:
-            raise ValueError("SWAP test method is not supported by the Braket backend.")
-
         n_qubits = bra_circuit.n_qubits
 
         # shift qubit support of bra and ket circuit blocks
@@ -128,9 +125,6 @@ def _select_overlap_method(
             return overlap
 
     elif method == OverlapMethod.HADAMARD_TEST:
-        if backend == BackendName.BRAKET:
-            raise ValueError("Hadamard test method is not supported by the Braket backend.")
-
         n_qubits = bra_circuit.n_qubits
 
         # construct controlled bra and ket blocks

qadence/parameters.py

@@ -228,7 +228,7 @@ def sympy_to_numeric(expr: Basic) -> TNumber:
         return float(expr)
 
 
-def evaluate(expr: Expr, values: dict = {}, as_torch: bool = False) -> TNumber | Tensor:
+def evaluate(expr: Expr, values: dict | None = None, as_torch: bool = False) -> TNumber | Tensor:
     """
     Arguments:
 
@@ -260,7 +260,8 @@ def evaluate(expr: Expr, values: dict = {}, as_torch: bool = False) -> TNumber |
     """
     res: Basic
     res_value: TNumber | Tensor
-    query: dict[Parameter, TNumber | Tensor] = {}
+    query: dict[Parameter, TNumber | Tensor] = dict()
+    values = values or dict()
     if isinstance(expr, Array):
         return Tensor(expr.tolist())
     else:

qadence/transpile/__init__.py

@@ -12,6 +12,7 @@ from .circuit import fill_identities
 from .digitalize import digitalize
 from .flatten import flatten
 from .invert import invert_endianness, reassign
+from .noise import set_noise
 from .transpile import blockfn_to_circfn, transpile
 
-__all__ = ["set_trainable", "invert_endianness"]
+__all__ = ["set_trainable", "invert_endianness", "set_noise"]

qadence/transpile/noise.py

@@ -0,0 +1,53 @@
+from __future__ import annotations
+
+from qadence.backend import ConvertedCircuit
+from qadence.blocks.abstract import AbstractBlock
+from qadence.circuit import QuantumCircuit
+from qadence.noise.protocols import NoiseHandler
+
+from .apply_fn import apply_fn_to_blocks
+
+
+def _set_noise(
+    block: AbstractBlock,
+    noise: NoiseHandler | None,
+    target_class: type[AbstractBlock] | None = None,
+) -> AbstractBlock:
+    """Changes the noise protocol of a given block in place."""
+    if target_class is not None:
+        if isinstance(block, target_class):
+            block._noise = noise  # type: ignore [attr-defined]
+    else:
+        block._noise = noise  # type: ignore [attr-defined]
+
+    return block
+
+
+def set_noise(
+    circuit: QuantumCircuit | AbstractBlock | ConvertedCircuit,
+    noise: NoiseHandler | None,
+    target_class: AbstractBlock | None = None,
+) -> QuantumCircuit | AbstractBlock:
+    """
+    Parses a `QuantumCircuit` or `CompositeBlock` to add noise to specific gates.
+
+    If `circuit` is a `ConvertedCircuit`, this is done within `circuit.abstract`.
+
+    Changes the input in place.
+
+    Arguments:
+        circuit: the circuit or block to parse.
+        noise: the NoiseHandler protocol to change to, or `None` to remove the noise.
+        target_class: optional class to selectively add noise to.
+    """
+    input_block: AbstractBlock
+    if isinstance(circuit, ConvertedCircuit):
+        input_block = circuit.abstract.block
+    elif isinstance(circuit, QuantumCircuit):
+        input_block = circuit.block
+    else:
+        input_block = circuit
+
+    apply_fn_to_blocks(input_block, _set_noise, noise, target_class)
+
+    return circuit

qadence/types.py

@@ -1,6 +1,7 @@
 from __future__ import annotations
 
 import importlib
+from dataclasses import dataclass
 from enum import Enum
 from typing import Callable, Iterable, Tuple, Union
 
@@ -8,7 +9,9 @@ import numpy as np
 import sympy
 from matplotlib.figure import Figure
 from numpy.typing import ArrayLike
-from pyqtorch.utils import SolverType
+from pyqtorch.noise import NoiseType as DigitalNoise
+from pyqtorch.noise.readout import WhiteNoise
+from pyqtorch.utils import DropoutMode, SolverType
 from torch import Tensor, pi
 from torch.nn import Module
 
@@ -26,6 +29,7 @@ TArray = Union[Iterable, Tensor, np.ndarray]
 TGenerator = Union[Tensor, sympy.Array, sympy.Basic]
 """Union of torch tensors and numpy arrays."""
 
+
 PI = pi
 
 # Modules to be automatically added to the qadence namespace
@@ -51,6 +55,9 @@ __all__ = [
     "SerializationFormat",
     "PI",
     "SolverType",
+    "DropoutMode",
+    "NoiseProtocol",
+    "WhiteNoise",
 ]  # type: ignore
 
 
@@ -221,8 +228,6 @@ class _BackendName(StrEnum):
 
     PYQTORCH = "pyqtorch"
     """The Pyqtorch backend."""
-    BRAKET = "braket"
-    """The Braket backend."""
     PULSER = "pulser"
     """The Pulser backend."""
     HORQRUX = "horqrux"
@@ -457,3 +462,34 @@ class ExperimentTrackingTool(StrEnum):
 
 
 LoggablePlotFunction = Callable[[Module, int], tuple[str, Figure]]
+
+
+class AnalogNoise(StrEnum):
+    """Type of noise protocol."""
+
+    DEPOLARIZING = "Depolarizing"
+    DEPHASING = "Dephasing"
+
+
+class ReadoutNoise(StrEnum):
+    """Type of readout protocol."""
+
+    INDEPENDENT = "Independent Readout"
+    """Simple readout protocols where each qubit is corrupted independently."""
+    CORRELATED = "Correlated Readout"
+    """Using a confusion matrix (2**n, 2**n) for corrupting bitstrings values."""
+
+
+@dataclass
+class NoiseProtocol:
+    """Type of noise protocol."""
+
+    ANALOG = AnalogNoise
+    """Noise applied in analog blocks."""
+    READOUT = ReadoutNoise
+    """Noise applied on outputs of quantum programs."""
+    DIGITAL = DigitalNoise
+    """Noise applied to digital blocks."""
+
+
+NoiseEnum = Union[DigitalNoise, AnalogNoise, ReadoutNoise]

{qadence-1.7.8.dist-info → qadence-1.9.0.dist-info}/METADATA

@@ -1,10 +1,9 @@
 Metadata-Version: 2.3
 Name: qadence
-Version: 1.7.8
+Version: 1.9.0
 Summary: Pasqal interface for circuit-based quantum computing SDKs
 Author-email: Aleksander Wennersteen <aleksander.wennersteen@pasqal.com>, Gert-Jan Both <gert-jan.both@pasqal.com>, Niklas Heim <niklas.heim@pasqal.com>, Mario Dagrada <mario.dagrada@pasqal.com>, Vincent Elfving <vincent.elfving@pasqal.com>, Dominik Seitz <dominik.seitz@pasqal.com>, Roland Guichard <roland.guichard@pasqal.com>, "Joao P. Moutinho" <joao.moutinho@pasqal.com>, Vytautas Abramavicius <vytautas.abramavicius@pasqal.com>, Gergana Velikova <gergana.velikova@pasqal.com>, Eduardo Maschio <eduardo.maschio@pasqal.com>, Smit Chaudhary <smit.chaudhary@pasqal.com>, Ignacio Fernández Graña <ignacio.fernandez-grana@pasqal.com>, Charles Moussa <charles.moussa@pasqal.com>, Giorgio Tosti Balducci <giorgio.tosti-balducci@pasqal.com>, Daniele Cucurachi <daniele.cucurachi@pasqal.com>
 License: Apache 2.0
-License-File: LICENSE
 Classifier: License :: OSI Approved :: Apache Software License
 Classifier: Programming Language :: Python
 Classifier: Programming Language :: Python :: 3
@@ -22,7 +21,7 @@ Requires-Dist: matplotlib
 Requires-Dist: nevergrad
 Requires-Dist: numpy
 Requires-Dist: openfermion
-Requires-Dist: pyqtorch==1.4.9
+Requires-Dist: pyqtorch==1.6.0
 Requires-Dist: pyyaml
 Requires-Dist: rich
 Requires-Dist: scipy
@@ -31,14 +30,11 @@ Requires-Dist: sympytorch>=0.1.2
 Requires-Dist: tensorboard>=2.12.0
 Requires-Dist: torch
 Provides-Extra: all
-Requires-Dist: braket; extra == 'all'
 Requires-Dist: libs; extra == 'all'
 Requires-Dist: mlflow; extra == 'all'
 Requires-Dist: protocols; extra == 'all'
 Requires-Dist: pulser; extra == 'all'
 Requires-Dist: visualization; extra == 'all'
-Provides-Extra: braket
-Requires-Dist: amazon-braket-sdk<1.71.2; extra == 'braket'
 Provides-Extra: dlprof
 Requires-Dist: nvidia-dlprof[pytorch]; extra == 'dlprof'
 Requires-Dist: nvidia-pyindex; extra == 'dlprof'
@@ -57,9 +53,9 @@ Requires-Dist: mlflow; extra == 'mlflow'
 Provides-Extra: protocols
 Requires-Dist: qadence-protocols; extra == 'protocols'
 Provides-Extra: pulser
-Requires-Dist: pasqal-cloud==0.12.2; extra == 'pulser'
-Requires-Dist: pulser-core==0.19.0; extra == 'pulser'
-Requires-Dist: pulser-simulation==0.19.0; extra == 'pulser'
+Requires-Dist: pasqal-cloud==0.12.5; extra == 'pulser'
+Requires-Dist: pulser-core==1.1.1; extra == 'pulser'
+Requires-Dist: pulser-simulation==1.1.1; extra == 'pulser'
 Provides-Extra: visualization
 Requires-Dist: graphviz; extra == 'visualization'
 Description-Content-Type: text/markdown