Merge branch 'publicrelease_25.11.1' into 'main'

sync with internal repo1 (commit 4bd6388bd)

See merge request cuda-hpc-libraries/cuquantum-sdk/cuquantum-public!43

Author: haidarazzam

python/README.md

@@ -53,7 +53,7 @@ Runtime dependencies of the cuQuantum Python package include:
 * Driver: Linux (525.60.13+ for CUDA 12, 580.65.06+ for CUDA 13)
 * CUDA Toolkit 12.x or 13.x
 * cuStateVec 1.11.0+
-* cuTensorNet 2.10.0+
+* cuTensorNet 2.10.1+
 * cuDensityMat >=0.3.2, <0.4.0
 * cuPauliProp  >=0.1.0, <0.2.0
 * cuStabilizer >=0.1.0, <0.2.0

python/cuquantum/_version.py

@@ -5,4 +5,4 @@
 # Note: cuQuantum Python follows the cuQuantum SDK version, which is now
 # switched to YY.MM and is different from individual libraries' (semantic)
 # versioning scheme.
-__version__ = '25.11.0'
+__version__ = '25.11.1'

python/setup.py

@@ -37,7 +37,7 @@
     'nvmath-python==0.6.0', # strict version before nvmath.internal module is stable
     # 'torch', # <-- PyTorch is optional; also, the PyPI version does not support GPU...
     f'custatevec-cu{utils.cuda_major_ver}~=1.11',  # ">=1.11.0,<2"
-    f'cutensornet-cu{utils.cuda_major_ver}~=2.10',  # ">=2.10.0,<3"
+    f'cutensornet-cu{utils.cuda_major_ver}>=2.10.1, <3',  # ">=2.10.0,<3"
     f'cudensitymat-cu{utils.cuda_major_ver}>=0.3.2, <0.4', # ">=0.3.2,<0.4.0"
     f'cupauliprop-cu{utils.cuda_major_ver}>=0.1.0, <0.2', # ">=0.1.0,<0.2.0"
     f'custabilizer-cu{utils.cuda_major_ver}>=0.1.0, <0.2', # ">=0.1.0,<0.2.0"

python/tests/cuquantum_tests/tensornet/experimental/_internal/mps_utils.py

@@ -213,6 +213,7 @@ def _swap(self, i, direction, exact=False):
         else:
             svd_options = self.svd_options.copy()
             svd_options['max_extent'] = min(self.max_extents[a], mid_extent)
+
         self[a], s, self[b] = tensor_decompose(decompose_expr, tmp, method='svd', **svd_options)
         if s is not None:
             self.gauges[b] = self.backend.asarray(s)
@@ -395,19 +396,18 @@ def _apply_gate_2q(self, i, j, operand):
         else:
             # insert swap gates recursively
             swaps = []
-            exact = self.is_exact_svd
             while (j != i+1):
-                self._swap(i, 'right', exact)
+                self._swap(i, 'right', False)
                 swaps.append([i, 'right'])
                 i += 1
                 if (j == i+1):
                     break
-                self._swap(j, 'left', exact)
+                self._swap(j, 'left', False)
                 swaps.append([j, 'left'])
                 j -= 1
             self._apply_gate_2q(i, j, operand)
             for x, direction in swaps[::-1]:
-                self._swap(x, direction=direction, exact=exact)
+                self._swap(x, direction=direction)
 
     def apply_gate(self, qudits, operand):
         gauge_option = self.gauge_option             
@@ -570,7 +570,6 @@ def from_factory(cls, factory, **kwargs):
             mps_tensors = None
         qudit_dims = factory.state_dims
         mps = cls(qudits, factory.backend, qudit_dims=qudit_dims, mps_tensors=mps_tensors, dtype=factory.dtype, **kwargs)
-        
         for op, modes, gate_info in factory.sequence:
             if gate_info is None:
                 if isinstance(op, (list, tuple)):

python/tests/cuquantum_tests/tensornet/experimental/_internal/state_matrix.py

@@ -155,7 +155,7 @@ def create_state_factory(*args, **kwargs):
 
 generic_states_L1 = [
     create_state_factory(5, "float32", "SDDS", next(rng_iterator)),
-    create_state_factory((2, 3, 4, 3, 2), "complex64", "SDDDS", next(rng_iterator), adjacent_double_layer=False),
+    create_state_factory((2, 3, 4, 3, 2), "complex64", "SDDS", next(rng_iterator)),
     create_state_factory(5, "complex64", "SDMS", next(rng_iterator), adjacent_double_layer=False, mpo_bond_dim=3),
     create_state_factory(4, "float64", "SDMDS", next(rng_iterator), mpo_bond_dim=2, mpo_num_sites=3, mpo_geometry="random-ordered"),
     create_state_factory((2, 5, 3, 2), "complex128", "SDMDS", next(rng_iterator), mpo_bond_dim=2, mpo_geometry="random"),

python/tests/cuquantum_tests/tensornet/experimental/test_network_state.py

@@ -341,6 +341,27 @@ def test_negative_complex_circuit(self, circuit, backend, dtype):
             with NetworkState.from_circuit(circuit, backend=backend, dtype=dtype) as state:
                 pass
             assert "imaginary part" in str(e.value)
+    
+    def test_mps_with_fixed_bond_truncation(self):
+        num_qubits = 6
+        num_double_layers = 2
+
+        config = MPSConfig(max_extent=3)
+        state = NetworkState((2, ) * num_qubits, dtype='float64', config=config)
+
+        for i in range(num_qubits):
+            state.apply_tensor_operator((i,), np.random.random([2, 2]))
+
+        for _ in range(num_double_layers):
+            for i in range(2):
+                for j in range(i, num_qubits-1, 2):
+                    state.apply_tensor_operator((j, j+1), np.random.random([2, 2, 2, 2]))
+
+        state.apply_tensor_operator((1, 3), np.random.random([2, 2, 2, 2]))
+
+        with state:
+            mps = state.compute_output_state()
+            assert mps is not None
 
 
 @pytest.fixture(params=CircuitStateMatrix.L1(), scope="class")

samples/custabilizer/custabilizer_demo.cpp

@@ -2,7 +2,7 @@
  * Copyright (c) 2025, NVIDIA CORPORATION & AFFILIATES.
  *
  * SPDX-License-Identifier: BSD-3-Clause
-*/
+ */
 
 #include <iostream>
 #include <vector>