tensorcircuit-nightly 1.3.0.dev20250908__py3-none-any.whl → 1.3.0.dev20250910__py3-none-any.whl

tensorcircuit/__init__.py

@@ -1,4 +1,4 @@
-__version__ = "1.3.0.dev20250908"
+__version__ = "1.3.0.dev20250910"
 __author__ = "TensorCircuit Authors"
 __creator__ = "refraction-ray"
 
@@ -23,8 +23,10 @@ from .cons import (
     runtime_contractor,
 )  # prerun of set hooks
 from . import gates
+from . import quditgates
 from . import basecircuit
 from .gates import Gate
+from .quditcircuit import QuditCircuit
 from .circuit import Circuit, expectation
 from .mpscircuit import MPSCircuit
 from .densitymatrix import DMCircuit as DMCircuit_reference

tensorcircuit/applications/layers.py

@@ -35,7 +35,7 @@ Symbol = Any  # sympy.Symbol
 
 def _resolve(symbol: Union[Symbol, Tensor], i: int = 0) -> Tensor:
     """
-    Make sure the layer is compatible with both multi-param and single param requirements。
+    Make sure the layer is compatible with both multi-param and single param requirements
 
     What could be the input: list/tuple of sympy.symbol, tf.tensor with 1D or 0D shape
     """

tensorcircuit/backends/abstract_backend.py

@@ -1039,8 +1039,8 @@ class ExtendedBackend:
         :type a: Tensor
         :param v: value to inserted
         :type v: Tensor
-        :param side:  If ‘left’, the index of the first suitable location found is given.
-            If ‘right’, return the last such index.
+        :param side:  If `left`, the index of the first suitable location found is given.
+            If `right`, return the last such index.
             If there is no suitable index, return either 0 or N (where N is the length of a),
             defaults to "left"
         :type side: str, optional

tensorcircuit/basecircuit.py

@@ -3,7 +3,7 @@ Quantum circuit: common methods for all circuit classes as MixIn
 
 Note:
   - Supports qubit (d = 2) and qudit (d >= 2) systems.
-  - For string-encoded samples/counts when d <= 36, digits use base-d characters 0–9A–Z (A = 10, …, Z = 35).
+  - For string-encoded samples/counts when d <= 36, digits use base-d characters 0-9A-Z (A = 10, ..., Z = 35).
 """
 
 # pylint: disable=invalid-name
@@ -362,7 +362,7 @@ class BaseCircuit(AbstractCircuit):
 
     def perfect_sampling(self, status: Optional[Tensor] = None) -> Tuple[str, float]:
         """
-        Sampling base-d strings (0–9A–Z when d <= 36) from the circuit output based on quantum amplitudes.
+        Sampling base-d strings (0-9A-Z when d <= 36) from the circuit output based on quantum amplitudes.
         Reference: arXiv:1201.3974.
 
         :param status: external randomness, with shape [nqubits], defaults to None
@@ -480,28 +480,12 @@ class BaseCircuit(AbstractCircuit):
 
     def amplitude_before(self, l: Union[str, Tensor]) -> List[Gate]:
         r"""
-        Returns the tensornetwor nodes for the amplitude of the circuit given a computational-basis label ``l``.
-        For a state simulator, it computes :math:`\langle l \vert \psi\rangle`;
-        for a density-matrix simulator, it computes :math:`\mathrm{Tr}(\rho \vert l\rangle\langle l\vert)`.
+        Returns the tensornetwor nodes for the amplitude of the circuit given the bitstring l.
+        For state simulator, it computes :math:`\langle l\vert \psi\rangle`,
+        for density matrix simulator, it computes :math:`Tr(\rho \vert l\rangle \langle 1\vert)`
         Note how these two are different up to a square operation.
 
-        :Example:
-
-        >>> c = tc.Circuit(2)
-        >>> c.X(0)
-        >>> c.amplitude("10")  # d=2, per-qubit digits
-        array(1.+0.j, dtype=complex64)
-        >>> c.CNOT(0, 1)
-        >>> c.amplitude("11")
-        array(1.+0.j, dtype=complex64)
-
-        For qudits (d>2, d<=36):
-        >>> c = tc.Circuit(3, dim=12)
-        >>> c.amplitude("0A2")  # base-12 string, A stands for 10
-
-        :param l: Basis label.
-            - If a string: it must be a base-d string of length ``nqubits``, using 0–9A–Z (A=10,…,Z=35) when ``d<=36``.
-            - If a tensor/array/list: it should contain per-site integers in ``[0, d-1]`` with length ``nqubits``.
+        :param l: The bitstring of 0 and 1s.
         :type l: Union[str, Tensor]
         :return: The tensornetwork nodes for the amplitude of the circuit.
         :rtype: List[Gate]
@@ -620,7 +604,7 @@ class BaseCircuit(AbstractCircuit):
                 "count_tuple": # (np.array([0]), np.array([2]))
 
                 "count_dict_bin": # {"00": 2, "01": 0, "10": 0, "11": 0}
-                    for cases d\in [11, 36], use 0–9A–Z digits (e.g., 'A' -> 10, …, 'Z' -> 35);
+                    for cases d\in [11, 36], use 0-9A-Z digits (e.g., 'A' -> 10, ..., 'Z' -> 35);
 
                 "count_dict_int": # {0: 2, 1: 0, 2: 0, 3: 0}
 

tensorcircuit/circuit.py

@@ -1,7 +1,7 @@
 """
 Quantum circuit: the state simulator.
 Supports qubit (dim=2) and qudit (3 <= dim <= 36) systems.
-For string-encoded samples/counts, digits use 0–9A–Z where A=10, …, Z=35.
+For string-encoded samples/counts, digits use 0-9A-Z where A=10, ..., Z=35.
 """
 
 # pylint: disable=invalid-name
@@ -768,7 +768,7 @@ class Circuit(BaseCircuit):
         Take measurement on the given quantum lines by ``index``.
 
         Return format:
-        - For d <= 36, the sample is a base-d string using 0–9A–Z (A=10,…).
+        - For d <= 36, the sample is a base-d string using 0-9A-Z (A=10,...).
 
         :Example:
 
@@ -874,7 +874,7 @@ class Circuit(BaseCircuit):
         :param nmc: repetition time for Monte Carlo sampling for noisfy calculation, defaults to 1000
         :type nmc: int, optional
         :param status: external randomness given by tensor uniformly from [0, 1], defaults to None,
-            used for noisfy circuit sampling
+            used for noisy circuit sampling
         :type status: Optional[Tensor], optional
         :raises ValueError: "Cannot measure two operators in one index"
         :return: Tensor with one element

tensorcircuit/densitymatrix.py

@@ -179,9 +179,9 @@ class DMCircuit(BaseCircuit):
 
     @staticmethod
     def check_kraus(kraus: Sequence[Gate]) -> bool:
-        """
+        r"""
         Check if Kraus operators satisfy the completeness relation:
-        sum_i K_i^† K_i = I
+        :math:`\sum_i K_i^\dagger K_i = I`
 
         :param kraus: Sequence of Kraus operators
         :type kraus: Sequence[Gate]

tensorcircuit/quantum.py

@@ -79,7 +79,7 @@ def _decode_basis_label(label: str, n: int, dim: int) -> List[int]:
     """
     Decode a string basis label into a list of integer digits.
 
-    The label is interpreted in base-``dim`` using characters ``0–9A–Z``.
+    The label is interpreted in base-``dim`` using characters ``0-9A-Z``.
     Only dimensions up to 36 are supported.
 
     :param label: basis label string, e.g. "010" or "A9F"
@@ -97,7 +97,7 @@ def _decode_basis_label(label: str, n: int, dim: int) -> List[int]:
     """
     if dim > 36:
         raise NotImplementedError(
-            f"String basis label supports d<=36 (0–9A–Z). Got dim={dim}. "
+            f"String basis label supports d<=36 (0-9A-Z). Got dim={dim}. "
             "Use an integer array/tensor of length n instead."
         )
     s = label.upper()
@@ -107,7 +107,7 @@ def _decode_basis_label(label: str, n: int, dim: int) -> List[int]:
     for ch in s:
         if ch not in _ALPHABET:
             raise ValueError(
-                f"Invalid character '{ch}' in basis label (allowed 0–9A–Z)."
+                f"Invalid character '{ch}' in basis label (allowed 0-9A-Z)."
             )
         v = _ALPHABET.index(ch)
         if v >= dim:
@@ -748,10 +748,10 @@ class QuOperator:
         return self.__mul__(other)
 
     def tensor_product(self, other: "QuOperator") -> "QuOperator":
-        """
+        r"""
         Tensor product with another operator.
         Given two operators `A` and `B`, produces a new operator `AB` representing
-        :math:`A ⊗ B`. The `out_edges` (`in_edges`) of `AB` is simply the
+        :math:`A \otimes B`. The `out_edges` (`in_edges`) of `AB` is simply the
         concatenation of the `out_edges` (`in_edges`) of `A.copy()` with that of
         `B.copy()`:
         `new_out_edges = [*out_edges_A_copy, *out_edges_B_copy]`
@@ -2403,13 +2403,13 @@ def free_energy(
 
 def renyi_entropy(rho: Union[Tensor, QuOperator], k: int = 2) -> Tensor:
     """
-    Compute the Rényi entropy of order :math:`k` by given density matrix.
+    Compute the Renyi entropy of order :math:`k` by given density matrix.
 
     :param rho: The density matrix in form of Tensor or QuOperator.
     :type rho: Union[Tensor, QuOperator]
-    :param k: The order of Rényi entropy, default is 2.
+    :param k: The order of Renyi entropy, default is 2.
     :type k: int, optional
-    :return: The :math:`k` th order of Rényi entropy.
+    :return: The :math:`k` th order of Renyi entropy.
     :rtype: Tensor
     """
     s = 1 / (1 - k) * backend.real(backend.log(trace_product(*[rho for _ in range(k)])))
@@ -2423,7 +2423,7 @@ def renyi_free_energy(
     k: int = 2,
 ) -> Tensor:
     """
-    Compute the Rényi free energy of the corresponding density matrix and Hamiltonian.
+    Compute the Renyi free energy of the corresponding density matrix and Hamiltonian.
 
     :Example:
 
@@ -2440,9 +2440,9 @@ def renyi_free_energy(
     :type h: Union[Tensor, QuOperator]
     :param beta: Constant for the optimization, default is 1.
     :type beta: float, optional
-    :param k: The order of Rényi entropy, default is 2.
+    :param k: The order of Renyi entropy, default is 2.
     :type k: int, optional
-    :return: The :math:`k` th order of Rényi entropy.
+    :return: The :math:`k` th order of Renyi entropy.
     :rtype: Tensor
     """
     energy = backend.real(trace_product(rho, h))